Lei Qian
Yanhui Wang2
Huaijian Liao- 1Institute of Leisure Agriculture, Jiangsu Academy of Agricultural Sciences, Nanjing, China
- 2College of Plant Protection, Nanjing Agricultural University, Nanjing, China
- 3College of Plant Protection, Yangzhou University, Yangzhou, China
- 4College of Biotechnology, Jiangsu University of Science and Technology, Zhenjiang, China
Introduction: The soybean hawkmoth, Clanis bilineata tsingtauica, is an edible insect that possesses high nutritional, medicinal and economic value. It has developed into a characteristic agricultural industry in China.
Methods: The dominant gut bacterium in diapause larvae of soybean hawkmoths was identified by metagenomics, and the effect of diapause time on gut microbiome composition, diversity and function was investigated.
Results: Enterococcus and Enterobacter were measured to be the dominant genera, with Enterococcus casseliflavus and Enterococcus pernyi being the dominant species. Compared to the controls, the relative abundance of E. casseliflavus and E. pernyi on day 14 was lower by 54.51 and 42.45%, respectively. However, the species richness (including the index of Chao and ACE) of gut microbiota increased on day 28 compared to controls. The gene function was mainly focused on carbohydrate and amino acid metabolism. Metabolic pathways annotated for amino acids on day 14 increased by 9.83% compared to controls. It is speculated that diapause soybean hawkmoths may up-regulate amino acid metabolism by reducing E. casseliflavus abundance to maintain their nutritional balance. Additionally, tetracycline, chloromycetin and ampicillin were screened as the top three antibiotics against E. casseliflavus.
Discussion: This study not only extends our knowledge of gut microbiome in soybean hawkmoths at the species level, but also provides an initial investigation of gene functionality in interaction with insect hosts.
1 Introduction
In China, there is a unique edible insect species Clanis bilineata tsingtauica, soybean hawkmoth, which belongs to the Sphingidae family of Lepidoptera (Zhu and Wang, 1997). The larvae of soybean hawkmoth have high nutritional, medicinal, and economic value with abundant protein, vitamins, amino acids (especially essential amino acids (EAAs)), unsaturated fatty acids, microelements, chitin, etc. (Tian and Zhang, 2012). These nutrients can help people promote brain development, prevent the cell deterioration and maintain the endocrine balance (Qian et al., 2022). Currently, in the rural areas of China, especially in the north of Jiangsu Province, the artificial breeding of soybean hawkmoth has become a characteristic agricultural industry with an output of about 30,000 tons and 4.5 billion yuan every year (Guo et al., 2021). Under natural conditions, one generation of the soybean hawkmoth occurs each year in China. The old fifth-instar larvae burrow into the soil and keep still in early September, entering the diapause stage (Guo et al., 2021; Qian et al., 2023a).
Many insect species adaptively choose diapause, similar to hibernation or dormancy in mammals, to increase survival under adverse environmental conditions [e.g., low temperature, lack of food and water (Dittmer and Brucker, 2021)]. It is critical to provide adequate nutrition for diapause insects, and many species build up nutrient reserves in advance. Despite entering a stage of developmental arrest, diapause insects are able to maintain a stable body weight and nutritional status (Santana et al., 2017). However, the nutrient composition of some insects also varies with the intensity of diapause. Examples include lipids in diapause Arimania comaroffi, and glycerol, trehalose and glycogen in diapause Aphidius gifuensis (Bemani et al., 2012). During diapause, insects can self-regulate the composition and levels of nutrients, such as protein, fatty acids, EAAs and carbohydrates (Hahn and Denlinger, 2011). In addition to nutrient storage, insects also reduce their metabolism in order to meet energy requirements for diapause and to provide energy sources for post-diapause growth and development (Poelchau et al., 2013; Dittmer and Brucker, 2021). Despite the growing interest of entomologists in the regulatory mechanisms of insect diapause, the interaction between gut microbes and diapause insects, especially in aspects of nutrient and metabolic regulation, is still poorly understood (Engel and Moran, 2013; Douglas, 2015).
In numerous insect species, gut microbes are involved in nutrient metabolism, digestion, gut homeostasis and endocrine signaling to promote their development and growth, as noted by Zheng et al. (2017) and van Opstal and Bordenstein (2019). For instance, gut microbiota can provide the host insects with a steady supply of essential nutrients, including EAAs, carbohydrates and lipids, to sustain their metabolic processes (Fraune and Bosch, 2010). The gut microbiota assists in fixing atmospheric nitrogen and promoting amino acid metabolism, leading to a partial shift in the host’s amino acid composition pattern (Storelli et al., 2011). In addition, gut bacteria contribute significantly to lipid and protein digestion in Tenebrio molitor and Anticarsia gemmatalis (Genta et al., 2006; Visotto et al., 2009). They also exhibit enzymatic action on starch, lipids, pectin and esters in Helicoverpa armigera, Bombyx mori, Plutella xylostella and Antheraea assamensis (Anand et al., 2010; Gandotra et al., 2018). In Riptortus pedestris, the presence of Burkholderia in the gut results in higher levels of juvenile hormone (JH), leading to increased fecundity and storage proteins (Lee et al., 2019). Moreover, it has been demonstrated that Klebsiella oxytoca can boost the developmental stage of Drosophila suzukii by activating the glycolysis pathway (Gao et al., 2023). The gut microbiota in Drosophila melanogaster modulates the insulin signaling pathway, a key regulator of insect growth and nutrient homeostasis, to regulate nutrient allocation in the host (Shin et al., 2011).
In our previous study, 16S rRNA sequencing was conducted to analyze the gut microbes of diapause soybean hawkmoths. The results showed that Proteobacteria and Firmicutes were the most prominent taxa at the phylum level in the larval gut microbiota, with Enterococcus and Stenotrophomonas being the dominant genera. In addition, non-target metabolomics analysis on the intestinal contents showed that the most enriched pathways was mainly focused on amino acid metabolism (Qian et al., 2023b). Our findings suggest that the dominant gut bacteria may be involved in the nutritional physiology of insect hosts. However, the specific dominant strains of diapause soybean hawkmoths are still unidentified, not to mention their regulatory mechanisms for nutrient biosynthesis and metabolism in the host. This study identified the predominant strains within the larval gut via the metagenomic sequencing, and investigated the impact of diapause time on the diversity, species composition and gene functions of the gut microbiome. This study reports the first discovery of the dominant gut bacteria of soybean hawkmoths, as well as their influence on the nutrition of the hosts. Additionally, antimicrobial susceptibility tests were conducted for the most effective antibiotic to target the dominant bacterium in diapause larvae. Understanding how dominant gut bacteria regulate nutrient synthesis and metabolism in the soybean hawkmoth will help to explore and utilize the nutritional substances of edible insects, promote the development of the soybean hawkmoth industry in China’s rural areas, and further advance the implementation of the rural revitalization strategy.
2 Materials and methods
2.1 Insects and experimental design
In this study, the population of soybean hawkmoth C. bilineata tsingtauica was sourced from Yuntai Farm of Jiangsu Agricultural Reclamation Group, situated in Lian Yun-gang City in Jiangsu Province, China. The soybean hawkmoths were reared on soybean plant leaves at 27°C/16 h light, and 25°C/8 h dark with 70% relative humidity (RH). Once the mature fifth instar larvae start to burrow and cease movement in the soil, they are deemed to have entered the diapause phase. The day when mature larvae entered the soil was recorded as day 0 (considered as the control). Based on our previous study utilizing 16S rRNA sequencing analysis (Qian et al., 2023b), the diapause larvae of soybean hawkmoth at day 0, 14 and 28 were selected for the following metagenomic analysis. Only diapause larvae showing signs of vitality or movement upon stimulation were considered viable. Live larvae were collected for subsequent DNA extraction. The efficient antibiotics were evaluated against the dominant strains identified through metagenomic analysis.
2.2 Midgut contents collection
Diapause larvae were collected at various time points (day 0 (D0), day 14 (D14) and day 28 (D28)) and dissected afterward. After cleansing the larvae in a 75% ethanol solution for 2 min, they were rinsed thrice with sterile water to eliminate surface contaminants, and then individually dissected using aseptic forceps and scissors. The midguts were aseptically sampled to obtain intestinal contents, rapidly frozen in liquid nitrogen, and stored at −80°C following the protocol described by Chen et al. (2023). Metagenomic analysis was performed on the intestinal contents collected from diapause larvae. One sample contained the intestinal contents from 20 larvae, with three experimental repeats.
2.3 Metagenome DNA extraction and shotgun sequencing
The total microbial genomic DNA samples were extracted using the OMEGA Mag-Bind Soil DNA Kit (M5635-02; Omega Bio-Tek, Norcross, GA, United States) and were kept stored at −20°C for subsequent evaluation. The quantity and quality of extracted DNA were analyzed with the Qubit™ 4 Fluorometer, with WiFi: Q33238 (Qubit™ Assay Tubes: Q32856; Qubit™ 1X dsDNA HS Assay Kit: Q33231; Invitrogen, USA) and agarose gel electrophoresis. The microbial DNA was used to establish metagenomics shotgun sequencing libraries with 400 bp inserts at Personal Biotechnology Co., Ltd., located in Shanghai, China, using the Illumina TruSeq Nano DNA LT Library Preparation Kit (Kuang et al., 2023). The Illumina NovaSeq platform (Illumina, United States) was employed for sequencing the libraries using a PE150 strategy.
2.4 Metagenomic analysis
Clean genomic assembly data were acquired from raw data using FASTP software (v18.0) (Chen S. et al., 2018). Megahit (v1.1.2) and MetaGeneMark (v3.38) were used for genetic assembly and prediction (Zhu et al., 2010; Li et al., 2015). To cluster the CDS sequences of all samples, mmseqs2 was employed in “easy-cluster” mode, with a protein sequence identity threshold of 0.90 and a coverage of 90% for the shorter contig (Steinegger and Söding, 2017). To analyze the gene abundance, high-quality reads from each sample were aligned to predicted gene sequences using Salmon in ‘quasi-mapping-based’ mode (Patro et al., 2015). Abundance values in the metagenomes were normalized using copies per kilobase per million mapped reads (CPM). Non-redundant gene function was determined by annotating with mmseqs2 in the search mode against the KEGG, EggNOG and CAZy protein databases, respectively. EggNOG and GO were obtained using EggNOG-mapper (v2) (Cantalapiedra et al., 2021). The GO ontology was procured via map2slim from www.metacpan.org, and KO was acquired through KOBAS (Bu et al., 2021).
2.5 Antimicrobial susceptibility tests
Enterococcus casseliflavus was identified as the dominant bacterium in diapause larvae of C. bilineata tsingtauica through metagenomic analysis. To explore the nutritional function of E. casseliflavus affecting on larvae, the antimicrobial susceptibility tests were carried out to screen the most effective antibiotic against E. casseliflavus. The following antibiotics were used in the experiment: ampicillin (AMP), vancomycin (VAN), ciprofloxacin (CIP), gentamicin (GM), erythromycin (EM), levofloxacin (LVX), chloromycetin (CHL) and tetracycline (TC) (Liang et al., 2022; Zhang and Ju, 2023). The strain of E. casseliflavus was obtained from China General Microbiological Culture Collection Center (CGMCC)1 and cultured according to the instructions. The solid media were prepared with 30 g trypticase soy broth, 15 g agar, and 1 L double distilled water, sterilized at 121°C for 15 min. Evenly applied the solution of E. casseliflavus onto the solid media, cultured at 37°C for 12 h. The colony count reached a period of logarithmic growth. Place the susceptibility paper bought from Shanghai Ebaiju Economic and Trade Co. LTD (Shanghai, China) at the center of the media. Incubate the media at 37°C for 24 h. Measure the bacteriostatic diameter in each medium using a vernier caliper, and carefully determine the clarity of the bacteriostatic zones. Each antibiotic treatment was repeated for six times.
2.6 Statistical analysis
Linear discriminant analysis effect size (LEfSe) was performed to detect differentially abundant taxa and functions across groups using the default parameters based on the taxonomic and functional profiles of non-redundant genes (Segata et al., 2011). Bray-Curtis distance metrics were utilized to perform beta diversity analysis, examining compositional and functional differences in microbial communities among samples. The results were presented through principal coordinate analysis (PCoA), non-metric multidimensional scaling (NMDS) and unweighted pair-group method with arithmetic means (UPGMA) hierarchical clustering (Bray and Curtis, 1957; Ramette, 2007). Antimicrobial susceptibility test results were reported as mean and standard errors of mean (Mean ± SEM). The diameters of inhibition zones treated with eight antibiotics were analyzed through the SPSS 25.0 software (Chicago, IL, United States), employing a one-way analysis of variance (ANOVA) test. Tukey’s test was used to determine statistically significant differences among the treatments at p < 0.05. All data were first analyzed by using Levene and Kolmogorov–Smirnov tests to determine the homogeneity of variance and normality, respectively. Transformed data was used for analysis if the raw data did not confirm the homogeneity.
3 Results
3.1 Quality analysis of metagenomic sequencing data from diapause larvae
Nine metagenomic DNA libraries were sequenced on the Illumina NovaSeq platform, constructed from the control group (D0) and the diapause groups (D14 and D28). As shown in Table 1, the number of total reads ranged from 41,811,872 to 50,374,154, with a base count ranging from 6,271,780,800 to 7,556,123,100. The GC (%) of the control and diapause groups ranged from 43.86 to 47.3% and from 42.94 to 49.59%, respectively (see Table 1). The Ns (%) presented in Table 1 indicate that fuzzy bases were found in both groups. Moreover, it is noteworthy that all Q20s (%) and Q30s (%) from both control and diapause groups exceeded 95 and 90%, respectively (refer to Table 1). The sequencing data met the prescribed criteria, as evidenced by statistical distributions for sequencing quality (Supplementary Figure S1A), mean mass (Supplementary Figure S1B), base species content (Supplementary Figure S1C), and GC content (Supplementary Figure S1D). Consequently, it is viable for further analysis.
Table 1. Sequencing data of intestinal microbial metagenomics in the dipause larvae of soybean hawkmoth Clanis bilineata tsingtauica.
3.2 Diapause time reduces the abundance of gut dominant bacteria Enterococcus
The predominant bacterial genera found in the midgut of diapause larvae of soybean hawkmoth were Enterococcus and Enterobacter, with the primary species being Enterococcus casseliflavus and Enterococcus pernyi (Figure 1 and Supplementary Figure S2). The relative abundance (RA) of Enterococcus significantly decreased by 65.74% on D14 compared to D0 controls (p < 0.05, F = 9.71, df = 2; Figure 1B). However, the RA of Enterobacter on D14 was significantly higher than the controls by 1031.10% (p < 0.05, F = 6.38, df = 2; Figure 1B). A similar trend was noticed in the RA of E. casseliflavus and E. pernyi over diapause time (Figure 1D). Specifically, both of E. casseliflavus and E. pernyi displayed a significantly lower RA on D14 in comparison to the control group (p < 0.05, F = 6.60, df = 2; p < 0.05, F = 6.81, df = 2; Figure 1D).
Figure 1. Effects of diapause time on the compositions of gut microbiota in the dipause larvae of soybean hawkmoth Clanis bilineata tsingtauica. The compositions (A) and the relative abundance (RA) (B) of Enterococcus and Enterobacter at the genus level. The compositions (C) and the RA (D) of Enterococcus casseliflavus and Enterococcus pernyi at the species level. Each histogram distinguishes taxa by color, and the ordinate represents the RA of each taxon. The longer the column, the higher the RA of the taxa. Different lowercase letters denote significant difference among diapause time by Tukey test at p < 0.05.
As shown in Figure 2, E. casseliflavus exhibited a positive association with Enterococcus faecium and Enterococcus gallinarum, and a negative association with Enterobacter cloacae and Enterobacter hormaechei. By contrast, E. pernyi displayed a positive correlation only with Burkholderia stabilis.
Figure 2. Map of the dominant microbial association network in the dipause larvae of soybean hawkmoth C. bilineata tsingtauica at species level. The dominant genera represented by nodes are marked with different colors. The connection between the nodes indicates a correlation between two species. Red lines indicate a positive correlation and green lines indicate a negative correlation.
3.3 Diapause time increases the diversity of larval gut microbiome
There was a significant increase in Chao and ACE (abundance-based coverage estimator) on D28 compared to D0 (p < 0.05, F = 10.03, df = 2; p < 0.05, F = 10.13, df = 2; Figure 3A). This suggests that the species richness of gut microbiota increased in diapause larvae on D28. PCoA demonstrated a clear separation of gut microbiota profiles between D14 and D0, D28 (Figure 3B). The cluster heatmap of the top 50 differential species, the LEfSe analysis and the UPGMA hierarchical clustering analysis all provided supportive evidence for PCoA (Supplementary Figure S3). A total of 99 common microbial species among groups were revealed in the Venn diagram presented in Figure 3C.
Figure 3. Effects of diapause time on the diversity and species difference of gut microbiota in the diapause larvae of C. bilineata tsingtauica. α-diversity (A), PCoA analysis (B) and Venn figure (C) of gut microbiota. Different lowercase letters denote significant differences among diapause time by Tukey test at p < 0.05.
3.4 Diapause time up-regulates larval amino acid metabolic pathways
Larval gut microbiota function was mainly focused on metabolism (see Figure 4A). The top five pathways included carbohydrate metabolism, amino acid metabolism, energy metabolism, metabolism of cofactors and vitamins, and nucleotide metabolism (Figure 4A). In addition, the annotated amino acid metabolic pathways were significantly higher in diapause larvae on D14 than in controls (p < 0.05, F = 11.98, df = 2; Figure 4B).
Figure 4. Analysis of annotated KEGG metabolic pathways in the diapause larvae of C. bilineata tsingtauica (A). Effects of diapause time on the KEGG metabolic pathways in the dipause larvae (B). The horizontal coordinates represent the number of proteins annotated to the corresponding pathways, and the vertical coordinates correspond to pathways at KEGG level 2. The classification of pathways at KEGG level 1 is shown on the right. Different lowercase letters denote significant differences among diapause time by Tukey test at p < 0.05.
Functional module annotation of carbohydrate-active enzymes (CAZy) primarily focused on glycosyl transferases (GT) and glycoside hydrolases (GH) (Figure 5A). The mean abundance of larval GT was significantly higher on D14 and D28 compared to the controls (p < 0.05, F = 7.89, df = 2; p < 0.05, F = 7.07, df = 2; Figure 5B). Additionally, the level of larval GH was significantly lower on D28 compared to the controls (p < 0.05, F = 8.53, df = 2; Figure 5B).
Figure 5. Analysis of the function module annotation of carbohydrate-active enzymes (CAZys) (A). Effects of diapause time on the mean abundance of CAZys (B). Different lowercase letters denote significant differences among diapause time by Tukey test at p < 0.05.
3.5 Tetracycline is the most effective antibiotic against E. casseliflavus
Eight antibiotics were evaluated for their effectiveness against E. casseliflavus, the dominant bacterium in the midgut of diapause soybean hawkmoths. TC, CHL and AMP have been determined as the most efficient antibiotics against E. casseliflavus in diapause larvae (see Figure 6). The mediums treated with these three antibiotics demonstrated more distinct bacteriostatic zones in comparison to others (p < 0.05, F = 60.75, df = 7; Figure 6A). The bacteriostatic diameter resulting from treatments with TC, CHL and AMP significantly exceeded that of the other five antibiotics (p < 0.05, F = 60.75, df = 7; Figure 6B).
Figure 6. Screening of antibiotics targeted on E. casseliflavus. The susceptibility paper tests (A). The bacteriostatic diameter in the media treated with different antibiotics (B). AMP, ampicillin; VAN, vancomycin; CIP, ciprofloxacin; GM, gentamicin; EM, erythromycin; LVX, levofloxacin; CHL, chloromycetin; TC, tetracycline. Different lowercase letters denote significant differences among different antibiotics by Tukey test at p < 0.05.
4 Discussion
As reported by Jing et al. (2020), the gut microbiota played a vital role in allocating and mobilizing nutrients in insects during diapause. This study employed metagenomic analysis to identify the dominant gut bacteria in the diapause larvae of soybean hawkmoths from Jiangsu Province. The dominant genera of gut bacteria found in diapause larvae were determined to be Enterococcus and Enterobacter. Our previous study on 16S rRNA sequencing analysis showed that Enterococcus and Stenotrophomonas were the predominant taxa at the genus level (Qian et al., 2023b). Thus, Enterococcus was certainly evidenced to be the predominant genus in the diapause larvae of soybean hawkmoths, while Enterobacter belongs to Proteobacteria, it was generally involved in the immunity and pathogenicity of insects. Xie et al. (2023) also reported that Enterococcus was the predominant genus found in soybean hawkmoths from Henan Province. Similar findings have been observed in other Lepidopteran species. For example, Kong et al. (2023) found that Enterococcus spp. were the predominant gut bacteria in Galleria mellonella larvae, with Enterobacter spp. being dominant in the pupae. However, Lv and Liu (2009) reported that Staphylococcus and Bacillus were the dominant bacterial genera in diapause soybean hawkmoths from Shandong Province by combining a combination of morphological identification and 16S DNA sequencing. Differences in results may be stem from variations in geographic populations and genetic differentiation, as well as advancements in next-generation sequencing (NGS) techniques. In recent years, metagenomic analysis has been widely used to study the gut microbiome of numerous species. The metagenomic analysis conducted in this study indicated that Enterococcus casseliflavus and Enterococcus pernyi, both belonging to the Enterococcus genus, were the dominant gut bacterial species detected in diapause larvae of soybean hawkmoth. It represents the first identification of such gut bacteria in soybean hawkmoths. This study advances our understanding of the interaction between gut microbes and insect hosts by identifying the dominant bacterium from genus to species levels.
The diapause stage is critical for insect development and presents challenges for both hosts and microbial partners (Dittmer and Brucker, 2021). This study examined how diapause time impacts composition, diversity and functionality of gut microbiomes present in soybean hawkmoths. The relative abundance of Enterococcus, including E. casseliflavus and E. pernyi, initially decreased on diapause D14 and slightly increased on D28. For the species diversity analysis, D14 was significantly separated from both D0 and D28, but D0 and D28 appear to overlap considerably. It suggested that the larvae may potentially survive harsh diapause conditions with the assistance of their dominant gut bacteria. Enterococcus is a predomiant gut commensal bacterium that plays a crucial role in the physiology and health of the major economic agricultural insect, B. mori (Zhang et al., 2022). Enterococcus faecium has been demonstrated to regulate the developmental genes in honey bees (Du et al., 2021). Kong et al. (2023) provided evidence that Enterococcus spp. significantly contributes to insect metamorphosis. Furthermore, Chao and ACE showed an increase with diapause time, indicating a greater abundance of gut microbiota in late diapause larvae.
Gut microbiota reportedly contributes significantly to insects’ nutritional metabolism, digestive function and endocrine signals (Fraune and Bosch, 2010; Dittmer and Brucker, 2021). Gut microbes are capable of consistently supplying insect hosts with some nutrients (e.g., EAAs, carbohydrates) which can boost metabolic production (Fraune and Bosch, 2010). During diapause, the axenic larvae of parasitic wasp Nasonia vitripennis consistently maintained lower levels of glucose and glycerol than conventional larvae (Dittmer and Brucker, 2021). In this study, the function of larval gut microbiota was primarily focused on metabolism, particularly in carbohydrate metabolism and amino acid metabolism. Chen B. et al. (2018) used shotgun metagenomic sequencing techniques to reveal that the gut microbiome of B. mori equally emphasizes the metabolism of both carbohydrates and amino acids. Additionally, the metabolic pathways of amino acids in the gut microbiota of soybean hawkmoths were up-regulated with extended diapause time, as indicated by the annotated data. An opposite trend was observed between the abundance of the dominant gut bacterium Enterococcus and the amino acid metabolic pathways with diapause time. Similarly, the relative abundance of Enterococcus in B. mori had a significant correlation with the concentration of amino acids (Li et al., 2022). Liang et al. (2022) reported that E. casseliflavus is capable of producing L-tryptophan in B. mori. It suggested that Enterococcus may upregulate amino acid metabolic pathways to maintain the nutritional balance in diapause soybean hawkmoths.
Based on the previous non-target metabolomics analysis on diapause soybean hawkmoths, it showed that the most enriched pathways were mainly those involved in amino acid metabolism. The differential metabolite D-glutamine was found to be involved in arginine biosynthesis pathways (Qian et al., 2023b). D-glutamine has been reported to assist in nutrient delivery (particularly amino acids) for the metabolism during insects’ immune response to adverse conditions, such as low temperatures or food restrictions (Dolezal et al., 2019). Moreover, the content of arginine in diapause larvae of soybean hawkmoths (Qian et al., 2023b) also exhibited an opposite alteration to the abundance of dominant gut bacterium E. casseliflavus. Therefore, it is speculated that soybean hawkmoths probably up-regulate arginine metabolism by reducing the abundance of E. casseliflavus, thus ensuring their nutritional balance during diapause.
E. casseliflavus, the dominant gut bacterium in diapause larvae of soybean hawkmoths, exhibits antibiotic resistance, complicating the investigation of its interaction with insect hosts. Tetracycline (TC), chloromycetin (CHL) and ampicillin (AMP) were screened as the most effective antibiotics for targeting it in this study. These results align with those found in other Lepidopteran insects. For example, Li (2022) applied TC and AMP to treat E. casseliflavus in Cnaphalocrocis medinalis, and TC was used to treat Tetranychus truncatus (Yang et al., 2020). Next, the mixed liquor of TC, CHL and AMP will be used to determine the nutritional function of E. casseliflavus affecting on C. bilineata tsingtauica. It will be helpful to the farming industry of C. bilineata tsingtauica via the manual intervention of intestinal microecology, and provide scientific basis for the resource exploitation and utilization of edible insects.
5 Conclusion
In summary, Enterococcus casseliflavus was identified as the dominant bacterium in diapause larvae of C. bilineata tsingtauica, with the gene function focusing on carbohydrate metabolism and amino acid metabolism. In addition, we demonstrated that diapause time reduced the relative abundance of E. casseliflavus, while concurrently enhancing the relative abundance of amino acid metabolic pathways. Building on our previous non-target metabolomics research, it is possible that diapause larvae of C. bilineata tsingtauica may up-regulate their arginine metabolism by reducing the levels of E. casseliflavus to maintain their nutritional balance. This not only extends our comprehension of dominant gut bacteria from the genus to the species level, but also offers an initial investigation of gene functionality in soybean hawkmoths. Furthermore, TC was screened as the most effective antibiotic target for E. casseliflavus in soybean hawkmoths. Nevertheless, further verification is required for the effect of TC in vivo experiments and the appropriate treatment concentrations. Next, we will systematically measure the gut microbiota of soybean hawkmoths from newly hatched larvae to adults, and investigate the interaction mechanism between E. casseliflavus and insect hosts.
Data availability statement
The original contributions presented in the study are publicly available. This data can be found here: https://ncbi.nlm.nih.gov/bioproject/PRJNA1043846.
Author contributions
LQ: Conceptualization, Funding acquisition, Investigation, Writing – original draft. YW: Investigation, Methodology, Writing – original draft. PD: Investigation, Methodology, Writing – original draft. JZ: Investigation, Writing – original draft. YQ: Data curation, Investigation, Writing – original draft. ZL: Data curation, Investigation, Writing – original draft. HL: Project administration, Writing – review & editing. FC: Supervision, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by the National Natural Science Foundation of China (NSFC) (32100398), and the Natural Science Foundation of Jiangsu Province (BK20231393).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2024.1343265/full#supplementary-material
Footnotes
References
Anand, A. A. P., Vennison, S. J., Sankar, S. G., Prabhu, D. I. G., Vasan, P. T., Raghuraman, T., et al. (2010). Isolation and characterization of bacteria from the gut of Bombyx mori that degrade cellulose, xylan, pectin and starch and their impact on digestion. J. Insect Sci. 10, 1–20. doi: 10.1673/031.010.10701
Bemani, M., Izadi, H., Mahdian, K., Khani, A., and Samih, M. A. (2012). Study on the physiology of diapause, cold hardiness and supercooling point of overwintering pupae of the pistachio fruit hull borer, Arimania comaroffi. J. Insect Physiol. 58, 897–902. doi: 10.1016/j.jinsphys.2012.04.003
Bray, J. R., and Curtis, J. T. (1957). An ordination of the upland forest communities of southern Wisconsin. Ecol. Monogr. 27, 325–349. doi: 10.2307/1942268
Bu, D., Luo, H., Huo, P., Wang, Z., Zhang, S., He, Z., et al. (2021). KOBAS-i: intelligent prioritization and exploratory visualization of biological functions for gene enrichment analysis. Nucleic Acids Res. 49, W317–W325. doi: 10.1093/nar/gkab447
Cantalapiedra, C. P., Hernández-Plaza, A., Letunic, I., Bork, P., and Huerta-Cepas, J. (2021). eggNOG-mapper v2: functional annotation, orthology assignments, and domain prediction at the metagenomic scale. Mol. Biol. Evol. 38, 5825–5829. doi: 10.1093/molbev/msab293
Chen, G., Li, Q., Yang, X., Li, Y., Liu, W., Chen, F., et al. (2023). Comparison of the co-occurrence patterns of the gut microbial community between Bt-susceptible and Bt-resistant strains of the rice stem borer, Chilo suppressalis. J. Pest. Sci. 96, 299–315. doi: 10.1007/s10340-022-01512-5
Chen, B., Yu, T., Xie, S., Du, K., Liang, X., Lan, Y., et al. (2018). Comparative shotgun metagenomic data of the silkworm Bombyx mori gut microbiome. Sci. Data 5:180285. doi: 10.1038/sdata.2018.285
Chen, S., Zhou, Y., Chen, Y., and Gu, J. (2018). Fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 34, i884–i890. doi: 10.1093/bioinformatics/bty560
Dittmer, J., and Brucker, R. M. (2021). When your host shuts down: larval diapause impacts host-microbiome interactions in Nasonia vitripennis. Microbiome 9:85. doi: 10.1186/s40168-021-01037-6
Dolezal, T., Krejcova, G., Bajgar, A., Nedbalova, P., and Strasser, P. (2019). Molecular regulations of metabolism during immune response in insects. Insect Biochem. Mol. Biol. 109, 31–42. doi: 10.1016/j.ibmb.2019.04.005
Douglas, A. E. (2015). Multiorganismal insects: diversity and function of resident microorganisms. Annu. Rev. Entomol. 60, 17–34. doi: 10.1146/annurev-ento-010814-020822
Du, Y., Luo, S., and Zhou, X. (2021). Enterococcus faecium regulates honey bee developmental genes. Int. J. Mol. Sci. 22:12105. doi: 10.3390/ijms222212105
Engel, P., and Moran, N. A. (2013). The gut microbiota of insects - diversity in structure and function. FEMS Microbiol. Rev. 37, 699–735. doi: 10.1111/1574-6976.12025
Fraune, S., and Bosch, T. C. (2010). Why bacteria matter in animal development and evolution? BioEssays 32, 571–580. doi: 10.1002/bies.200900192
Gandotra, S., Bhuyan, P. M., Gogoi, D. K., Kumar, A., and Subramanian, S. (2018). Screening of nutritionally important gut bacteria from the lepidopteran insects through qualitative enzyme assays. Proc. Natl Acad. Sci. India Section B Biol. Sci. 88, 329–337. doi: 10.1007/s40011-016-0762-7
Gao, H. H., Zhao, S., Wang, R. J., Qin, D. Y., Chen, P., Zhang, A. S., et al. (2023). Gut bacterium promotes host fitness in special ecological niche by affecting sugar metabolism in Drosophila suzukii. Insect Sci. 30, 1713–1733. doi: 10.1111/1744-7917.13189
Genta, F. A., Dillon, R. J., Terra, W. R., and Ferreira, C. (2006). Potential role for gut microbiota in cell wall digestion and glucoside detoxification in Tenebrio molitor larvae. J. Insect Physiol. 52, 593–601. doi: 10.1016/j.jinsphys.2006.02.007
Guo, M. M., Li, X. F., Deng, P., Li, D. W., Li, J. L., Fan, J. W., et al. (2021). Diapause termination and post-diapause of overwintering Clanis bilineata tsingtauica larvae. Chin. J. Appl. Entomol. 58, 966–972.
Hahn, D. A., and Denlinger, D. L. (2011). Energetics of insect diapause. Annu. Rev. Entomol. 56, 103–121. doi: 10.1146/annurev-ento-112408-085436
Jing, T. Z., Qi, F. H., and Wang, Z. Y. (2020). Most dominant roles of insect gut bacteria: digestion, detoxification, or essential nutrient provision? Microbiome 8:38. doi: 10.1186/s40168-020-00823-y
Kong, H. G., Son, J. S., Chung, J. H., Lee, S., Kim, J. S., and Ryu, C. M. (2023). Population dynamics of intestinal Enterococcus modulate galleria mellonella metamorphosis. Microbiol. Spectr. 11:e0278022. doi: 10.1128/spectrum.02780-22
Kuang, J., Wang, J., Li, Y., Li, M., Zhao, M., Ge, K., et al. (2023). Hyodeoxycholic acid alleviates non-alcoholic fatty liver disease through modulating the gut-liver axis. Cell Metab. 35, 1752–1766.e8. doi: 10.1016/j.cmet.2023.07.011
Lee, J., Kim, C. H., Jang, H. A., Kim, J. K., Kotaki, T., Shinoda, T., et al. (2019). Burkholderia gut symbiont modulates titer of specific juvenile hormone in the bean bug Riptortus pedestris. Dev. Comp. Immunol. 99:103399. doi: 10.1016/j.dci.2019.103399
Li, C. (2022). Effect of food nutrition an gut microbiota on growth performance of the rice leaffolder Cnaphalocrocis medinalis (Güenée), vol. 2. Yangzhou, China: Yangzhou University, 60.
Li, D., Liu, C. M., Luo, R., Sadakane, K., and Lam, T. W. (2015). MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics 31, 1674–1676. doi: 10.1093/bioinformatics/btv033
Li, G., Zheng, X., Zhu, Y., Long, Y., and Xia, X. (2022). In-depth insights into the disruption of the microbiota-gut-blood barrier of model organism (Bombyx mori) by fluoride. Sci. Total Environ. 838:156220. doi: 10.1016/j.scitotenv.2022.156220
Liang, X., He, J., Zhang, N., Muhammad, A., Lu, X., and Shao, Y. (2022). Probiotic potentials of the silkworm gut symbiont Enterococcus casseliflavus ECB140, a promising L-tryptophan producer living inside the host. J. Appl. Microbiol. 133, 1620–1635. doi: 10.1111/jam.15675
Lv, F., and Liu, Y. (2009). The identification of dominant intestinal bacteriae of Clanis bilineata tsingtauica Mell larva. Chin. J. Microecol. 21, 435–439.
Patro, R., Duggal, G., and Kingsford, C. (2015). Salmon: accurate, versatile and ultrafast quantification from RNA-seq data using lightweight-alignment. BioRxiv. doi: 10.1101/021592
Poelchau, M. F., Reynolds, J. A., Elsik, C. G., Denlinger, D. L., and Armbruster, P. A. (2013). Deep sequencing reveals complex mechanisms of diapause preparation in the invasive mosquito, Aedes albopictus. Proc. R. Soc. B Biol. Sci. 280:20130143. doi: 10.1098/rspb.2013.0143
Qian, L., Chen, B. J., Gui, F. R., Qin, Y., Deng, P., and Liao, H. J. (2023a). Nutritional and feeding adaptability of Clanis bilineata tsingtauica larvae to different cultivars of soybean, (Glycine max). Food Secur. 12:1721. doi: 10.3390/foods12081721
Qian, L., Deng, P., Chen, F., Cao, Y., Sun, H., and Liao, H. (2022). The exploration and utilization of functional substances in edible insects: a review. Food Prod. Process. Nutr. 4:11. doi: 10.1186/s43014-022-00090-4
Qian, L., Chen, B. J., Deng, P., Gui, F. R., Cao, Y., Qin, Y., et al. (2023b). TM7 (Saccharibacteria) regulates the synthesis of linolelaidic acid and tricosanoic acid, and alters the key metabolites in diapause Clanis bilineata tsingtauica. Front. Physiol. 14:1093713. doi: 10.3389/fphys.2023.1093713
Ramette, A. (2007). Multivariate analyses in microbial ecology. FEMS Microbiol. Ecol. 62, 142–160. doi: 10.1111/j.1574-6941.2007.00375.x
Santana, C. C., Barbosa, L. A., Junior, I. D. B., Nascimento, T. G., Dornelas, C. A., and Grillo, L. A. M. (2017). Lipase activity in the larval midgut of Rhynchophorus palmarum: biochemical characterization and the effects of reducing agents. Insects 8:100. doi: 10.3390/insects8030100
Segata, N., Izard, J., Waldron, L., Gevers, D., Miropolsky, L., Garrett, W. S., et al. (2011). Metagenomic biomarker discovery and explanation. Genome Biol. 12:R60. doi: 10.1186/gb-2011-12-6-r60
Shin, S. C., Kim, S. H., You, H., Kim, B., Kim, A. C., Lee, K. A., et al. (2011). Drosophila microbiome modulates host developmental and metabolic homeostasis via insulin signaling. Science 334, 670–674. doi: 10.1126/science.1212782
Steinegger, M., and Söding, J. (2017). MMseqs2 enables sensitive protein sequence searching for the analysis of massive data sets. Nat. Biotechnol. 35, 1026–1028. doi: 10.1038/nbt.3988
Storelli, G., Defaye, A., Erkosar, B., Hols, P., Royet, J., and Leulier, F. (2011). Lactobacillus plantarum promotes Drosophila systemic growth by modulating hormonal signals through TOR-dependent nutrient sensing. Cell Metab. 14, 403–414. doi: 10.1016/j.cmet.2011.07.012
Tian, H., and Zhang, Y. M. (2012). The nutritional components analysis and evaluation of Clanis bilineata tsingtauica Mell. Acta Nutrimenta Sin. 34, 289–291.
van Opstal, E. J., and Bordenstein, S. R. (2019). Phylosymbiosis impacts adaptive traits in nasonia wasps. MBio 10:e0088719. doi: 10.1128/mBio.00887-19
Visotto, L. E., Oliveira, M. G. A., Guedes, R. N. C., Ribon, A. O. B., and Good-God, P. I. V. (2009). Contribution of gut bacteria to digestion and development of the velvetbean caterpillar, Anticarsia gemmatalis. J. Insect Physiol. 55, 185–191. doi: 10.1016/j.jinsphys.2008.10.017
Xie, J., Li, X., Shi, Q., Wu, J., Zhu, J., Li, R., et al. (2023). Effects of Pueraria lobata leaf on growth performance and intestinal microflora diversity of Clanis bilineata tingtauica larvae. J. Northern Agric. 51, 107–117.
Yang, K., Chen, H., Bing, X. L., Xia, X., and Hong, X. Y. (2020). Antibiotic-induced changes in tetranychus truncatus bacterial community alter its fecundity, longevity and sex-ratio. Syst. Appl. Acarol. 25, 1668–1682. doi: 10.11158/saa.25.9.12
Zhang, X., Feng, H., He, J., Muhammad, A., Zhang, F., and Lu, X. (2022). Features and colonization dtrategies of Enterococcus faecalis in the gut of Bombyx mori. Front. Microbiol. 13:921330. doi: 10.3389/fmicb.2022.921330
Zhang, Y., and Ju, F. (2023). Uninheritable but widespread bacterial symbiont mediates insecticide detoxification of an agricultural invasive pest Spodoptera frugiperda. bioRxiv :559648. doi: 10.1101/2023.09.26.559648
Zheng, H., Powell, J. E., Steele, M. I., Dietrich, C., and Moran, N. A. (2017). Honeybee gut microbiota promotes host weight gain via bacterial metabolism and hormonal signaling. Proc. Natl. Acad. Sci. USA 114, 4775–4780. doi: 10.1073/pnas.1701819114
Zhu, W., Alexandre, L., and Mark, B. (2010). Ab initio gene identification in metagenomic sequences. Nucleic Acids Res. 38:e132. doi: 10.1093/nar/gkq275
Keywords: edible insect, Clanis bilineata tsingtauica, metagenomics, dominant gut bacterium, nutrient metabolism
Citation: Qian L, Wang Y, Deng P, Zhang J, Qin Y, Li Z, Liao H and Chen F (2024) Enterococcus casseliflavus regulates amino acid metabolism in edible insect Clanis bilineata tsingtauica: a functional metagenomics study. Front. Microbiol. 15:1343265. doi: 10.3389/fmicb.2024.1343265
Edited by:
Paul Dyson, Swansea University Medical School, United KingdomReviewed by:
Ningning Li, Sun Yat-Sen University, ChinaAbhay K. Pandey, Tea Research Association, India
Copyright © 2024 Qian, Wang, Deng, Zhang, Qin, Li, Liao and Chen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Huaijian Liao, bGhqQGphYXMuYWMuY24=; Fajun Chen, ZmFqdW5jaGVuMThAc2luYS5jb20=