Skip to main content

Metagenomic profiling of gut microbiota in Fall Armyworm (Spodoptera frugiperda) larvae fed on different host plants

Abstract

Background

The fall armyworm (FAW, Spodoptera frugiperda) is a polyphagous pest known for causing significant crop damage. The gut microbiota plays a pivotal role in influencing the biology, physiology and adaptation of the host. However, understanding of the taxonomic composition and functional characteristics of the gut microbiota in FAW larvae fed on different host plants remains limited.

Methods

This study utilized metagenomic sequencing to explore the structure, function and antibiotic resistance genes (ARGs) of the gut microbiota in FAW larvae transferred from an artificial diet to four distinct host plants: maize, sorghum, tomato and pepper.

Results

The results demonstrated significant variations in gut microbiota structure among FAW larvae fed on different host plants. Firmicutes emerged as the dominant phylum, with Enterococcaceae as the dominant family and Enterococcus as the prominent genus. Notably, Enterococcus casseliflavus was frequently observed in the gut microbiota of FAW larvae across host plants. Metabolism pathways, particularly those related to carbohydrate and amino acid metabolism, played a crucial role in the adaptation of the FAW gut microbiota to different host plants. KEGG orthologs associated with the regulation of the peptide/nickel transport system permease protein in sorghum-fed larvae and the 6-phospho-β-glucosidase gene linked to glycolysis/gluconeogenesis as well as starch and sucrose metabolism in pepper-fed larvae were identified. Moreover, the study identified the top 20 ARGs in the gut microbiota of FAW larvae fed on different host plants, with the maize-fed group exhibiting the highest abundance of vanRC.

Conclusions

Our metagenomic sequencing study reveals significant variations in the gut microbiota composition and function of FAW larvae across diverse host plants. These findings underscore the intricate co-evolutionary relationship between hosts and their gut microbiota, suggesting that host transfer profoundly influences the gut microbiota and, consequently, the adaptability and pest management strategies for FAW.

Peer Review reports

Background

The fall armyworm (FAW), Spodoptera frugiperda (J.E. Smith) (Lepidoptera: Noctuidae), poses a formidable threat as an invasive and highly polyphagous agricultural pest that is native to tropical and subtropical regions of the Americas [1,2,3]. Its rapid reproduction in local maize crops has led to population outbreaks, posing severe challenges to crop production [4,5,6].

The microbial community residing within the gut of insects encompasses a diverse array of organisms, including viruses, archaea, bacteria, fungi, and protozoa [7,8,9]. As part of the co-evolutionary dynamics between microbes and their hosts, these entities engage in complex interactions that significantly impact host metabolism and physiological functions [10]. Studies have elucidated the profound influence of gut microbiota on various aspects of insect biology, such as longevity, development, fecundity and behavior regulation. They play critical roles in nutrient provisioning, digestion facilitation, as well as immune defense and detoxification mechanisms [11,12,13,14].

The gut microbiota of insects is influenced by the host plants they feed on; this impact is a critical aspect of understanding insect–microbe interactions and their ecological and evolutionary implications. The feeding habits of insects significantly impact the composition, structure, diversity and function of their gut bacterial communities. Zhu et al. observed distinct bacterial diversity and community structures in Chilo suppressalis corresponding to different diets [15]. Similarly, Yang et al. established a significant correlation between the gut microbiota of Cnaphalocrocis medinalis and its respective host plant [16]. Meng et al. highlighted the influence of host diet on the bacterial community structure of Diaphorina citri [17]. Wang et al. demonstrated significant variations in the gut bacterial community of Locusta migratoria manilensis in response to different host plants [18]. Furthermore, Yuan et al. conducted a study on the gut community structure of Grapholita molesta, revealing noteworthy shifts when transferring from an artificial diet to various host fruits [19]. Collectively, these findings highlight the critical importance of understanding the impact of host plants on insect gut microbiota, given its profound ecological and evolutionary implications.

Several studies have demonstrated the significant influence of host plants on the gut microbiota of FAW larvae. Ma et al. observed that FAW larvae fed on hairy vetch displayed reduced gut microbiota diversity compared to those fed on maize, with specific microbial taxa exhibiting a notable correlation with larval performance [20]. Similarly, Wang et al. found that the richness and diversity of gut bacteria in fifth-instar FAW larvae peaked when fed on corn, whereas sixth-instar larvae showed higher measures when fed on sorghum, highland barley, and citrus [21]. De Oliveira and Cônsoli revealed genus-level differences in FAW larval gut microbiota when feeding on different corn and rice strains, with notable phylum-level differences observed when larvae fed on cotton plants [22]. Furthermore, Ugwu et al. identified distinct patterns in the structure and community of gut bacteria in FAW larvae fed on sugarcane, maize, and onion [23]. Lü et al. compared the gut microbial community structure and diversity in FAW larvae feeding on various hosts, including corn with and without seed coating agent, wild oat, oilseed rape, pepper, and an artificial diet [24]. Their results highlighted that the gut microbiota of FAW larvae is not only influenced by host species but also by various host treatments. These findings underscore the significant impact of host plant selection on the gut microbiota of FAW larvae, emphasizing the importance of considering plant species in understanding insect–microbe interactions. They offer valuable insights into the intricate relationship between host plants, gut microbiota, and FAW development, underscoring the necessity for further research in this area.

Previous studies have extensively investigated the gut microbiota of FAW larvae using both culture-dependent methods [25,26,27,28,29,30] and amplicon-oriented DNA sequencing [29, 31,32,33,34,35]. Metagenomic sequencing offers distinct advantages over traditional culture techniques, as it can detect microbial species that are uncultivable. Moreover, compared to 16S rRNA amplification sequencing, metagenomic sequencing allows for the identification of microbial species at various taxonomic levels, facilitating a more precise analysis of microbial community structure and function [36].

Metagenomic sequencing is a key method for studying the impact of pesticides on soil microbes. It allows researchers to understand the effects on microbial diversity and function, The complex effects of pesticide mixtures on soil microbial communities and functions provide a new perspective for understanding the ecological risks associated with pesticides [37]. This technique has also been effective in water treatment, analyzing microbial dynamics, tracking antibiotic resistance genes, and identifying human pathogens, aiding in strategies to combat antibiotic resistance [38]. When discussing the role of metagenomics in plant disease management, Jyoti Taunk and Umesh Goutam (2021) detailed in Chapter 9 of the book "Role of Metagenomics in Plant Disease Management" how metagenomic techniques can be used to identify novel biological control agents. Through high-throughput sequencing and bioinformatics analysis, they successfully identified a variety of microorganisms and functional genes that have significant effects on the prevention and control of plant diseases, providing an important reference for the innovation of biological control strategies [39]. Additionally, metagenomic has provided insights into insect gut microbiomes, Herbivorous insects like the diamondback moth, Plutella xylostella, have evolved mechanisms to counteract plant defenses, and their gut microbiota, dominated by a few key bacterial species, plays a crucial role in aiding digestion, detoxifying plant compounds, and synthesizing essential amino acids, presenting possible targets for innovative pest control strategies [40]. The gut microbiota of the rice pest Cnaphalocrocis medinalis plays a vital role in its nutritional intake and development, with the core community remaining stable across different diets, and the study's findings on the symbiotic relationship between the insect and its microbes offer new perspectives for pest management strategies [41]. As a powerful tool, metagenomics has demonstrated its unique value across various domains. From soil microbial communities to insect gut microbiota, and from plant diseases to antibiotic resistance gene research, metagenomics provides a comprehensive and in-depth perspective. This aids in better understanding and managing the complex ecological and health issues related to microorganisms.

In this study, we employed metagenomic sequencing to explore the gut microbiota communities of FAW larvae reared in the laboratory using two Poaceae crops (maize and sorghum), two Solanaceous plants (pepper and tomato), and an artificial diet. The findings of the current research revealed that the gut microbial structure of FAW larvae can be influenced by the host plants, subsequently impacting the functional genes associated with gut microbiota and antibiotic resistance. These findings suggest that the gut microbiota of FAW may play a pivotal role in host adaptation.

Methods

Insects rearing

FAW individuals were collected from a maize field near the Guizhou Academy of Agricultural Sciences in June 2019 and were fed an artificial diet, as described by Di et al. [42] to establish a multigenerational laboratory population at the Institute of Entomology, Guizhou University, located in Huaxi, Guiyang, Guizhou, China. For our study, we collected eggs from the FAW population that was fed the artificial diet. Newly hatched larvae were individually placed in 6-well culture plates. These larvae were then provided with fresh leaves from four different host plants: pepper, tomato, maize, sorghum, and the artificial diet on a daily basis. The rearing conditions were maintained in an incubator at 25°C ± 1°C and 60% ± 5% relative humidity, under a photoperiod of 16-h light and 8-h dark. Adult FAWs were provided with a 10% (w/v) honey solution diluted in water.

Host plants and artificial diet formula

The host plants used in this study included pepper, tomato, maize, and sorghum (Table 1). These plants were cultivated in pots filled with soil sourced from a vegetable experimental base and were maintained in a greenhouse at 25°C ± 2°C in Guizhou University, Guizhou, China. Tender leaves from maize, pepper, tomato, and sorghum plants were harvested at the specific 3–4-true-leaf stages for all subsequent experiments. These leaves were then provided daily to feed the newly hatched FAW larvae until they reached the fifth instar stage.

Table 1 Host plants used for feeding FAW larvae

According to Di et al. [42] an improved artificial diet was used, which consisted of 400 mL of distilled water, 40 g soybean powder, 40 g wheat powder, 16 g yeast powder, 8 g casein, 3.2 g ascorbic acid, 8 g agar powder, 0.4 g choline chloride, 0.8 g sorbic acid, 0.08 g inositol, and 0.08 g cholesterol. In this experiment, two factors were considered that constituted the treatment: host plants and artificial diet.

Sample collection

FAW larvae were fed on various host plants until they reached the fifth instar before being collected for experimentation [29]. Upon collection, the larvae were washed with sterile water and surface-sterilized in 75% ethanol for 90 s, followed by three rinses in sterile water for 30 s each. The entire gut was aseptically dissected using sterile forceps in a sterile petri dish (placed on ice) with a diameter of 90 mmwithin an ultraclean bench. Each of the five treatments (pepper, tomato, maize, sorghum and artificial diet) consisted of 60 individual guts, which were then stored in sterile 1.5-mL centrifuge tubes [43, 44], with three replicates for each treatment. In total, 15 whole gut samples were collected. All gut samples were immediately stored in liquid nitrogen after collection and preserved at − 80°C until DNA extraction.

Sample test

The cetyltrimethylammonium bromide (CTAB) method was used for total metagenomic DNA extraction from each sample according to the protocol. DNA concentration was determined using a NanoDrop 2000 spectrophotometer, and the purity and integrity of the samples were verified via agarose gel electrophoresis prior to sequencing. Sequencing of the total DNA was performed by Wekemo Tech Group Co., Ltd. in Shenzhen, Guangdong, China.

Library construction

A DNA library was constructed using the NEBNext® UltraTM DNA Library Prep Kit according to the manufacturer’s instructions (Illumina, San Diego, CA, USA). Genomic DNA samples were fragmented by sonication to a size of 350 bp. The DNA fragments were subjected to end polishing, A-tailing, and ligation with full-length adapters for Illumina sequencing. Subsequently, polymerase chain reaction (PCR) was performed. The PCR products were purified using the AMPure XP system (Beverly, CA, USA). The quality of the libraries was assessed using the Agilent 5400 system (Agilent, Santa Clara, CA, USA), and quantification was performed using QPCR, achieving a concentration of 1.5 nM.

Sequencing

We opted for a paired-end (PE) library with an insert size of 350 bp. Subsequently, all samples were subjected to sequencing using the PE150 (2 × 150) protocol on the Illumina Novaseq 6000 platform (Illumina, San Diego, CA, USA).

Bioinformatics analysis

Data quality control and de-host sequence

Metagenomic sequencing was conducted using the Illumina Novaseq high-throughput sequencing platform to capture raw data encompassing bacteria, fungi, and viruses across 15 samples. Quality control and sequence trimming were performed using the KneadData toolkit (version 0.7.4) to ensure the reliability of the data. Illumina adapters and low-quality sequences were removed from raw reads using Trimmomatic (version 0.39) [45]. Furthermore, to eliminate any potential contamination from the host S. frugiperda genome, KneadData integrated with the Bowtie2 tool (version 2.3.4.1) was used [46]. Subsequently, the efficacy and impact of the quality control measures were assessed using FastQC to generate clean and reliable data.

Taxonomy annotation

The Kraken2 program was utilized in conjunction with a self-constructed microbial database [47]. This database comprised sequences from bacteria, fungi, archaea, and viruses, which were filtered from the NT nucleic acid database and the RefSeq whole genome database of NCBI. These resources were employed to classify the species present in all samples. Furthermore, the Bayesian Re-estimation of Abundance after Classification with Kraken (BRACKen) was used to estimate the actual abundance of species in the samples [48].

Identification of taxa with differential abundance

The phylogenetic diversity among 15 samples was examined based on alpha diversity, calculated using the vegan (version 2.6.0). This provided insights into the richness and evenness of each sample [49]. For the analysis of beta diversity, principal coordinate analysis (PCoA) was employed based on Bray–Curtis metrics to compare dissimilarities among different groups [50]. The R package Venn diagram (ggplot2 version 3.3.5.9000) was utilized to visualize common and unique taxa at the species level among the five groups. To detect significantly differentially enriched taxa among these groups, linear discriminant analysis (LDA) effect size (LEfSe) was used, LEfSe is a statistical method that combines the Kruskal–Wallis rank sum test with LDA to identify and quantify the differentially abundant microbial taxa and functional genes, offering a way to reveal significant and biologically relevant distinctions within microbial communities, and set thresholds of LDA > 2–4 to differentiate species of high abundance [51].

Common function database annotation of metagenomic DNA sequences

The metagenomic DNA sequences were aligned against the UniProt Reference Clusters 90 (UniRef90) database using the DIAMOND algorithm [52]. Subsequently, a tiered search strategy was implemented in the HMP Unified Metabolic Analysis Network (HUMAnN2) [53]. This strategy was used to quantify the abundance of Kyoto Encyclopedia of Genes and Genomes (KEGG) orthologs (KOs) against the KEGG database [54]. Using LEfSe, significantly differentially abundant KOs were identified among the groups. This process also facilitated the identification of potential functional species associated with different KOs.

Antibiotic resistance gene (ARG) analysis

After quality control, host sequences in each sample were removed using the ARG database. The DIAMOND software was then employed for blasting and annotating clean reads [52]. The Comprehensive Antibiotic Resistance Database (CARDS) algorithm was utilized with specific parameters set to − e 0.001 (e-value < le − 3) and − i 80 (percent identity > 80%) for each sample.

Availability of data and materials

All metagenome sequence reads can be accessed in the Sequence Read Archive (SRA) of NCBI under the BioProject with the accession number PRJNA1123873, available at https://www.ncbi.nlm.nih.gov/sra/PRJNA1123873.

Results

Summary of metagenomic sequencing results

Metagenomic sequencing was employed to analyze 15 samples, with the sequencing performed using the Illumina Novaseq platform. A total of 327,709,326 raw reads were collected, which, after filtering and quality control, were processed to generate 100,370,155 clean reads. These clean reads were distributed across samples, ranging from 3,018,373 to 14,229,565 per sample, with an average of 6,691,344 clean reads per sample (Table 2). Cluster analysis resulted in 1,904 operational taxonomic units (OTUs) at a 97% sequence similarity, encompassing 14 phyla, 40 classes, 92 orders, 182 families, 604 genera, and 1,835 species. The sample rarefaction curve (Fig. 1a) and Shannon index rarefaction curve (Fig. 1b) confirmed that the sequencing depth was sufficient, indicating that increasing the sample volume would not generate additional OTUs.

Table 2 Summary of metagenome sequence statistics for the Illumina Novaseq runs of all samples
Fig. 1
figure 1

The sample rarefaction curve (a) and Shannon index rarefaction curve (b); Cluster heatmap of the 20 abundant genera in the gut microbial community of FAW larvae fed on an artificial diet (AD), pepper (CA), tomato (SL), sorghum (SM) and maize (ZM) (c)

Composition of gut microbiota of FAW fed on different host plants

The taxonomic units detected in this study included archaea (119), bacteria (44,039,175), phages (3,997), plants, and fungi (122,422). The proportions of each category were 0.00%, 99.73%, 0.00%, and 0.27%, respectively (Table 3).

Table 3 The relative abundance of 15 samples at the Kingdom level

At the phylum level, the three dominant phyla were Firmicutes, Proteobacteria, and Actinobacteria, all of which belong to bacteria. These were followed by Microsporidia and Ascomycota, which belong to fungi. Firmicutes was the dominant phylum among the gut microbiota of the majority of the groups, with the notable exception of the sorghum-fed group, especially in the tomato-fed group, accounting for 98.73% ± 0.37%. Proteobacteria was the dominant phylum in the sorghum-fed group, accounting for 59.64% ± 11.16% (Fig. 2a, Table 4).

Fig. 2
figure 2

The ten most abundant gut microbiota of FAW larvae fed on an artificial diet (AD), pepper (CA), tomato (SL), sorghum (SM) and maize (ZM) at the phylum level (c); the family level (d); the genus level (e); and the species level (f)

Table 4 The relative abundance of five dominant phyla in FAW fed on the four host plants and an artificial diet

At the family level, the five most prevalent families were Enterococcaceae, Lactobacillaceae, Enterobacteriaceae, Brucellaceae, and Rhizobiaceae (Fig. 2b). Enterococcaceae was present in all samples, with the highest relative abundance observed in samples collected from the tomato-fed group at 98.64% ± 0.36%. Conversely, the lowest relative abundance of Enterococcaceae was found in samples collected from the sorghum-fed group at 37.51% ± 11.96%. Lactobacillaceae only exhibited high abundance (26.42%% ± 6.76%) when FAW larvae were fed on the artificial diet. Additionally, Enterobacteriaceae (12.61% ± 2.54%), Brucellaceae (21.04% ± 2.13%), and Rhizobiaceae (16.73% ± 8.22%) were more abundant when FAW larvae were feeding on sorghum compared to the other groups (Table 5).

Table 5 The relative abundance of five most prevalent families in FAW fed on the four host plants and an artificial diet

At the genus level, Enterococcus, Lactobacillus, and Ochrobactrum were the most dominant in the gut microbiota across all samples of multiple groups, including the tomato, pepper, maize, and artificial diet-fed groups, followed closely by Klebsiella and Agrobacterium. In the artificial diet-fed group, Lactobacillus was found to be the most abundant genus (Fig. 2c).

Among the five groups observed, the most abundant species were Enterococcus casseliflavus, Enterococcus mundtii, Enterococcus sp. CR-Ec1, Enterococcus gallinarum, and Enterococcus sp. FDAARGOS 375 (Fig. 2d).

To elucidate the changes in the gut microbiota of FAW larvae when fed on different host plants, we selected the top 20 genera based on relative abundances and created a relative abundance cluster heatmap (Fig. 1c). Among the samples, Enterococcus, belonging to the phylum Firmicutes, emerged as the most abundant genus. Ochrobactrum, from the phylum Proteobacteria, was the primary component of the bacterial communities in sorghum-fed FAW. Additionally, the genera Klebsiella (5.64), Agrobacterium (5.53), and Devosia (4.94), all from the phylum Proteobacteria, were also abundant in the gut microbiota of sorghum-fed FAW. Lactobacillus, from the phylum Firmicutes, was the major component of the bacterial communities in FAW larvae fed an artificial diet. Interestingly, Agrobacterium and Rhizobium were exclusively present in the gut microbiota of FAW feeding on Poaceae crops (maize and sorghum). Moreover, within the gut microbiota of FAW larvae feeding on sorghum, the abundance of Agrobacterium and Rhizobium was higher compared to those feeding on maize. However, these two genera were not found in the gut microbiota of FAW fed on an artificial diet or pepper.

Comparison of gut microbiota diversity among FAW fed on different host plants

Alpha diversity

The alpha diversity index was used to identify differences in the gut microbiota of FAW larvae when fed on different host plants. The richness of the gut microbiota was evaluated using the ACE and Chao 1 indices. The FAW larvae feeding on sorghum had the highest ACE and Chao 1 values, indicating a rich gut microbiota. In contrast, the maize-fed FAW larvae had the lowest values for both ACE and Chao 1, suggesting a less rich gut microbiota (Fig. 3a, b). The diversity of the gut microbiota was evaluated using the Simpson and Shannon indices. The sorghum-fed group had the highest values in both Simpson and Shannon, indicating high diversity. Conversely, the pepper-fed group had the lowest Simpson and Shannon diversity, indicating lower diversity (Fig. 3c, 3d).

Fig. 3
figure 3

The alpha diversity index was used to identify differences in the gut microbiota of FAW larvae fed on different host plants. ACE indices (h); Chao 1 indices (i); Simpson indices (j); and Shannon indices (k). The abbreviations in the figure represent different diets, AD, artificial diet-fed; CA, pepper-fed; SL, tomato-fed; SM, sorghum-fed; ZM, maize-fed. Different lowercase letters indicate significant differences (one-way ANOVA, Tukey post-hoc test, P < 0.05) in the mean values

Differentially enriched taxa associated with FAW fed on host plants

The beta diversity, evaluated by the Bray–Curtis distance, clearly delineated the microbiota composition between samples (Fig. 4a). The gut microbiota of FAW larvae fed on sorghum exhibited a significant difference in composition compared to those fed on other host plants. Additionally, the gut microbiota of FAW larvae fed on an artificial diet and sorghum was distinct from that of FAW larvae fed on maize and pepper. The PCoA plot revealed distinct clustering of samples from pepper-fed and maize-fed larvae, although some overlap was observed.

Fig. 4
figure 4

Differential taxa were enriched in the AD, CA, SL, SM and ZM groups. The PCoA plot showed the beta diversity assessed by the Bray–Curtis distance (a); Venn diagram represented common and unique taxa at the species level (b); Taxonomic cladogram from LEfSe showed different taxa enriched in the AD, CA, SL, and SM group (LDA > 4) (c); LEfSe identified significantly differentially abundant taxa (LDA > 4) (d). The abbreviations in the figure represent different diets, AD, artificial diet-fed; CA, pepper-fed; SL, tomato-fed; SM, sorghum-fed; ZM, maize-fed. PCoA, Principal coordinate analysis; LEfSe, linear discriminant analysis effect size; LDA, linear discriminant analysis

A Venn diagram was used to represent the core microbiome, consisting of shared taxa among samples within certain groups. At the species level, 84 common species were shared among the five groups. There were 80, 37, 36, 1417, and 2 unique species in the artificial diet-fed, pepper-fed, tomato-fed, sorghum-fed, and maize-fed groups, respectively (Fig. 4b).

To identify biomarkers with significant variations across different samples, we utilized the LEfSe method to analyze various taxonomic orders, including kingdom, phylum, class, order, family, genus, and species. Our approach involved setting a minimum LDA value of four to screen these categories, thereby establishing a benchmark for statistical significance. Nineteen taxa were clustered within Proteobacteria, with the dominant Alpha-proteobacteria being the most highly discriminating taxon in the sorghum-fed group. The gut microbiota of the tomato-fed group showed significant differences across seven taxa, primarily featuring Firmicutes, Enterococcaceae, Enterococcus, and Lactobacillus. Only one microbe, belonging to the species Enterococcus casseliflavus in the phylum Firmicutes, was found in the larvae that fed on pepper. Eight different taxa showed a significant difference in the gut microbiota of the larvae fed on the artificial diet, with a predominant in Enterococcus mundtii, Lactobacillaceae, and Lactobacillus. The absence of microbiome features in the image regarding the larvae feeding on maize indicates no distinguishable species in the gut microbiome between the larvae fed on maize and the other host plants (Fig. 4c, d).

Functional analysis in gut microbiota of FAW larvae fed on different host plants

To investigate the functionality of the gut microbiota across five groups, we annotated the gene catalogs with metagenomic sequencing data using the KEGG database. The KEGG PATHWAY database, which includes areas such as metabolism, genetic information processing, cellular processes, human diseases, environmental information processing, and organismal systems, served as the foundation of the KEGG database. Among these pathways, those related to metabolism were the most extensively annotated, accounting for 70.36% of the total (Fig. 5a, 6a, Supplementary Table S1).

Fig. 5
figure 5

Functional annotation of the gut microbiome across AD, CA, SL, SM and ZM groups at the KEGG first level pathway analysis (a); KEGG second level pathway analysis (b); KEGG third-level pathway analysis (c); and KOs level analysis (d)

Fig. 6
figure 6

Linear discriminant analysis (LDA) effect size analysis (LEfSe) was used to comparison of KEGG functions differences in gut microbiota of FAW fed on different host plants. KEGG first level analysis (LDA > 2, P < 0.05) (a); KEGG second level analysis (LDA > 2, P < 0.05) (b); KEGG third level pathway analysis (LDA > 3, P < 0.05) (c) and different KOs enriched in the AD, CA, SL and SM group identified by LEfSe (LDA > 3, P < 0.05) (d). The abbreviations in the figure represent different diets, AD, artificial diet-fed; CA, pepper-fed; SL, tomato-fed; SM, sorghum-fed; ZM, maize-fed. KEGG, Kyoto Encyclopedia of Genes and Genomes; KOs, KEGG Orthologous groups. LDA, linear discriminant analysis; LEfSe, linear discriminant analysis effect size

Our analysis identified 57 pathway categories at the second level. The most prevalent was carbohydrate metabolism, constituting 13.56% of the total relative abundance, with the tomato-fed group exhibiting the highest relative abundance. Additionally, we found that amino acid metabolism made up the second largest proportion, with the sorghum-fed group having the highest relative abundance. The metabolism of cofactors and vitamins constituted the third largest proportion in all pathways (Fig. 5b, 6b, Supplementary Table S2). The top 10 functional annotations at the third-level pathways primarily involved metabolism pathways (terpenoid and steroid biosynthesis; D-alanine metabolism; starch and sucrose metabolism; valine, leucine, and isoleucine biosynthesis), genetic information processing (ribosome), environmental information processing (phosphotransferase system [PTS]), human disease (vancomycin resistance), and cellular processes (bacterial chemotaxis; peptidoglycan biosynthesis; and flagellar assembly) (Fig. 5c, Supplementary Table S3). Using the KEGG Orthology database and performing BLAST analysis, we obtained the top 10 most abundant KOs (Fig. 5d).

Further examination of annotated genes using the KEGG databases revealed significant differences between the groups across various pathways and functions. These differences were particularly evident in the third-level components of KEGG among the five groups. Specifically, the biosynthesis of terpenoid backbone and zeatin was associated with the metabolism of terpenoids and polyketides in the maize-fed group. Valine, leucine, and isoleucine biosynthesis, histidine metabolism associated with amino acid metabolism, and C5-branched dibasic acid metabolism pertain to carbohydrate metabolism. Furthermore, cell cycle-caulobacter, 2-oxocarboxylic acid metabolism, and limonene and pinene degradation were intricately associated with the sorghum-fed group. The PTS and fructose and mannose metabolism were notably linked to the tomato-fed group. Moreover, flagellar assembly, starch and sucrose metabolism, vancomycin resistance, other glycan degradation, and peptidoglycan biosynthesis were closely aligned with the pepper-fed group, whereas D-alanine metabolism, lipoic acid metabolism, streptomycin biosynthesis, thiamine metabolism, and pentose phosphate pathway exhibited strong correlations with the artificial diet-fed group (Fig. 6c, Supplementary Table S4).

Figure 6d illustrates several KOs with significant differences in relative abundance among the sorghum-fed, pepper-fed, and artificial diet-fed groups. Utilizing LEfSe (LDA > 3), we identified 13 KOs enriched in the pepper-fed group, 3 KOs enriched in the sorghum-fed group, and only 1 KO enriched in the artificial diet-fed group. Among the most abundant KOs in the sorghum-fed group were K02033, K02034, and K02032, which are involved in regulating peptide/nickel transport system permease protein and peptide/nickel transport system ATP-binding protein. These KOs are related to quorum sensing. Additionally, K01223, which regulates 6-phospho-β-glucosidase associated with glycolysis/gluconeogenesis and starch and sucrose metabolism, was enriched in the pepper-fed group. Notably, only one KO (K03704: cold shock protein) was enriched in the artificial diet-fed group but was not associated with any specific pathway. To further explore the function of different KOs, we initially annotated the KO gene, clarified the pathway, and then selected representative pathway diagrams for presentation, as depicted in Fig. 7.

Fig. 7
figure 7

The biomarker genes across the AD, CA, SL and SM group in the Glycolysis / Gluconeogenesis pathway (map00010). The abbreviations in the figure represent different diets, AD, artificial diet-fed; CA, pepper-fed; SL, tomato-fed; SM, sorghum-fed; ZM, maize-fed. The gray rectangular box indicates the detection of the gene in the sample, while the colored rectangular box represents characteristic genes corresponding to specific color-coded groups

Following the selection of feature functions across different groups, we identified functional species using the results from HUMAnN2 (Fig. 8). E. mundtii was found to contribute to the artificial diet-fed group (Fig. 8a). E. casseliflavus was linked to the metabolic dynamics of the pepper-fed group (Fig. 8b) and was similarly associated with metabolism and environmental information processing in the tomato-fed group (Fig. 8c). Agrobacterium tumefaciens was found to be associated with quorum sensing pathways, which were correlated with cellular processes in the sorghum-fed group (Fig. 8d). However, it is important to note that while this association is observed, further experimental evidence is needed to establish a causal link between A. tumefaciens and quorum sensing in this context. Similarly, Enterococcus casseliflavus was identified to have an association with the regulation of the energy-coupling factor transport system permease protein, which is linked to ABC transporters that were more abundant in the maize-fed group (Fig. 8e). Again, this finding points to an association rather than a direct causal relationship, emphasizing the need for additional experimental verification.

Fig. 8
figure 8

Bacterial species contributed to key functions. The ten bacterial species contributed to cold shock protein (K03704) enrichment in the artificial diet-fed (AD) group (a); The ten bacterial species contributed to 6-phospho-β-glucosidase (K01223) enrichment in the pepper-fed (CA) group (b); The ten bacterial species contributed to fructose PTS system EIIBC or EIIC component (K02770) enrichment in the tomato-fed (SL) group (c); The ten bacterial species contributed to peptide/nickel transport system permease protein (K02033) enrichment in the sorghum-fed (SM) group (d). The ten bacterial species contributed to energy-coupling factor transport system permease protein (K16785) enrichment in the maize-fed (ZM) group (e)

To identify the top 20 ARGs based on relative abundance among 15 samples, we compared our data with the Comprehensive Antibiotic Resistance Database. We calculated the content and percentage of antibiotic resistance ontology (ARO) in each sample and selected the top 20 ARO results with the highest abundance for display (Supplementary Figure S1). ARO:3002922 (vanRC) and ARO:3003948 (efrA) were found to be predominant across the 15 samples, indicating a high level of resistance to antibiotics. Based on the LEfSe analysis, the dominant ARGs of FAW feeding on maize were vanRC, vanSC, and vanC, which confer resistance to glycopeptides. For FAW feeding on pepper, the dominant ARG was TetM, which encodes a ribosomal protection protein that confers resistance to tetracycline. The dominant gut ARG of FAW larvae feeding on sorghum was CAT, a chloramphenicol resistance gene. When the larvae were fed an artificial diet, the dominant ARG was dfrE, a chromosome-encoded dihydrofolate reductase that confers resistance to diaminopyrimidine antibiotics (Fig. 6b). Supplementary Table S5 lists each resistance gene type and its corresponding antibiotics. The study found that the number of glycopeptide antibiotics was the highest among all samples.

Discussion

In this study, we examined the structure, function, and ARGs of the gut microbiota in FAW. The larvae were fed on four host plants: two Poaceae crops (maize and sorghum), two Solanaceae plants (pepper and tomato), and an artificial diet under laboratory conditions. Our findings suggest that the diversity and abundance of the FAW larval gut microbial community, its function, and ARGs are influenced by the host plants. This study is, to our knowledge, the first to use metagenomic analysis, as opposed to the commonly used 16S rRNA amplicon sequencing, to study the gut microbiota in FAW influenced by different host plants.

Previous studies have indicated that the structure of insect gut bacteria can be affected by the host diet and species [43, 44, 55]. Our study supports these findings, showing that the diversity of gut microbiota changes to varying degrees when transitioning from an artificial diet to different host plants. The gut microbiota richness of FAW increased after transitioning from an artificial diet to sorghum. However, when transferred to maize, tomatoes, and especially pepper, the gut microbial diversity was notably lower. This suggests that Poaceae and Solanaceae crops provide different nutritional substances for the growth and development of FAW larvae.

Our previous experimental findings suggest that among the Solanaceae crops consumed, the consumption of pepper by the FAW leads to the shortest adult stage duration, the longest preadult stage, and the smallest pupal weight. In contrast, when compared to the consumption of the Poaceae crop sorghum, there are significant differences in both the adult stage, the preadult stage, and pupal weight [56]. These results align with a previous finding that the microbial diversity of the FAW larvae fed on pepper was the lowest [24].

Significant variations in gut microbiota diversity at the phylum level have been observed. Firmicutes is the dominant phylum in the gut microbiota of FAW larvae feeding on maize, tomato, pepper, and an artificial diet. This aligns with previous FAW studies [21, 22, 24, 34, 57,58,59,60], indicating that Firmicutes could be crucial for the development and fitness of FAW. However, in FAW larvae feeding on sorghum, the dominant phylum is Proteobacteria, followed by Firmicutes. This is in contrast to a previous study suggesting that differences in sorghum varieties could lead to this discrepancy. Furthermore, discrepancies may arise from various growth conditions of sorghum, including temperature, humidity, and light cycles, as well as differing physiological states. Similarly, variations in the fall armyworm's provenance, genetic lineage, and rearing environments are also likely to contribute to these inconsistencies [21]. Studies have shown that the bacterial diversity and abundance in the guts of FAW larvae are influenced by six different corn varieties [21]. Different host plants have varying nutritional substances, digestion facilitation and secondary metabolites, which may affect the growth and development of the fall armyworm [60]. The significant differences in the gut microbiota of fall armyworm larvae feeding on different host plants may be due to the different nutritional components and inhibitory secondary metabolites of these four host plants. Exploring the complex dynamics between crop diversity and insect gut microbiota is crucial for the development of advanced pest management strategies. Targeting the gut microbiota of sorghum-fed pests with specialized biopesticides presents a refined approach to enhance pest control efficacy [61]. Develop effective pest control strategies by understanding the preferences and selection behaviors of pests towards their host plants. Strategically designed crop rotation and intercropping systems, when precisely implemented, can effectively mitigate pest issues by disrupting their life cycles and foraging behaviors, thereby diminishing pest prevalence and damage [62]. Firmicutes, a dominant phylum in the gut microbiota, are known for their ability to encode enzymes that are integral to energy metabolism and potentially the biosynthesizing of vitamin B. They play pivotal roles in nutritional supplementation, energy absorption, and host immunity [63,64,65]. On the other hand, Proteobacteria, another significant phylum, are associated with a broad spectrum of metabolic processes. These include the decomposition and fermentation of complex sugars and the production of vitamins. They also play an essential role in the degradation of various aromatic compounds, thereby enhancing the absorption of nutrients by the host [64, 66]. Both these phyla, Firmicutes and Proteobacteria, are crucial in regulating carbohydrate and amino acid metabolism. They also facilitate host membrane transport pathways [30, 67, 68]. The stable colonization of gut bacteria, particularly Proteobacteria, is known to play a functionally crucial role in the adaptation of insects to specific host plants [43]. However, the role of Proteobacteria in the metabolic pathway of FAW larvae and their adaptation to hosts requires further research.

Enterococcaceae, a unique bacterial family, exhibits a significant presence in the gut microbiota of FAW larvae that feed on tomato, pepper, and maize. In contrast, its presence is relatively low in the gut microbiota of FAW larvae that feed on an artificial diet and sorghum, particularly in the sorghum-fed group. Enterococcaceae is known to play a crucial role in degrading alkaloids and latex, thereby stabilizing resistant plants [44]. This bacterial family is also found in other insects such as Hyphantria cunea (Lepidoptera: Erebidae) [69], Galleria mellonella (Lepidoptera: Pyralidae) [70], and Paederus fuscipes (Coleoptera: Staphylinidae) [71].

On the other hand, Lactobacillaceae was detected exclusively in the gut microbiota samples of the artificial diet-fed group but not in the FAW larvae fed on four different host plants: maize, sorghum, pepper, and tomato. Lactobacillaceae has also been detected in other insects, such as Grapholita molesta (Lepidoptera: Tortricidae) [19], and Drosophila melanogaster [72, 73]. A previous study found that Lactobacillus plantarum promotes the growth and development of germ-free fruit flies, an effect that is related to the Tor and insulin signaling pathways [74]. Schretter et al. reported that monocolonization of germ-free Drosophila with Lactobacillus brevis restores their overactivated behavior to normal levels, similar to those of conventional flies [75]. The xylose isomerase expressed by L. brevis plays a key role in this process. Further research is required to determine if L. plantarum can promote growth and oviposition in FAW and if L. brevis can influence the locomotor behavior of FAW.

Enterobacteriaceae have been consistently observed in the gut microbiota of FAW larvae as well as in other insect species. Huang et al. demonstrated its role in accelerating reproduction in Harmonia axyridis, particularly by reducing the preoviposition duration and incubation period [64]. These findings highlight the importance of microbial communities for insect reproductive biology. Research indicates that Enterobacteriaceae, a symbiotic microorganism prevalent in both human and insect guts, plays pivotal roles in vitamin biosynthesis, pheromone production, and degradation of plant compounds [76]. Additionally, it contributes to enzyme metabolism through activities such as superoxide dismutase or catalase enzyme activity [40]. Moreover, Enterobacteriaceae have been linked to sugar metabolism in hosts and may aid in various physiological processes, including digestion, protection, courtship, and reproduction [77].

Enterococcus emerges as the most abundant genus observed in the gut microbiota of FAW larvae, serving as a stable insect symbiotic bacterium throughout metamorphosis. Its protective effects against pathogens [78], capability to neutralize toxic molecules in plants [79], and facilitation of host adaptation [80] contribute to increasing tolerance to harmful foods [81]. Additionally, Enterococcus plays roles in vitamin and pheromone synthesis, degradation of plant compounds, nitrogen fixation and insecticide [61, 82, 83].

Significant differences in gut microbial community structures of fifth-instar FAW larvae were observed based on the host plant they fed on, as revealed via PCoA. Notably, there was some overlap in the PCoA plots between pepper and maize, suggesting that these plants may offer unique nutrients crucial for FAW larvae growth and development. Ugwu et al. reported that the observed variations in the gut of FAW larvae could be attributed to the diverse chemical components of the different diet groups [23]. These findings underscore the importance of understanding the relationship between host plant nutrition and insect microbiomes for effective FAW infestation management.

In the investigation of the gut microbiota of FAW larvae feeding on various host plants, researchers commonly utilize 16S amplicon sequencing technology for analysis. Additionally, functional annotation is predicted and analyzed using Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt) based on OTU representative sequences. Previous studies have indicated variations in the gut bacteria of FAW larvae fed on different host plants, with the most dominant putative function prediction categories in metabolic pathways [21, 24, 57, 58]. In this study, metagenomic sequencing was employed to accurately decipher the microbial community functions and obtain functional abundance and differences using the KEGG database. Results revealed that the gut microbiota of FAW larvae was predominantly associated with metabolism, followed by genetic information processing and cellular processes. At the KEGG second-level pathway, the FAW gut microbiota was primarily involved in carbohydrate metabolism, amino acid metabolism, and metabolism of cofactors and vitamins. This suggests a significant role of the microbiota in food digestion, degradation, and nutrient absorption. Additionally, the microbiota plays a crucial role in vitamin metabolism, providing essential vitamins for host growth and development.

At the KEGG third-level pathway, various cellular processes were observed, including the biosynthesis of terpenoids and steroids; ribosomal functions; PTS; metabolism of D-alanine; starch and sucrose metabolism; biosynthesis of valine, leucine, and isoleucine; vancomycin resistance; bacterial chemotaxis; biosynthesis of peptidoglycan; flagellar assembly; and the biosynthesis of amino acids. Tang et al. explored the variations in gut microbiota composition and functional annotation between S. litura and S. frugiperda, predicting that S. frugiperda had a significantly higher relative abundance of gut microbiota involved in detoxification compared to S. litura [84]. Additionally, pathways related to detoxification were found to be more abundant in S. frugiperda.

Furthermore, we identified several ARGs, indicating that the gut microbiota of FAW may serve as a reservoir for these genes. Our analysis revealed the top 20 ARGs, with vanRC being the most abundant in the gut microbiota of maize-fed FAW. These ARGs confer resistance to antibiotics such as glycopeptide, macrolide, fluoroquinolone, rifamycin, and rifampicin antibiotics, underscoring the potential for antibiotic resistance development within the FAW gut microbiota. This has significant implications for pest management strategies. Interestingly, Tang et al. detected different types of FAW resistance genes, with the relative abundances of AbaF, DfrA42, and MexD being most abundant, contrasting with our study results [84]. These findings highlight the prevalence and diversity of ARGs in the gut microbiota of FAW larvae and underscore the importance of host plant selection in shaping antibiotic resistance profiles. Further research in this area is crucial for a comprehensive understanding of the factors influencing antibiotic resistance in agricultural ecosystems.

In this study, we only observed one generation and one instar (fifth instar) after the host switch. Long-term observation of successive generations of FAW larvae fed on host plants would provide a more comprehensive understanding of the dynamic changes in the gut microbiome during the process of host transfer. Additionally, employing multiomics technologies and conducting bacterial removal and re-inoculation experiments are necessary to elucidate the specific functions of distinct strains of gut microbiota in FAW adaptation to different host plants. This would ultimately enhance our understanding of the relationship between FAW fitness and its host plants.

Conclusion

In conclusion, our study utilized metagenomic sequencing to explore the gut microbiota of FAW larvae fed on various host plants and an artificial diet. We observed significant variations in gut microbiota composition among different host plants, with Firmicutes being the dominant phylum, Enterococcaceae as the dominant family, and Enterococcus as the predominant genus. Notably, Enterococcus casseliflavus emerged as the most abundant species. Functional analysis indicated a predominance of metabolic pathways, particularly those related to carbohydrate and amino acid metabolism. We identified vanRC as the most abundant ARG, conferring resistance to glycopeptide antibiotics. These findings enhance our understanding of gut microbiota dynamics and functional profiles in FAW larvae across different host plants, offering valuable insights into pest management strategies.

Availability of data and materials

All metagenome sequence reads can be accessed in the Sequence Read Archive (SRA) of NCBI under the BioProject with the accession number PRJNA1123873, available at https://www.ncbi.nlm.nih.gov/sra/PRJNA1123873.

Abbreviations

FAW:

Fall Armyworm

ARG:

Antibiotic Resistance Gene

PCR:

Polymerase Chain Reaction

qPCR:

Quantitative Polymerase Chain Reaction

PE:

Paired-End

NCBI:

National Center of Biotechnology Information

PCoA:

Principal Coordinate Analysis

BRACKen:

Bayesian Re-estimation of Abundance after Classification with Kraken

LDA:

Linear Discriminant Analysis

LEfSe:

Linear Discriminant Analysis (LDA) Effect Size

UniRef90:

UniProt Reference Clusters 90

HUMAnN2:

HMP Unified Metabolic Analysis Network

KEGG:

Kyoto Encyclopedia of Genes and Genomes

KOs:

KEGG Orthologs

CARDS:

Comprehensive Antibiotic Resistance Database

OTUs:

Operational Taxonomic Units

ARO:

Antibiotic Resistance Ontology

PICRUSt:

Phylogenetic Investigation of Communities by Reconstruction of Unobserved States

PTS:

Phosphotransferase System

References

  1. Sparks AN. A review of the biology of the fall armyworm. Fla Entomo. 1979;62(2):82–7. https://doi.org/10.2307/3494083.

    Article  Google Scholar 

  2. Todd EL, Poole RW. Keys and illustrations for the armyworm moths of the noctuid genus Spodoptera Guenée from the Western Hemisphere. Ann Entomol Soc Am. 1980;73(6):722–38. https://doi.org/10.1093/aesa/73.6.722.

    Article  Google Scholar 

  3. Martinelli S, Barata RM, Zucchi MI, Silva-Filho MC, Omoto C. Molecular variability of Spodoptera frugiperda (Lepidoptera: Noctuidae) populations associated to maize and cotton crops in Brazil. J Econ Entomol. 2006;99(2):519–26. https://doi.org/10.1603/0022-0493-99.2.519.

    Article  CAS  PubMed  Google Scholar 

  4. Goergen G, Kumar PL, Sankung SB, Togola A, Tamò M. First report of outbreaks of the fall armyworm Spodoptera frugiperda (J E Smith) (Lepidoptera, Noctuidae), a new alien invasive pest in West and Central Africa. PLoS ONE. 2016;11(10):e0165632. https://doi.org/10.1371/journal.pone.0165632.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Cock MJW, Beseh PK, Buddie AG, Cafá G, Crozier J. Molecular methods to detect Spodoptera frugiperda in Ghana, and implications for monitoring the spread of invasive species in developing countries. Sci Rep-UK. 2017;7(1):4103. https://doi.org/10.1038/s41598-017-04238-y.

    Article  CAS  Google Scholar 

  6. Midega CAO, Pittchar JO, Pickett JA, Hailu GW, Khan ZR. A climate-adapted push-pull system effectively controls fall armyworm, Spodoptera frugiperda (J E Smith), in maize in East Africa. Crop Prot. 2018;105:10–5. https://doi.org/10.1016/j.cropro.2017.11.003.

    Article  Google Scholar 

  7. Gurung K, Wertheim B, Falcao Salles J. The microbiome of pest insects: It is not just bacteria. Entomol Exp Appl. 2019;167(3):156–70. https://doi.org/10.1111/eea.12768.

    Article  Google Scholar 

  8. Paniagua Voirol LR, Frago E, Kaltenpoth M, Hilker M, Fatouros NE. Bacterial symbionts in Lepidoptera: their diversity, transmission, and impact on the host. Front Microbiol. 2018;9:556. https://doi.org/10.3389/fmicb.2018.00556.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Rangberg A, Diep DB, Rudi K, Amdam GV. Paratransgenesis: an approach to improve colony health and molecular insight in honey bees (Apis mellifera)? Integr Comp Biol. 2012;52(1):89–99. https://doi.org/10.1093/icb/ics089.

    Article  PubMed  Google Scholar 

  10. Engel P, Moran NA. The gut microbiota of insects-diversity in structure and function. FEMS Microbiol Rev. 2013;37(5):699–735. https://doi.org/10.1111/1574-6976.12025.

    Article  CAS  PubMed  Google Scholar 

  11. Fraune S, Bosch TC. Why bacteria matter in animal development and evolution. BioEssays. 2010;32(7):571–80. https://doi.org/10.1002/bies.200900192.

    Article  CAS  PubMed  Google Scholar 

  12. Janson EM, Stireman JO 3rd, Singer MS, Abbot P. Phytophagous insect-microbe mutualisms and adaptive evolutionary diversification. Evolution. 2008;62(5):997–1012. https://doi.org/10.1111/j.1558-5646.2008.00348.x.

    Article  PubMed  Google Scholar 

  13. Hansen AK, Moran NA. The impact of microbial symbionts on host plant utilization by herbivorous insects. Mol Ecol. 2014;23(6):1473–96. https://doi.org/10.1111/mec.12421.

    Article  PubMed  Google Scholar 

  14. He M, Chen H, Yang X, Gao Y, Lu Y, Cheng D. Gut bacteria induce oviposition preference through ovipositor recognition in fruit fly. Commun Biol. 2022;5(1):973. https://doi.org/10.1038/s42003-022-03947-z.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Zhu P, Mamdouh B, Magdy M, Chen J, Ma W, Xu L. Diet induced variation in gut microbiota is linked to the growth performance of an agricultural pest Chilo suppressalis. Agronomy-basel. 2023;13(2):304. https://doi.org/10.3390/agronomy13020304.

    Article  CAS  Google Scholar 

  16. Yang Y, Liu X, Xu H, Liu Y, Lu Z. Effects of host plant and insect generation on shaping of the gut microbiota in the rice leaffolder, Cnaphalocrocis medinalis. Front Microbiol. 2022;13:824224. https://doi.org/10.3389/fmicb.2022.824224.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Meng L, Xia C, Jin Z, Zhang H. Investigation of gut bacterial communities of Asian citrus psyllid (Diaphorina citri) reared on different host plants. Insects. 2022;13(8):694. https://doi.org/10.3390/insects13080694.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Wang Q, Liu Y, Yin X. Comparison of gut bacterial communities of Locusta migratoria manilensis (Meyen) reared on different food plants. Biology-Basel. 2022;11(9):1347. https://doi.org/10.3390/biology11091347.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Yuan XQ, Zhang X, Liu X, Dong YL, Yan ZZ, Lv DB, Wang P, Li YP. Comparison of gut bacterial communities of Grapholita molesta (Lepidoptera: Tortricidae) reared on different host plants. Int J Mol Sci. 2021;22(13):6843. https://doi.org/10.3390/ijms22136843.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Ma L, Wang DT, Ren QL, Zhang L, Cheng YX, Jiang XF. Feeding experience on hairy vetch induces a loss of host-specifc adaptation to maize and changes in the gut microbiota communities of the fall armyworm. J Pest Sci on line. 2023. https://doi.org/10.1007/s10340-023-01682-w.

    Article  Google Scholar 

  21. Wang YP, Liu X, Yi CY, Chen XY, Liu CH, Zhang CC, et al. The adaptive evolution in the fall armyworm Spodoptera frugiperda (Lepidoptera: Noctuidae) revealed by the diversity of larval gut bacteria. Genes-Basel. 2023;14(2):321. https://doi.org/10.3390/genes14020321.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. de Oliveira NC, Cônsoli FL. Dysbiosis of the larval gut microbiota of Spodoptera frugiperda strains feeding on different host plants. Symbiosis. 2023;89(2):197–211. https://doi.org/10.1007/s13199-023-00907-x.

    Article  Google Scholar 

  23. Ugwu JA, Wenzi R, Asiegbu FO. Monocot diet sources drive diversity of gut bacterial communities in Spodoptera frugiperda (Lepidoptera: Noctuidae) larvae. J Appl Entomol. 2022;146(8):942–56. https://doi.org/10.1111/jen.13022.

    Article  CAS  Google Scholar 

  24. LÇš DB, Liu XY, Dong YL, Yan ZZ, Zhang X, Wang P, et al. Comparison of gut bacterial communities of fall armyworm (Spodoptera frugiperda) reared on different host plants. Int J Mol Sci. 2021;22(20):11266. https://doi.org/10.3390/ijms222011266.

    Article  CAS  Google Scholar 

  25. Acevedo FE, Peiffer M, Tan CW, Stanley BA, Stanley A, Wang J, Jones AG, Hoover K, Rosa C, Luthe D, et al. Fall armyworm-associated gut bacteria modulate plant defense responses. Mol Plant Microbe Interact. 2017;30(2):127–37. https://doi.org/10.1094/MPMI-11-16-0240-R.

    Article  CAS  PubMed  Google Scholar 

  26. Almeida LG, Moraes LA, Trigo JR, Omoto C, Cônsoli FL. The gut microbiota of insecticide-resistant insects houses insecticide-degrading bacteria: A potential source for biotechnological exploitation. Plos one. 2017;12(3):e0174754. 82. https://doi.org/10.1371/journal.pone.0174754.

    Article  CAS  Google Scholar 

  27. Mason CJ, Ray S, Shikano I, Peiffer M, Jones AG, Luthe DS, Hoover K, Felton GW. Plant defenses interact with insect enteric bacteria by initiating a leaky gut syndrome. P Natl Acad Sci USA. 2019;116(32):15991–6. https://doi.org/10.1073/pnas.1908748116.

    Article  CAS  Google Scholar 

  28. Mason CJ, St Clair A, Peiffer M, Gomez E, Jones AG, Felton GW, et al. Diet influences proliferation and stability of gut bacterial populations in herbivorous lepidopteran larvae. PLoS ONE. 2020;15(3):e0229848. https://doi.org/10.1371/journal.pone.0229848.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Gomes AFF, Omoto C, Cônsoli FL. Gut bacteria of field-collected larvae of Spodoptera frugiperda undergo selection and are more diverse and active in metabolizing multiple insecticides than laboratory-selected resistant strains. J Pest Sci. 2020;93(2):833–51. https://doi.org/10.1007/s10340-020-01202-0.

    Article  Google Scholar 

  30. Wang XL, Sun SJ, Yang XL, Cheng J, Wei HS, Li Z, et al. Variability of gut microbiota across the life cycle of Grapholita molesta (Lepidoptera: Tortricidae). Front Microbiol. 2020;11:1366. https://doi.org/10.3389/fmicb.2020.01366.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Jones AG, Mason CJ, Felton GW, Hoover K. Host plant and population source drive diversity of microbial gut communities in two polyphagous insects. Sci Rep-UK. 2019;9(11):2792. https://doi.org/10.1038/s41598-019-39163-9.

    Article  CAS  Google Scholar 

  32. Rozadilla G, Cabrera NA, Virla EG, Greco NM, McCarthy CB. Gut microbiota of Spodoptera frugiperda (J. E. Smith) larvae as revealed by metatranscriptomic analysis. J Appl Entomol. 2020;144(19):351–63. https://doi.org/10.1111/jen.12742.

    Article  CAS  Google Scholar 

  33. Gichuhi J, Sevgan S, Khamis F, Van den Berg J, du Plessis H, Ekesi S, Herren JK. Diversity of fall armyworm, Spodoptera frugiperda and their gut bacterial community in Kenya. Peer J. 2020;8:e8701. https://doi.org/10.7717/peerj.8701.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Ugwu JA, Liu MX, Sun H, Asiegbu FO. Microbiome of the larvae of Spodoptera frugiperda (JE Smith) (Lepidoptera: Noctuidae) from maize plants. J Appl Entomol. 2020;114(9):764–76. https://doi.org/10.1111/jen.12821.

    Article  CAS  Google Scholar 

  35. Mason CJ, Hoover K, Felton GW. Effects of maize (Zea mays) genotypes and microbial sources in shaping fall armyworm (Spodoptera frugiperda) gut bacterial communities. Sci Rep-UK. 2021;11(1):4429. https://doi.org/10.1038/s41598-021-83497-2.

    Article  CAS  Google Scholar 

  36. Rausch P, Rühlemann M, Hermes BM, Doms S, Dagan T, Dierking K, et al. Comparative analysis of amplicon and metagenomic sequencing methods reveals key features in the evolution of animal metaorganisms. Microbiome. 2019;7(1):133. https://doi.org/10.1186/s40168-019-0743-1.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Yk Z, Ke Mj YuZT, Ct L, Liu M, Yang Yh. Combined effects of azoxystrobin and oxytetracycline on rhizosphere microbiota of Arabidopsis thaliana. Environ Int. 2024;186:108655. https://doi.org/10.1016/j.envint.108655.

    Article  Google Scholar 

  38. Shi LM, Zhang JY, Lu TD, Zhang KC. Metagenomics revealed the mobility and hosts of antibiotic resistance genes in typical pesticide wastewater treatment plants. Sci Total Environ. 2022;817:153033. https://doi.org/10.1016/j.scitotenv.153033.

    Article  CAS  PubMed  Google Scholar 

  39. Taunk J, & Goutam U. Role of metagenomics in plant disease management. In R. Prasad, V. Kumar, J. Singh, & C. P. Upadhyaya (Eds.), Recent developments in microbial technologies. Environmental and Microbial Biotechnology. Springer Nature Singapore Pte Ltd. 203–220;2021. https://doi.org/10.1007/978-981-15-4439-2

  40. Xia X, Gurr GM, Vasseur L, Zheng D, Zhong H, Qin B. Metagenomic sequencing of diamondback moth gut microbiome unveils key holobiont adaptations for herbivory. Front Microbiol. 2017;8:663. https://doi.org/10.3389/fmicb.2017.00663.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Li C, Han G, Sun J, Huang L, Lu Y, Xia Y, Liu Q, Xu J. The gut microbiota composition of cnaphalocrocis medinalis and their predicted contribution to larval nutrition. Front Microbiol. 2022;13:909863. https://doi.org/10.3389/fmicb.2022.909863.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Di XY, Yan B, Wu CX, Yu XF, Liu JF, Yang MF. Does larval rearing diet lead to premating isolation in Spodoptera litura (Fabricius) (Lepidoptera: Noctuidae)? Insects. 2021;12(3):203. https://doi.org/10.3390/insects12030203.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Yang FY, Saqib HSA, Chen JH, Ruan QQ, Vasseur L, He WY, You MS. Differential profiles of gut microbiota and metabolites associated with host shift of Plutella xylostella. Int J Mol Sci. 2020;21(17):6283. https://doi.org/10.3390/ijms21176283.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Liu Y, Shen Z, Yu J, Li Z, Liu X, Xu H. Comparison of gut bacterial communities and their associations with host diets in four fruit borers. Pest Manag Sci. 2020;76(4):1353–62. https://doi.org/10.1002/ps.5646.

    Article  CAS  PubMed  Google Scholar 

  45. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114-2120. 82. https://doi.org/10.1093/bioinformatics/btu170.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9(D1):357–9. https://doi.org/10.1038/nmeth.1923.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Wood DE, Lu J, Langmead B. Improved metagenomic analysis with Kraken2. Genome Biol. 2019;20(1):257. https://doi.org/10.1186/s13059-019-1891-0.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Lu J, Breitwieser FP, Thielen P, Salzberg SL. Bracken: Estimating species abundance in metagenomics data. Peerj Comput Sci. 2017;3:e104. https://doi.org/10.7717/peerj-cs.104.

    Article  Google Scholar 

  49. Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, McGlinn D et al. Package ‘Vegan’. community ecology package, version 2.5–7. 2019;2. https://github.com/vegandevs/vegan

  50. Liu YX, Qin Y, Chen T, Lu MP, Qian XB, Guo XX, et al. A practical guide to amplicon and metagenomic analysis of microbiome data. Protein Cell. 2021;12(5):315–30. https://doi.org/10.1007/s13238-020-00724-8.

    Article  PubMed  Google Scholar 

  51. Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, et al. Metagenomic biomarker discovery and explanation. Genome Biol. 2011;12(6):R60. https://doi.org/10.1186/gb-2011-12-6-r60.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2015;12(1):59–60. https://doi.org/10.1038/nmeth.3176.

    Article  CAS  PubMed  Google Scholar 

  53. Franzosa EA, McIver LJ, Rahnavard G, Thompson LR, Schirmer M, Weingart G, Lipson KS, Knight R, Caporaso JG, Segata N, Huttenhower C. Species-level functional profiling of metagenomes and metatranscriptomes. Nat Methods. 2018;15(11):962–8. https://doi.org/10.1038/s41592-018-0176-y.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Kanehisa M, Furumichi M, Tanabe M, Sato Y, Morishima K. KEGG: new perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res. 2017;45(D1):D353–61. https://doi.org/10.1093/nar/gkw1092.

    Article  CAS  PubMed  Google Scholar 

  55. Chen BS, Du KQ, Sun C, Vimalanathan A, Liang XL, Li Y, Wang BH, Lu XM, Li LJ, Shao YQ. Gut bacterial and fungal communities of the domesticated silkworm (Bombyx mori) and wild mulberry-feeding relatives. Isme J. 2018;12(9):2252–62. https://doi.org/10.1038/s41396-018-0174-1.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Wu LH, 2022. Effects of host plants on the development, reproduction, mating behavior and gut microbiota of fall armyworm (Spodoptera frugiperda). Master Thesis, Guizh University, Guiyang, Guizhou, China. https://doi.org/10.27047/d.cnki.ggudu.2022.001768

  57. Han WK, Tang FX, Yan YY, Wang Y, Liu ZW. Plasticity of the gene transcriptional level and microbiota in the gut contributes to the adaptability of the fall armyworm to rice plants. J Agric Food Chem. 2023;71(47):18546–56. https://doi.org/10.1021/acs.jafc.3c05506.

    Article  CAS  PubMed  Google Scholar 

  58. Li DD, Li JY, Hu ZQ, Liu TX, Zhang SZ. Fall armyworm gut bacterial diversity associated with different developmental stages, environmental habitats, and diets. Insects. 2022;13(9):762. https://doi.org/10.3390/insects13090762.

    Article  PubMed  PubMed Central  Google Scholar 

  59. Chen YP, Li YH, Sun ZX, Du EW, Lu ZH, Li H, Gui FR. Effects of host plants on bacterial community structure in larvae midgut of Spodoptera frugiperda. Insects. 2022;13(4):373. https://doi.org/10.3390/insects13040373.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Wu LH, Zhou C, Long GY, Yang XB, Wei ZY, Liao YJ, Yang H, Hu CX. Fitness of fall armyworm, Spodoptera frugiperda to three solanaceous vegetables. J Integr Agr. 2021;20(3):755–63. https://doi.org/10.1016/S2095-3119(20)63476-1.

    Article  Google Scholar 

  61. Gomes AFF, de Almeida LG, Cônsoli FL. Comparative genomics of pesticide-degrading enterococcus symbionts of Spodoptera frugiperda (Lepidoptera: Noctuidae) leads to the identification of two new species and the reappraisal of insect-associated enterococcus species. Microb Ecol. 2023;86(4):2583–605. https://doi.org/10.1007/s00248-023-02264-0.

    Article  CAS  PubMed  Google Scholar 

  62. Jeon J, Rahman MM, Han C, Shin J, Sa KJ, Kim J. Spodoptera frugiperda (Lepidoptera: Noctuidae) life table comparisons and gut microbiome analysis reared on corn varieties. Insects. 2023;14(4):358. https://doi.org/10.3390/insects14040358.

    Article  PubMed  PubMed Central  Google Scholar 

  63. Flint HJ, Scott KP, Duncan SH, Louis P, Forano E. Microbial degradation of complex carbohydrates in the gut. Gut Microbes. 2012;3(4):289–306. https://doi.org/10.4161/gmic.19897.

    Article  PubMed  PubMed Central  Google Scholar 

  64. Colston TJ, Jackson CR. Microbiome evolution along divergent branches of the vertebrate tree of life: what is known and unknown. Mol Ecol. 2016;25(16):3776–800. https://doi.org/10.1111/mec.13730.

    Article  PubMed  Google Scholar 

  65. Huang ZD, Zhu L, Lv J, Pu ZX, Zhang LP, Chen GQ, et al. Dietary effects on biological parameters and gut microbiota of Harmonia axyridis. Front Microbiol. 2022;12: 818787. https://doi.org/10.3389/fmicb.2021.818787.

    Article  PubMed  PubMed Central  Google Scholar 

  66. Reid NM, Addison SL, Macdonald LJ, Lloyd-Jones G. Biodiversity of active and inactive bacteria in the gut flora of wood-feeding Huhu beetle larvae (Prionoplus reticularis). Appl Environ Microb. 2010;77(19):7000–6. https://doi.org/10.1128/AEM.05609-11.

    Article  CAS  Google Scholar 

  67. Genta FA, Dillon RJ, Terra WR, Ferreira C. Potential role for gut microbiota in cell wall digestion and glucoside detoxification in Tenebrio molitor larvae. J Insect Physiol. 2006;52(6):593–601. https://doi.org/10.1016/j.jinsphys.2006.02.007.

    Article  CAS  PubMed  Google Scholar 

  68. Warnecke F, Luginbühl P, Ivanova N, Ghassemian M, Richardson TH, Stege JT, et al. Metagenomic and functional analysis of hindgut microbiota of a wood-feeding higher termite. Nature. 2007;450(7169):560–5. https://doi.org/10.1038/nature06269.

    Article  CAS  PubMed  Google Scholar 

  69. Chen HJ, Hao DJ, Chen CY, Sun YH, Yu XH. Effects of midgut bacteria in Hyphantria cunea (Lepidoptera: Erebidae) on nuclear polyhedrosis virus and Bacillus thuringiensis (Bacillales: Bacillaceae). J Insect Sci. 2023;23(2):1. https://doi.org/10.1093/jisesa/iead009.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Noël G, Serteyn L, Sare AR, Massart S, Delvigne F, Francis F. Co-diet supplementation of low density polyethylene and honeybee wax did not influence the core gut bacteria and associated enzymes of Galleria mellonella larvae (Lepidoptera: Pyralidae). Int Microbiol. 2023;26(2):397–409. https://doi.org/10.1007/s10123-022-00303-3.

    Article  CAS  PubMed  Google Scholar 

  71. Chamankar B, Maleki-Ravasan N, Karami M, Forouzan E, Karimian F, Naeimi S, Choobdar N. The structure and diversity of microbial communities in Paederus fuscipes (Coleoptera: Staphylinidae): from ecological paradigm to pathobiome. Microbiome. 2023;11(1):11. https://doi.org/10.1186/s40168-022-01456-z.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  72. Shin SC, Kim SH, You H, Kim B, Kim AC, Lee KA, et al. Drosophila microbiome modulates host developmental and metabolic homeostasis via insulin signaling. Science. 2011;334(6056):670–4. https://doi.org/10.1126/science.1212782.

    Article  CAS  PubMed  Google Scholar 

  73. Tefit MA, Leulier F. Lactobacillus plantarum favors the early emergence of fit and fertile adult Drosophila upon chronic undernutrition. J Exp Biol. 2017;220(5):900–7. https://doi.org/10.1242/jeb.151522.

    Article  PubMed  PubMed Central  Google Scholar 

  74. Storelli G, Defaye A, Erkosar B, Hols P, Royet J, Leulier F. Lactobacillus plantarum promotes Drosophila systemic growth by modulating hormonal signals through TOR-dependent nutrient sensing. Cell Metab. 2011;14(3):403–14. https://doi.org/10.1016/j.cmet.2011.07.012.

    Article  CAS  PubMed  Google Scholar 

  75. Schretter CE, Vielmetter J, Bartos I, Marka Z, Marka S, Argade S, Mazmanian SK. A gut microbial factor modulates locomotor behaviour in Drosophila. Nature. 2018;563(7731):402–6. https://doi.org/10.1038/s41586-018-0634-9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  76. Schloss PD, Delalibera IJ, Handelsman Jo, Raffa KF. Bacteria associated with the guts of two wood-boring beetles: Anoplophora glabripennis and Saperda vestita (Cerambycidae). Environ Entomol. 2006;35(3):625–9. https://doi.org/10.1603/0046-225X-35.3.625.

    Article  Google Scholar 

  77. Zhao XF, Zhang XY, Chen ZS, Wang Z, Lu YY, Cheng DF. The divergence in bacterial components associated with Bactrocera dorsalis across developmental stages. Front Microbiol. 2018;9:114. https://doi.org/10.3389/fmicb.2018.00114.

    Article  PubMed  PubMed Central  Google Scholar 

  78. Johnston PR, Rolff J. Host and symbiont jointly control gut microbiota during complete metamorphosis. Plos Pathog. 2015;11(11):e1005246. https://doi.org/10.1371/journal.ppat.1005246.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  79. Shao YQ, Spiteller D, Tang XS, Ping LY, Colesie C, Munchberg U, et al. Crystallization of alpha- and beta-carotene in the foregut of Spodoptera larvae feeding on a toxic food plant. Insect Biochem Molec. 2011;41(4):273–81. https://doi.org/10.1016/j.ibmb.2011.01.004.

    Article  CAS  Google Scholar 

  80. Shao YQ, Chen BS, Sun C, Ishida K, Hertweck C, Boland W. Symbiont-derived antimicrobials contribute to the control of the lepidopteran gut microbiota. Cell Chem Biol. 2017;24(1):66–75. https://doi.org/10.1016/j.chembiol.2016.11.015.

    Article  CAS  PubMed  Google Scholar 

  81. Vilanova C, Baixeras J, Latorre A, Porcar M. The generalist inside the specialist: gut bacterial communities of two insect species feeding on toxic plants are dominated by Enterococcus sp. Front Microbiol. 2016;7:1005. https://doi.org/10.3389/fmicb.2016.01005.

    Article  PubMed  PubMed Central  Google Scholar 

  82. Lilburn TG, Kim KS, Ostrom NE, Byzek KR, Leadbetter JR, Breznak JA. Nitrogen fixation by symbiotic and free-living spirochetes. Science. 2001;292(5526):2495–8. https://doi.org/10.1126/science.1060281.

    Article  CAS  PubMed  Google Scholar 

  83. Morales-Jiménez J, Zúñiga G, Ramírez-Saad HC, Hernández-Rodríguez C. Gut-associated bacteria throughout the life cycle of the bark beetle Dendroctonus rhizophagus Thomas and Bright (Curculionidae: Scolytinae) and their cellulolytic activities. Microb Ecol. 2012;64(1):268–78. https://doi.org/10.1007/s00248-011-9999-0.

    Article  PubMed  Google Scholar 

  84. Tang RX, Liu FY, Lan Y, Wang J, Wang L, Li J, et al. Transcriptomics and metagenomics of common cutworm (Spodoptera litura) and fall armyworm (Spodoptera frugiperda) demonstrate differences in detoxification and development. BMC Genomics. 2022;23(1):388. https://doi.org/10.1186/s12864-022-08613-6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

This study was supported by the Guizhou Provincial Science and Technology Projects ([2020]1Y105), the Guizhou Province Science and Technology Innovation Talent Team Project ([2021]004), and the Training Program of Guizhou University ([2019]32 and [2019]45).

Funding

This study was supported by the Supported by Guizhou Provincial Science and Technology Projects ([2020]1Y105), the Guizhou Province Science and Technology Innovation Talent Team Project ([2021]004), and the Training Program of Guizhou University ([2019]32 and [2019]45).

Author information

Authors and Affiliations

Authors

Contributions

LHW performed the experiments, analyzed the data, prepared figures and tables, authored drafts of the paper, and approved the final draft; CXH and TXL conceived the ideas, designed the experiments and selected the appropriate methodologies. All authors contributed to the manuscript and gave their final approval for the publication of the current manuscript.

Corresponding authors

Correspondence to Chao-xing Hu or Tong-xian Liu.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wu, Lh., Hu, Cx. & Liu, Tx. Metagenomic profiling of gut microbiota in Fall Armyworm (Spodoptera frugiperda) larvae fed on different host plants. BMC Microbiol 24, 337 (2024). https://doi.org/10.1186/s12866-024-03481-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12866-024-03481-4

Keywords