
==== Front
BMC Microbiol
BMC Microbiol
BMC Microbiology
1471-2180
BioMed Central London

39256682
3481
10.1186/s12866-024-03481-4
Research
Metagenomic profiling of gut microbiota in Fall Armyworm (Spodoptera frugiperda) larvae fed on different host plants
Wu Li-hong 123
Hu Chao-xing hcx333@foxmail.com

123
Liu Tong-xian tx.liu@gzu.edu.cn

123
1 https://ror.org/02wmsc916 grid.443382.a 0000 0004 1804 268X Institute of Entomology, Guizhou University, Guiyang, China
2 Guizhou Provincial Key Laboratory for Agricultural Pest Management of the Mountainous Region, Scientific Observing and Experimental Station of Crop Pests in Guiyang, Ministry of Agriculture and Rural Affairs of the People’s Republic of China, Guiyang, China
3 https://ror.org/02wmsc916 grid.443382.a 0000 0004 1804 268X Institute of Plant Health and Medicine, Guizhou University, Guiyang, 550025 China
11 9 2024
11 9 2024
2024
24 33731 3 2024
27 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ 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/.
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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-024-03481-4.

Keywords

Antibiotic resistance genes
Co-evolutionary relationship
Gut microbiota
Host plants
Metagenomic sequencing
issue-copyright-statement© BioMed Central Ltd., part of Springer Nature 2024
==== Body
pmcBackground

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–3]. Its rapid reproduction in local maize crops has led to population outbreaks, posing severe challenges to crop production [4–6].

The microbial community residing within the gut of insects encompasses a diverse array of organisms, including viruses, archaea, bacteria, fungi, and protozoa [7–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–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–30] and amplicon-oriented DNA sequencing [29, 31–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

Host plant’s abbreviation	Host plant’s common name	Scientific name	Varieties	Family	Stage of plant used in feeding the larvae	
CA	pepper	Capsicum annuum L	Guiyan 13	Solanaceous	3–4-true-leaf stages	
SL	tomato	Solanum lycopersicum Mill	Jingang	Solanaceous	3–4-true-leaf stages	
ZM	maize	Zea mays L	Qiannuo 168	Poaceae	3–4-true-leaf stages	
SM	Sorghum	Sorghum bicolor (L.) Moench	Hongyingzi No.1	Poaceae	3–4-true-leaf stages	

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

Sample ID	Insert Size (bp)	Seq Strategy	Raw Reads	Raw Base (GB)	%GC	Raw Q20 (%)	Raw Q30 (%)	Clean Reads	Cleaned (%)	Clean Q20 (%)	Clean Q30 (%)	
AD1	350	(150:150)	22,127,797	6.64	36	96.97	91.84	3,469,319	15.68	98.27	93.82	
AD2	350	(150:150)	20,968,284	6.29	36	97.20	92.39	3,846,087	18.34	98.44	94.31	
AD3	350	(150:150)	21,400,017	6.42	36	97.16	92.31	3,347,801	15.64	98.40	94.19	
CA1	350	(150:150)	22,554,767	6.77	39	97.53	93.14	12,388,909	54.93	98.67	95.02	
CA2	350	(150:150)	22,245,738	6.67	38	97.27	92.56	8,454,381	38	98.55	94.65	
CA3	350	(150:150)	21,966,559	6.59	38	97.07	92.15	7,006,173	31.89	98.41	94.27	
SL1	350	(150:150)	22,194,540	6.66	36	97.45	93.01	4,866,414	21.93	98.66	94.93	
SL2	350	(150:150)	21,950,213	6.59	37	97.51	93.14	5,397,615	24.59	98.70	95.03	
SL3	350	(150:150)	20,769,577	6.23	37	97.51	93.13	6,512,698	31.36	98.74	95.18	
ZM1	350	(150:150)	19,805,025	5.94	36	96.97	91.90	3,018,373	15.24	98.27	93.84	
ZM2	350	(150:150)	20,793,787	6.24	37	97.46	93.04	3,599,988	17.31	98.69	95.00	
ZM3	350	(150:150)	23,226,950	6.97	36	97.11	92.21	3,417,473	14.71	98.36	94.10	
SM1	350	(150:150)	22,350,535	6.71	45	97.59	93.57	14,229,565	63.67	98.86	95.72	
SM2	350	(150:150)	22,243,831	6.67	38	97.41	93.19	7,314,913	32.89	98.81	95.48	
SM3	350	(150:150)	23,111,706	6.93	44	97.79	93.89	13,500,446	58.41	98.90	95.80	
Sample ID represents the sample name indicating different diets: AD for artificial diet-fed, CA for pepper-fed, SL for tomato-fed, SM for sorghum-fed, and ZM for maize-fed. Insert Size is the length selected during library construction and gel cutting. Suitable Insert Size can avoid sequencing adaptor contamination. The sequencing strategy is generally paired-end, with each being 150 bp. Raw reads are the number of original reads that were sequenced using Illumina NovaSeq. Raw Base (GB) refers to the number of raw reads in GB, which is the total number of bases in the sequencing raw data, calculated by multiplying the number of raw reads by the sequencing length. GC (%) is the percentage of G and C type bases in the total base. Clean reads refer to the high-quality reads that remain after quality control and splicing. Cleaned(%) refers to the percentage of remaining sequences after filtering compared to the raw reads. Cleaned Q20 (%) refers to the percentage of bases with a quality score higher than 20 after filtering. Cleaned Q30 (%) refers to the percentage of bases with a quality score higher than 30 after filtering

Fig. 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

taxonomy	AD1	AD2	AD3	CA1	CA2	CA3	SL1	SL2	SL3	SM1	SM2	SM3	ZM1	ZM2	ZM3	Mean	
Archaea	0.00%	0.00%	0.00%	0.00%	0.00%	0.00%	0.00%	0.00%	0.00%	0.00%	0.00%	0.00%	0.00%	0.00%	0.00%	0.00%	
Bacteria	99.97%	99.85%	99.98%	99.80%	99.83%	99.92%	99.99%	99.96%	99.70%	99.74%	99.79%	98.35%	99.85%	99.46%	99.71%	99.73%	
Phages	0.00%	0.00%	0.00%	0.00%	0.00%	0.00%	0.00%	0.00%	0.00%	0.04%	0.01%	0.01%	0.00%	0.00%	0.00%	0.00%	
Plants and Fungi	0.03%	0.15%	0.02%	0.20%	0.17%	0.08%	0.01%	0.04%	0.29%	0.22%	0.20%	1.64%	0.15%	0.54%	0.29%	0.27%	
The abbreviations in the table represent different diets, AD artificial diet-fed, CA pepper-fed, SL tomato-fed, SM sorghum-fed, ZM maize-fed

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 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

Sample ID	Firmicutes	Proteobacteria	Actinobacteria	Microsporidia	Ascomycota	
AD	92.98 ± 5.00 a	5.77 ± 5.39 b	1.18 ± 0.96 ab	0.03 ± 0.02 a	0.03 ± 0.02 a	
CA	97.83 ± 0.33 a	1.97 ± 0.29 b	0.04 ± 0.01 b	0.10 ± 0.02 a	0.04 ± 0.01 a	
SL	98.73 ± 0.37 a	0.98 ± 0.21 b	0.07 ± 0.01 b	0.09 ± 0.07 a	0.04 ± 0.01 a	
SM	37.76 ± 11.97 b	59.64 ± 11.16 a	1.76 ± 0.56 a	0.25 ± 0.21 a	0.27 ± 0.18 a	
ZM	94.58 ± 2.29 a	4.90 ± 2.16 b	0.19 ± 0.12 ab	0.23 ± 0.11 a	0.06 ± 0.02 a	
Data are mean ± SE. Sample ID represents the sample name indicating different diets: AD artificial diet-fed, CA pepper-fed, SL tomato-fed, SM sorghum-fed, ZM maize-fed

The different lowercase letter indicates that the significant differences in relative abundance in the same column in the mean values. The same lowercase letter indicates that there is no significant difference between the groups, and different letters indicate that the differences between the groups is significant (one-way ANOVA, Tukey post-hoc test, P < 0.05)

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

Sample ID	Enterococcaceae	Lactobacillaceae	Enterobacteriaceae	Brucellaceae	Rhizobiaceae	
AD	59.99 ± 11.98 b	26.42 ± 6.76 a	5.69 ± 5.34 ab	0.01 ± 0.01 b	0.00 ± 0.00 b	
CA	97.77 ± 0.32 a	0.00 ± 0.00 b	1.44 ± 0.13 b	0.43 ± 0.28 b	0.00 ± 0.00 b	
SL	98.64 ± 0.36 a	0.00 ± 0.00 b	0.91 ± 0.20 b	0.03 ± 0.01 b	0.00 ± 0.00 b	
SM	37.51 ± 11.96 b	0.00 ± 0.00 b	12.61 ± 2.54 a	21.04 ± 2.13 a	16.73 ± 8.22 a	
ZM	94.52 ± 2.30 a	0.00 ± 0.00 b	3.58 ± 1.61 ab	1.20 ± 0.54 b	0.04 ± 0.02 b	
Note: Data are mean ± SE. Sample ID represents the sample name indicating different diets: AD for artificial diet-fed, CA for pepper-fed, SL for tomato-fed, SM for sorghum-fed, and ZM for maize-fed. The different lowercase letter indicates that the significant differences in relative abundance in the same column in the mean values. The same lowercase letter indicates that there is no significant difference between the groups, and different letters indicate that the differences between the groups is significant (one-way ANOVA, Tukey post-hoc test, P < 0.05)

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 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 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 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 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 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 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–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–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.

Supplementary Information

Supplementary Material 1.

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

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).

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.

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).

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.

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.

Publisher’s Note

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

1. Sparks AN A review of the biology of the fall armyworm Fla Entomo 1979 62 2 82 87 10.2307/3494083
Sparks AN. A review of the biology of the fall armyworm. Fla Entomo. 1979;62(2):82–7. 10.2307/3494083.10.2307/3494083
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 738 10.1093/aesa/73.6.722
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. 10.1093/aesa/73.6.722.10.1093/aesa/73.6.722
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 526 10.1603/0022-0493-99.2.519 16686155
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. 10.1603/0022-0493-99.2.519.16686155 10.1603/0022-0493-99.2.519
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 10.1371/journal.pone.0165632 27788251
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. 10.1371/journal.pone.0165632.27788251 10.1371/journal.pone.0165632
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 10.1038/s41598-017-04238-y
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. 10.1038/s41598-017-04238-y.10.1038/s41598-017-04238-y
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 15 10.1016/j.cropro.2017.11.003
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. 10.1016/j.cropro.2017.11.003.10.1016/j.cropro.2017.11.003
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 170 10.1111/eea.12768
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. 10.1111/eea.12768.10.1111/eea.12768
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 10.3389/fmicb.2018.00556 29636736
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. 10.3389/fmicb.2018.00556.29636736 10.3389/fmicb.2018.00556
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 10.1093/icb/ics089 22659204
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. 10.1093/icb/ics089.22659204 10.1093/icb/ics089
10. Engel P Moran NA The gut microbiota of insects-diversity in structure and function FEMS Microbiol Rev 2013 37 5 699 735 10.1111/1574-6976.12025 23692388
Engel P, Moran NA. The gut microbiota of insects-diversity in structure and function. FEMS Microbiol Rev. 2013;37(5):699–735. 10.1111/1574-6976.12025.23692388 10.1111/1574-6976.12025
11. Fraune S Bosch TC Why bacteria matter in animal development and evolution BioEssays 2010 32 7 571 580 10.1002/bies.200900192 20544735
Fraune S, Bosch TC. Why bacteria matter in animal development and evolution. BioEssays. 2010;32(7):571–80. 10.1002/bies.200900192.20544735 10.1002/bies.200900192
12. Janson EM Stireman JO 3rd Singer MS Abbot P Phytophagous insect-microbe mutualisms and adaptive evolutionary diversification Evolution 2008 62 5 997 1012 10.1111/j.1558-5646.2008.00348.x 18298649
Janson EM, Stireman JO 3rd, Singer MS, Abbot P. Phytophagous insect-microbe mutualisms and adaptive evolutionary diversification. Evolution. 2008;62(5):997–1012. 10.1111/j.1558-5646.2008.00348.x.18298649 10.1111/j.1558-5646.2008.00348.x
13. Hansen AK Moran NA The impact of microbial symbionts on host plant utilization by herbivorous insects Mol Ecol 2014 23 6 1473 1496 10.1111/mec.12421 23952067
Hansen AK, Moran NA. The impact of microbial symbionts on host plant utilization by herbivorous insects. Mol Ecol. 2014;23(6):1473–96. 10.1111/mec.12421.23952067 10.1111/mec.12421
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 10.1038/s42003-022-03947-z 36109578
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. 10.1038/s42003-022-03947-z.36109578 10.1038/s42003-022-03947-z
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 10.3390/agronomy13020304
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. 10.3390/agronomy13020304.10.3390/agronomy13020304
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 10.3389/fmicb.2022.824224 35479615
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. 10.3389/fmicb.2022.824224.35479615 10.3389/fmicb.2022.824224
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 10.3390/insects13080694 36005319
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. 10.3390/insects13080694.36005319 10.3390/insects13080694
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 10.3390/biology11091347 36138826
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. 10.3390/biology11091347.36138826 10.3390/biology11091347
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 10.3390/ijms22136843 34202141
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. 10.3390/ijms22136843.34202141 10.3390/ijms22136843
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 10.1007/s10340-023-01682-w
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. 10.1007/s10340-023-01682-w.10.1007/s10340-023-01682-w
21. Wang YP Liu X Yi CY Chen XY Liu CH Zhang CC 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 10.3390/genes14020321 36833248
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. 10.3390/genes14020321.36833248 10.3390/genes14020321
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 10.1007/s13199-023-00907-x
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. 10.1007/s13199-023-00907-x.10.1007/s13199-023-00907-x
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 956 10.1111/jen.13022
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. 10.1111/jen.13022.10.1111/jen.13022
24. Lǚ DB Liu XY Dong YL Yan ZZ Zhang X Wang P Comparison of gut bacterial communities of fall armyworm (Spodoptera frugiperda) reared on different host plants Int J Mol Sci 2021 22 20 11266 10.3390/ijms222011266 34681926
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. 10.3390/ijms222011266.34681926 10.3390/ijms222011266
25. Acevedo FE Peiffer M Tan CW Stanley BA Stanley A Wang J Jones AG Hoover K Rosa C Luthe D Fall armyworm-associated gut bacteria modulate plant defense responses Mol Plant Microbe Interact 2017 30 2 127 137 10.1094/MPMI-11-16-0240-R 28027025
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. 10.1094/MPMI-11-16-0240-R.28027025 10.1094/MPMI-11-16-0240-R
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 10.1371/journal.pone.0174754 28358907
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. 10.1371/journal.pone.0174754.28358907 10.1371/journal.pone.0174754
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 15996 10.1073/pnas.1908748116
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. 10.1073/pnas.1908748116.10.1073/pnas.1908748116
28. Mason CJ St Clair A Peiffer M Gomez E Jones AG Felton GW Diet influences proliferation and stability of gut bacterial populations in herbivorous lepidopteran larvae PLoS ONE 2020 15 3 e0229848 10.1371/journal.pone.0229848 32168341
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. 10.1371/journal.pone.0229848.32168341 10.1371/journal.pone.0229848
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 851 10.1007/s10340-020-01202-0
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. 10.1007/s10340-020-01202-0.10.1007/s10340-020-01202-0
30. Wang XL Sun SJ Yang XL Cheng J Wei HS Li Z Variability of gut microbiota across the life cycle of Grapholita molesta (Lepidoptera: Tortricidae) Front Microbiol 2020 11 1366 10.3389/fmicb.2020.01366 32714300
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. 10.3389/fmicb.2020.01366.32714300 10.3389/fmicb.2020.01366
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 10.1038/s41598-019-39163-9
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. 10.1038/s41598-019-39163-9.10.1038/s41598-019-39163-9
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 363 10.1111/jen.12742
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. 10.1111/jen.12742.10.1111/jen.12742
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 10.7717/peerj.8701 32185109
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. 10.7717/peerj.8701.32185109 10.7717/peerj.8701
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 776 10.1111/jen.12821
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. 10.1111/jen.12821.10.1111/jen.12821
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 10.1038/s41598-021-83497-2
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. 10.1038/s41598-021-83497-2.10.1038/s41598-021-83497-2
36. Rausch P Rühlemann M Hermes BM Doms S Dagan T Dierking K Comparative analysis of amplicon and metagenomic sequencing methods reveals key features in the evolution of animal metaorganisms Microbiome 2019 7 1 133 10.1186/s40168-019-0743-1 31521200
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. 10.1186/s40168-019-0743-1.31521200 10.1186/s40168-019-0743-1
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 10.1016/j.envint.108655 38626494
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. 10.1016/j.envint.108655.38626494 10.1016/j.envint.108655
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 10.1016/j.scitotenv.153033 35026253
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. 10.1016/j.scitotenv.153033.35026253 10.1016/j.scitotenv.153033
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. 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 10.3389/fmicb.2017.00663 28491055
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. 10.3389/fmicb.2017.00663.28491055 10.3389/fmicb.2017.00663
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 10.3389/fmicb.2022.909863 35668757
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. 10.3389/fmicb.2022.909863.35668757 10.3389/fmicb.2022.909863
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 10.3390/insects12030203 33673724
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. 10.3390/insects12030203.33673724 10.3390/insects12030203
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 10.3390/ijms21176283 32872681
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. 10.3390/ijms21176283.32872681 10.3390/ijms21176283
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 1362 10.1002/ps.5646 31605420
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. 10.1002/ps.5646.31605420 10.1002/ps.5646
45. Bolger AM Lohse M Usadel B Trimmomatic: a flexible trimmer for Illumina sequence data Bioinformatics 2014 30 15 2114 2120. 82 10.1093/bioinformatics/btu170 24695404
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114-2120. 82. 10.1093/bioinformatics/btu170.24695404 10.1093/bioinformatics/btu170
46. Langmead B Salzberg SL Fast gapped-read alignment with bowtie 2 Nat Methods 2012 9 D1 357 359 10.1038/nmeth.1923 22388286
Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9(D1):357–9. 10.1038/nmeth.1923.22388286 10.1038/nmeth.1923
47. Wood DE Lu J Langmead B Improved metagenomic analysis with Kraken2 Genome Biol 2019 20 1 257 10.1186/s13059-019-1891-0 31779668
Wood DE, Lu J, Langmead B. Improved metagenomic analysis with Kraken2. Genome Biol. 2019;20(1):257. 10.1186/s13059-019-1891-0.31779668 10.1186/s13059-019-1891-0
48. Lu J Breitwieser FP Thielen P Salzberg SL Bracken: Estimating species abundance in metagenomics data Peerj Comput Sci 2017 3 e104 10.7717/peerj-cs.104
Lu J, Breitwieser FP, Thielen P, Salzberg SL. Bracken: Estimating species abundance in metagenomics data. Peerj Comput Sci. 2017;3:e104. 10.7717/peerj-cs.104.10.7717/peerj-cs.104
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 A practical guide to amplicon and metagenomic analysis of microbiome data Protein Cell 2021 12 5 315 330 10.1007/s13238-020-00724-8 32394199
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. 10.1007/s13238-020-00724-8.32394199 10.1007/s13238-020-00724-8
51. Segata N Izard J Waldron L Gevers D Miropolsky L Garrett WS Metagenomic biomarker discovery and explanation Genome Biol 2011 12 6 R60 10.1186/gb-2011-12-6-r60 21702898
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. 10.1186/gb-2011-12-6-r60.21702898 10.1186/gb-2011-12-6-r60
52. Buchfink B Xie C Huson DH Fast and sensitive protein alignment using DIAMOND Nat Methods 2015 12 1 59 60 10.1038/nmeth.3176 25402007
Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2015;12(1):59–60. 10.1038/nmeth.3176.25402007 10.1038/nmeth.3176
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 968 10.1038/s41592-018-0176-y 30377376
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. 10.1038/s41592-018-0176-y.30377376 10.1038/s41592-018-0176-y
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 D361 10.1093/nar/gkw1092 27899662
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. 10.1093/nar/gkw1092.27899662 10.1093/nar/gkw1092
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 2262 10.1038/s41396-018-0174-1 29895989
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. 10.1038/s41396-018-0174-1.29895989 10.1038/s41396-018-0174-1
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. 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 18556 10.1021/acs.jafc.3c05506 37963218
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. 10.1021/acs.jafc.3c05506.37963218 10.1021/acs.jafc.3c05506
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 10.3390/insects13090762 36135463
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. 10.3390/insects13090762.36135463 10.3390/insects13090762
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 10.3390/insects13040373 35447815
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. 10.3390/insects13040373.35447815 10.3390/insects13040373
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 763 10.1016/S2095-3119(20)63476-1
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. 10.1016/S2095-3119(20)63476-1.10.1016/S2095-3119(20)63476-1
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 2605 10.1007/s00248-023-02264-0 37433981
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. 10.1007/s00248-023-02264-0.37433981 10.1007/s00248-023-02264-0
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 10.3390/insects14040358 37103173
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. 10.3390/insects14040358.37103173 10.3390/insects14040358
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 10.4161/gmic.19897 22572875
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. 10.4161/gmic.19897.22572875 10.4161/gmic.19897
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 3800 10.1111/mec.13730 27297628
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. 10.1111/mec.13730.27297628 10.1111/mec.13730
65. Huang ZD Zhu L Lv J Pu ZX Zhang LP Chen GQ Dietary effects on biological parameters and gut microbiota of Harmonia axyridis Front Microbiol 2022 12 818787 10.3389/fmicb.2021.818787 35154044
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. 10.3389/fmicb.2021.818787.35154044 10.3389/fmicb.2021.818787
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 7006 10.1128/AEM.05609-11
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. 10.1128/AEM.05609-11.10.1128/AEM.05609-11
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 10.1016/j.jinsphys.2006.02.007 16600286
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. 10.1016/j.jinsphys.2006.02.007.16600286 10.1016/j.jinsphys.2006.02.007
68. Warnecke F Luginbühl P Ivanova N Ghassemian M Richardson TH Stege JT Metagenomic and functional analysis of hindgut microbiota of a wood-feeding higher termite Nature 2007 450 7169 560 565 10.1038/nature06269 18033299
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. 10.1038/nature06269.18033299 10.1038/nature06269
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 10.1093/jisesa/iead009 36916277
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. 10.1093/jisesa/iead009.36916277 10.1093/jisesa/iead009
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 10.1007/s10123-022-00303-3 36484909
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. 10.1007/s10123-022-00303-3.36484909 10.1007/s10123-022-00303-3
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 10.1186/s40168-022-01456-z 36670494
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. 10.1186/s40168-022-01456-z.36670494 10.1186/s40168-022-01456-z
72. Shin SC Kim SH You H Kim B Kim AC Lee KA Drosophila microbiome modulates host developmental and metabolic homeostasis via insulin signaling Science 2011 334 6056 670 674 10.1126/science.1212782 22053049
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. 10.1126/science.1212782.22053049 10.1126/science.1212782
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 907 10.1242/jeb.151522 28062579
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. 10.1242/jeb.151522.28062579 10.1242/jeb.151522
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 414 10.1016/j.cmet.2011.07.012 21907145
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. 10.1016/j.cmet.2011.07.012.21907145 10.1016/j.cmet.2011.07.012
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 406 10.1038/s41586-018-0634-9 30356215
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. 10.1038/s41586-018-0634-9.30356215 10.1038/s41586-018-0634-9
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 629 10.1603/0046-225X-35.3.625
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. 10.1603/0046-225X-35.3.625.10.1603/0046-225X-35.3.625
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 10.3389/fmicb.2018.00114 29449838
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. 10.3389/fmicb.2018.00114.29449838 10.3389/fmicb.2018.00114
78. Johnston PR Rolff J Host and symbiont jointly control gut microbiota during complete metamorphosis Plos Pathog 2015 11 11 e1005246 10.1371/journal.ppat.1005246 26544881
Johnston PR, Rolff J. Host and symbiont jointly control gut microbiota during complete metamorphosis. Plos Pathog. 2015;11(11):e1005246. 10.1371/journal.ppat.1005246.26544881 10.1371/journal.ppat.1005246
79. Shao YQ Spiteller D Tang XS Ping LY Colesie C Munchberg U 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 281 10.1016/j.ibmb.2011.01.004
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. 10.1016/j.ibmb.2011.01.004.10.1016/j.ibmb.2011.01.004
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 10.1016/j.chembiol.2016.11.015 28107652
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. 10.1016/j.chembiol.2016.11.015.28107652 10.1016/j.chembiol.2016.11.015
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 10.3389/fmicb.2016.01005 27446044
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. 10.3389/fmicb.2016.01005.27446044 10.3389/fmicb.2016.01005
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 2498 10.1126/science.1060281 11431569
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. 10.1126/science.1060281.11431569 10.1126/science.1060281
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 278 10.1007/s00248-011-9999-0 22234511
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. 10.1007/s00248-011-9999-0.22234511 10.1007/s00248-011-9999-0
84. Tang RX Liu FY Lan Y Wang J Wang L Li J 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 10.1186/s12864-022-08613-6 35596140
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. 10.1186/s12864-022-08613-6.35596140 10.1186/s12864-022-08613-6
