Front. Microbiol. Frontiers in Microbiology Front. Microbiol. 1664-302X Frontiers Media S.A. 10.3389/fmicb.2023.1251627 Microbiology Original Research Variation of Helicoverpa armigera symbionts across developmental stages and geographic locations Zhao Chenchen 1 2 Wang Li 1 3 Zhang Kaixin 1 3 4 Zhu Xiangzhen 1 3 4 Li Dongyang 1 3 4 Ji Jichao 1 3 4 * Luo Junyu 1 4 * Cui Jinjie 1 4 * 1National Key Laboratory of Cotton Bio-breeding and Integrated Utilization, Institute of Cotton Research, Chinese Academy of Agricultural Sciences, Anyang, Henan, China 2Henan International Laboratory for Green Pest Control, College of Plant Protection, Henan Agricultural University, Zhengzhou, Henan, China 3Zhengzhou Research Base, State Key Laboratory of Cotton Biology, Zhengzhou University, Zhengzhou, Henan, China 4Western Agricultural Research Center, Chinese Academy of Agricultural Sciences, Changji, China

Edited by: Kiran Pawar, Shivaji University, India

Reviewed by: Guan-Hong Wang, Chinese Academy of Sciences (CAS), China; J. P. Michaud, Kansas State University, United States

 *Correspondence: Jichao Ji hnnydxjc@163.com Junyu Luo luojunyu1818@126.com Jinjie Cui aycuijinjie@163.com 07 09 2023 2023 14 1251627 02 07 2023 18 08 2023 Copyright © 2023 Zhao, Wang, Zhang, Zhu, Li, Ji, Luo and Cui. 2023 Zhao, Wang, Zhang, Zhu, Li, Ji, Luo and Cui

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.

Cotton bollworm (Helicoverpa armigera) poses a global problem, causing substantial economic and ecological losses. Endosymbionts in insects play crucial roles in multiple insect biological processes. However, the interactions between H. armigera and its symbionts have not been well characterized to date. We investigated the symbionts of H. armigera in the whole life cycle from different geographical locations. In the whole life cycle of H. armigera, Proteobacteria, Firmicutes, Bacteroidetes, and Actinobacteria were the dominant bacteria at the phylum level, while Enterococcus, Enterobacter, Glutamicibacter, and Bacillus were the four dominant bacteria at the genus level. Furthermore, high similarity in symbiotic bacterial community was observed in different stages of H. armigera, which were dominated by Enterococcus and Enterobacter. In fields, the dominant bacteria were Proteobacteria and Bacteroidetes, whereas, in the laboratory, the dominant bacteria were Proteobacteria. At the genus level, the dominant bacteria in cotton bollworm eggs of wild populations were Enterobacter, Morganella, Lactococcus, Asaia, Apibacter, and Enterococcus, and the subdominant bacteria were Bartonella, Pseudomonas, and Orbus. Moreover, the symbionts varied with geographical locations, and the closer the geographical distance, the more similar the microbial composition. Taken together, our study identifies and compares the symbiont variation along with geographical gradients and host development dynamic and reveals the high flexibility of microbiome communities in H. armigera, which probably benefits for the successful survival in a complicated changing environment.

 cotton bollworm development stage symbionts geographical population dynamic change section-at-acceptance Microbial Symbioses Introduction

Host physiological and biochemical functions (development, physiology, ecological interactions, and evolutionary diversity) are closely associated with microbes (Zhu et al., 2022). Host-associated microbiota studies have indicated that microbes can influence host fitness (Jia et al., 2021). In addition, symbionts play key roles in insect development (Hammer et al., 2017). For instance, the nitrogen-recycling activity of gut bacteria in turtle ants (Cephalotes) is involved in several biosynthetic pathways contributing to host cuticle formation (Duplais et al., 2021), and burying beetle (Nicrophorus vespilloides) preserves carcasses by regulating its microbial growth, which is important for optimal larval development (Shukla et al., 2018). Interestingly, gut microbes of caterpillars can extend to not only just insects but also plant phenotypes and alter outcomes of plant–herbivore interaction (Mason et al., 2019). Deciphering the composition and function of bacterial symbionts and their effects on the hosts poses a significant challenge (Zhao et al., 2019).

The cotton bollworm Helicoverpa armigera (Lepidoptera: Noctuidae) is a cosmopolitan pest with a wide host range, and it usually causes huge economic and ecological losses each year. Helicoverpa armigera has historically been controlled by transgenic cotton-containing toxic protein from Bacillus thuringiensis (Bt) (Downes et al., 2017; Wilson et al., 2018). The Bt cotton planting area has increased considerably in many countries, including China, since 1996 (Cui et al., 2011; Lu et al., 2012). However, in the recent 20 years, Bt cotton planting area contraction in North China Plain has caused the regional resurgence of cotton bollworm, huge crop yield loss, and excess pesticide use (Lu et al., 2022). Considerable studies have demonstrated that symbiont of H. armigera can protect the host from baculovirus and Bt biopesticides (Xu et al., 2014; Yuan et al., 2021). However, as a highly polyphagous worldwide lepidopteran pest and an important model system used in various biological research studies, the symbiont profiles and dynamics of H. armigera at developmental stages and different geographical locations remain largely unknown.

Helicoverpa armigera undergoes complete metamorphosis, and its larvae and adults have distinct morphological characteristics, habitat, feeding behavior, and function. To explore the response of symbionts to host ages and geographical populations, we systematically investigated the bacterial community by 16S ribosomal RNA (16S rRNA) gene sequencing. Our results provide an insight into the interactions between symbionts and H. armigera in different developmental stages or from multiple cotton planting regions and offer valuable references for further exploring the potential roles of microbiota in host cotton bollworm and developing new management strategies by disturbing bacterial communities.

 Materials and methods Insect rearing

Larvae of H. armigera were purchased from Henan Jiyuan Baiyun Industry Co., Ltd., Henan, China. Then, H. armigera were reared for over three generations without insecticide exposure in the laboratory. Adults were maintained in an environmental chamber and supplied with a 10% honey solution, whereas larvae were maintained on an artificial diet as previously described (Qin, 2015). Rearing conditions were as follows: at 26 ± 1°C, relative humidity of 70 ± 5%, and a 14-h:10-h light: dark cycle. Helicoverpa armigera at eight different development stages (eggs, 1st−5th instar larvae, pupae, and adults; Figure 1A) were randomly sampled from the laboratory population. The eight developmental stage experiments of H. armigera were performed with three replicates per stage and 20 individuals or 100 eggs per replicate. The samples were immediately flash-frozen in liquid nitrogen and stored at −80°C before DNA extraction.

Life cycle (A) and a map illustarating regions for sample collection of Helicoverpa armigera (B).

 Collection of insect field populations

Helicoverpa armigera were sampled from six geographical locations in the Yellow River cotton planting region of China (mainly in Henan province, Hebei province, and Shandong province; Figure 1B). Adults were captured with 1,000 W light traps during the second generation.

Field-collected female adults from six locations were placed separately in 250 ml plastic cups covered with cotton gauze on which the female adult laid eggs. A piece of absorbent cotton was placed in each cup containing 10% honey solution for feeding the adults. Female adults were kept at 27–30°C, 70%−80% relative humidity (RH), and light:dark (L:D) = 14:10. We placed the single mated female in a small plastic cup (diameter1 = 52 mm, diameter2 = 75 mm, and height = 85 mm, covered by sterile plastic paper) for egg collection.

The eggs were collected daily and used for the subsequent experiment. The cotton bollworms captured from each geographical location were divided into five groups, and eggs were collected from each group five times with 300 eggs per time.

 DNA extraction, PCR amplification, and sequencing

Total genomic DNA was extracted from cotton bollworm samples at different development stages or from different geographic populations, respectively, using the TIANamp Genomic DNA Kit (Tiangen Biotech Inc., China), according to the manufacturer's instructions. The surface of the larvae was cleaned with 75% ethanol and rinsed thrice with sterile water before DNA extraction. The lysozyme (50 mg/ml) was added to samples and incubated for 30 min at 37°C, to break up gram-positive bacterial cells. The quantity and quality of the DNA were detected with a NanoDrop 2000C spectrophotometer (Thermo Scientific, USA) and agarose gel electrophoresis, respectively. The V3–V4 variable region of the 16S rRNA gene was amplified using 338F/806R primers (338F: 5′- ACTCCTACGGGAGGCAGCAG-3′, 806R: 5′- GGACTACHVGGGTWTCTAAT-3′) (Xu et al., 2016). The quantification, qualification, and purification of PCR products (Yeasen, China), library preparation, and sequencing were conducted as previously reported (Zhang et al., 2019). Purified amplicons were pooled in equimolar concentration and paired-end sequenced (2 × 300) on an Illumina MiSeq platform (Illumina, San Diego, USA), according to the standard protocols by Shanghai Majorbio Bio-pharm Technology Co., Ltd. The raw reads were submitted to the NCBI Sequence Read Archive (SRA) database with an accession number PRJNA591375.

 Bioinformatics analysis

Data were analyzed as described in previous studies (Zhao et al., 2019, 2020). Based on their unique barcode, paired-end reads were assigned to the library and truncated according to the barcode and primer sequence, to obtain a standard sequence. The obtained sequence reads were processed using QIIME (Caporaso et al., 2010). Partial 16S rRNA bacterial sequences were filtered using Mothur (Schloss et al., 2009), with the inclusion criteria of mean quality score ≥20 and length ≥250 bp. To classify the entire microbial community in H. armigera, reads were clustered into operational taxonomic units (OTUs) with a 97% similarity cutoff. Sequence analysis was performed by UPARSE (Edgar, 2013). The rarefied OTU tables were generated to prevent possible sample heterogeneity due to their different sequence numbers. Chao1 and Ace (abundance-based coverage estimator) indices and Shannon and Simpson indices were used to evaluate microbial community richness and diversity, respectively. A Venn diagram was drawn using Venn diagram plotter software (https://omics.pnl.gov/software/venn-diagram-plotter), to visualize unique and shared OTUs across all evaluated samples. Bray–Curtis distance matrices were visualized through principal coordinate analysis (PCoA). UniFrac analysis compared microbial diversity with different samples (Hamady et al., 2010). To investigate community composition differences, we conducted a non-parametric multivariate analysis of variance using the R vegan package (Oksanen et al., 2009). Differences in dominant bacteria were analyzed via one-way analysis of variance (ANOVA), followed by Tukey's honestly significant difference (HSD) test using SAS software. Non-parametric Kruskal–Wallis test was performed when data were not normally distributed (P < 0.05). These analyses were carried out by Statistics Analysis System software (SAS 9.4, SAS Institute Inc.).

We used Moran's I to analyze the spatial autocorrelation of the dominant bacteria including Enterobacter, Morganella, Lactococcus, Asaia, and Apibacter in field-collected cotton bollworms. A structural equation model (SEM) (Tenenhaus et al., 2005) with a Satorra–Bentler correction was used to evaluate the effects of geographic or climatic factors on the predominant symbionts in cotton bollworm (Jiang et al., 2023). The temperature and humidity data were provided by the China Meteorological Data Service Center (CMDC) and the global climate database (WorldClim; https://www.worldclim.org/).

 Quantification of bacteria communities

The copies of 16S rRNA genes of 32 dominant bacteria were determined by a previously reported method (Zhang et al., 2019). The standard vector was serially diluted to generate standard curves. The target sequence in the standard vectors was identified by sequencing. The target sequence was amplified through PCR. PCR was performed on a StepOnePlus™ Real-Time PCR System (Applied Biosystems, Foster City, CA, USA). The 20 μl PCR system contained 10 μl 2 × TransStart Green qPCR SuperMix (TransGen Biotech Co., Ltd, China), 0.4 μl each of 10 mM forward and reverse primers, 1.0 μl template DNA (12.5 ng/μl and 2 μl DNA extraction control DNA considered as negative controls), 0.4 μl 50 × ROX, and 7.8 μl H2O. The PCR program was as follows: 95°C for 3 min, followed by 40 cycles of a two-step PCR (95°C for 5 s and 60°C for 30 s). Primers used for qPCR analysis are presented in Supplementary Table 1.

 Results Overview of <italic>H. armigera</italic> microbiota

Microbial communities in cotton bollworm at different developmental stages or from multiple geographical areas were determined through 16S rRNA sequencing. A total of 1,671 and 1,443 OTUs were identified from cotton bollworm samples of developmental stages or geographical areas, respectively, with an average length of 428 bp (Supplementary Tables 1, 2). These OTUs were clustered into 25 phyla, 548 genera, and 898 species. The rarefaction curves for all samples approached saturation after the number of sequences reached 18,000 (developmental stages) and 35,000 (geographical areas; Supplementary Figure 1), indicating most microbial species were captured in our study. The average sequencing coverage rate of each sample was more than 99.6% (Supplementary Tables 1, 3), which further confirmed that the experimental data accurately reflected the composition of most of the bacterial community.

 Effects of host development stages on microbiome community

Taxonomic analysis showed that Proteobacteria was the most prevalent phylum. Proteobacteria and Firmicutes were the most abundant phyla across all samples (Figure 2A). The relative abundance of Proteobacteria was high in eggs (68.56%), 1st instar larvae (62.46%), 5th instar larvae (85.44%), pupae (97.63%), and female adults (47.89%) and very low in 2nd instar larvae (10.48%), 3rd instar larvae (8.56%), and 4th instar larvae (6.47%). In contrast, the predominant phyla in 2nd, 3rd, and 4th instar larvae were Firmicutes with the relative abundance of 87.73%, 91.04%, and 93.41%, respectively (Figure 2A). At the genus level, there were significant differences in microbial composition and relative abundance at different developmental stages of cotton bollworm (Figure 2B). In contrast to the irregularity of bacteria composition in other developmental stages of H. armigera, the 2nd, 3rd, and 4th instar (intermediate instar) larvae exhibited similar microbial community composition, and their dominant bacteria were Enterococcus (>87%), which was also dominant in eggs (29.91%). In egg, pupa, female adult, and 1st instar larvae, the dominant bacterium was Enterobacter (accounting for 67.09%, 91.26%, 30.22%, and 40.50%). Interestingly, the dominant bacterial genera of the male adult were Glutamicibacter (52.26%) and Bacillus (11.12%), respectively. The dominant bacteria of the 5th instar larvae were Acinetobacter (37.85%), Ralstonia (22.62%), and Sphingomonas (15.66%; Figures 2B, C).

Bacterial community dynamics among different developmental stages in Helicoverpa armigera. Relative abundance of bacteria communities at the phylum level (A) and genus level (B) in different groups (Non-parametric Kruskal–Wallis test *0.01 < P ≤ 0.05, **0.001 < P ≤ 0.01, and ***P ≤ 0.001). (C) Heat map of major taxa over the life cycle at the genus level. Cluster analysis using the Bray–Curtis distance and the complete-linkage method. Each column consisted of top 20 bacteria genus from the eight treatments. Columns are clustered according to the similarity of bacterial abundance profiles. Color gradient represents the proportion of species. Plotting scale, from red to blue, indicates the decrease in richness of bacterial communities. Alpha diversity comparisons of bacterial communities across life stages of H. armigera (D, E). Different lowercase labels above each group indicate significant differences (one-way ANOVA, Tukey's HSD test, P < 0.05) of group mean value. PCoA plots (F) and NMDS plot (G) based on weighted-UniFrac of bacterial communities across life stages of H. armigera.

To further investigate the microbial community present in H. armigera larvae, we compared multiple indices of community richness (OTUs, Chao1, ACE, and PD whole tree), diversity (Shannon and Simpson), and sequencing depth (Good's coverage which represents the percentage of the total microbial species detected in a sample; Figures 2D, E, Supplementary Table 1). The average of Good's coverage was estimated to be 99.39%, suggesting that the overwhelming majority of microbial species present in H. armigera larvae were included in this study (Supplementary Table 1). The OTU number and Shannon index data showed that the microbiome exhibited large dynamic changes across the eight stages, indicating that H. armigera tended to harbor a relatively dynamic microbial diversity. One-way ANOVA showed that the ACE, Chao1, Simpson, and PD tree indices significantly differed across the eight stages (Figures 2C–G). The 1st instar larva and adult displayed higher species richness and diversity than other stages. The microbial community richness in eight stages ranked as follows: 1st instar larvae > male adult > female adult > 5th instar larvae > 2nd instar larvae > egg > pupa> 3rd instar larvae > 4th instar larvae (Figure 2D). The microbial community diversity in eight stages ranked as follows: 1st instar larvae > male adult > female adult > 2nd instar larvae >egg > 3rd instar larvae > pupa > 5th instar larvae > 4th instar larvae (Figure 2E). Taken together, these results indicated that the bacterial community structures in H. armigera exhibited dynamic variation across the eight stages.

We performed unconstrained principal coordinate analyses (PCoAs) of UniFrac distances (weighted and unweighted) and Bray–Curtis metrics to investigate patterns of separation among microbial communities in H. armigera at different developmental stages (Figures 2C, F, G, Supplementary Table 1). The results showed that the microbiome from cotton bollworm at each of the eight developmental stages clustered separately and exhibited a significant difference in bacterial composition (Adonis, P < 0.05; Figures 2F, G). The microbial community composition of the intermediate instar larvae (instars 2–4) was similar, and that of the pupa, egg, and female adult was similar, while that of the 1st instar larvae, 5th instar larvae, and male adult was significantly different from that at other developmental stages (Figure 2F). Complementary analyses were conducted by non-metric multidimensional scaling (NMDS; stress < 0.1, R = 0.6008, P = 0.001, by ADONIS), which was consistent with those results through PCoA, further indicating the reliability of sample clustering results based on microbial community composition (Figure 2G). The microorganism communities at the 1st instar larvae, the 5th instar larvae, and a male adult significantly differed from those at other instars. The intermediate (2nd to 4th instar) larval microbial communities were similar (Figure 2G). These results suggested significant differences in the microorganism communities of cotton bollworm across different developmental stages.

 Effects of different geographical populations on bacterial community diversity in <italic>H. armigera</italic>

In this study, the richness and diversity of microbial communities in cotton bollworm eggs laid by six field populations and one laboratory population were compared (Figure 3, Supplementary Table 2), with the laboratory population as the control. The results showed that the laboratory population (Lab) exhibited the lowest Sob (observed species) index, followed by the Qiuxian population (QX), which was lower than that of the other three geographical populations, whereas the Binzhou population (BZ) displayed the highest Sob index (P < 0.05; Figures 3A, B). The community diversity and species richness of the laboratory population (Lab) were the lowest, followed by the Qiuxian population. Microbial richness increased progressively across the Cangzhou, Anyang, Yanggu, and Binzhou populations but was not significant. Except for the Anyang population, the species richness of geographic populations was inversely proportional to longitude with an order of Qiuxian population < Yanggu population < Cangzhou population < Binzhou population (Figure 3). Simpson, Ace, Chao1, and PD tree indices of the microbial community from these four geographical populations exhibited similar pattern (Supplementary Table 2), which further confirmed the association between the microbial species richness and longitude of host geographic population. The microbial species richness and diversity indices of microbiotas in eggs laid by the laboratory population were also significantly lower than those of field populations (Figures 3D, E, Supplementary Table 2).

Bacterial community dynamics among different developmental stages in Helicoverpa armigera. Relative abundance of bacteria communities at the phylum level (A, B) in different groups (Non-parametric Kruskal–Wallis test *0.01 < P ≤ 0.05, **0.001 < P ≤ 0.01, and ***P ≤ 0.001). Heat map of major taxa over the life cycle at the genus level (C). Cluster analysis using the Bray–Curtis distance and the complete-linkage method. Each column represents top 35 bacteria genus from a certain geographical populations of cotton bollworm and is clustered according to the similarity of bacterial abundance profiles. Boxplot of species richness (number of OTUs) (D) and community diversity measured by the Shannon index (E). Different lowercase labels above each group indicate significant differences (one-way ANOVA, LSD post-hoc test, P < 0.05) of group mean value. PCoA plots (F) and NMDS plot (G) based on weighted-UniFrac of bacterial communities across life stages of H. armigera.

To visualize differences in the bacterial community, the histograms of the top nine bacterial genera were constructed among the different sample groups using the QIIME toolkit (Figures 3A, B). According to the OTU representative sequences, Firmicutes, Proteobacteria, Actinobacteria, Bacteroidetes, Cyanobacteria, and Fusobacteria were identified. Firmicutes, Proteobacteria, and Actinobacteria were identified as the dominant community (accounting for >87%) in H. armigera at the phylum level (Figures 3A, B).

The heatmap showed the relative abundance of the top 50 genera from six geographical populations. At the genus level, the primary dominant bacteria in cotton bollworm eggs were Enterobacter, Morganella, Lactococcus, Asaia, Apibacter, and Enterococcus. Their total relative abundance reached 69.13% (Anyang), 54.80% (Binzhou), 68.13% (Cangzhou), 84.66% (Qiuxian), 95.13% (Lab), and 79.97% (Yanggu), respectively. At the genus level, the secondary dominant bacteria were Bartonella, Pseudomonas, and Orbus. The relative abundance of Pseudomonas in the Yanggu population was the highest (9.6%), and the relative abundance of Bartonella in the Cangzhou and Binzhou populations (6.4 and 7.99%) was higher than that in other populations (< 3%). The relative abundance of Acinetobacter in the Anyang population (5.46%) was higher than that in other geographical populations (< 0.5%). Except for Enterobacter, Enterococcus, and Asaia, the relative abundance of other bacteria in the laboratory population was extremely low, compared with that of the field population. Apibacter, Bartonella, Orbaceae, and Orbus were not even detected in the laboratory population (Figure 3C). The cluster analysis showed that the dominant microbial composition in the field populations was similar, the closer the distance was, the higher the similarity was (Anyang–Qiuxian and Cangzhou–Binzhou), and the dominant genus composition in the laboratory population was different. Most of the dominant bacteria in the field populations belong to Proteobacteria and Bacteroidetes, while most of the bacteria in the laboratory population belong to Proteobacteria. In summary, the composition of the microbiome varied greatly across geographic populations.

OTU information of each sample was used to establish the Bray–Curtis distance matrices and draw the PCoA plot, to analyze the differences in microbial composition among six geographic populations (Figure 3F). The results showed that two principal components could explain 35.88 and 11.40% of the total variation. ANOSIM showed that the intra-group differences were smaller than the inter-group differences. PCoA could significantly distinguish laboratory population from field population. The samples of the field populations were separated, indicating that the microbial community structure of the field populations differed greatly within the population. NMDS analysis results (stress < 0.152, R = 0.2875, P = 0.001) indicated that the grouping was reasonable, and the reliability of this test was high. These results of NMDS analysis were in line with those of PCoA. The laboratory population was relatively clustered and could be distinguished significantly from other populations. The microbial communities of different field populations were interwoven, indicating that the microbial community structure among wild populations was similar (Figure 3G).

The linear discriminant analysis effect size (LEfSe) was performed to analyze the microbiota of eggs laid by cotton bollworm from different geographical populations (with an LDA threshold of 4; Figure 4A). The result showed that Bacteroidetes, Flavobacteriia, Flavobacteriales, Flavobacteriaceae, Apibacter, and Chryseobacterium were dominant bacteria in the Anyang population. Betaproteobacteria, Rhizobiales, Neisseriaceae, and Commensalibacter mainly existed in the Binzhou population. Orbales, Orbaceae, Bartonellaceae, Bartonella, and Orbus mainly existed in the Cangzhou population. Bacteroidales, Streptococcaceae, Porphyromonadaceae, and Dysgonomonas mainly existed in the Qiuxian population. Proteobacteria, Gammaproteobacteria, Enterobacteriales, Enterobacteriaceae, and Enterobacter mainly existed in the laboratory population. Pseudomonadales, Pseudomonadaceae, Morganella, and Pseudomonas mainly existed in the Yanggu population.

Bacteria of Helicoverpa armigera populations with an LDA score of > 4 (A) and inferred functions of bacterial communities associated with H. armigera in different geographical populations (B) (*0.05 < P ≤ 0.01, **0.05 < P ≤ 0.01, and ***P ≤ 0.001).

Based on the 16S rRNA sequences and the KEGG results, we carried out a phylogenetic investigation of communities by reconstructing unobserved states (PICRUSt), to predict the gene functions of different geographic populations. The results showed that the genes exhibited differences in metabolism, environmental information processing, and genetic information processing (Figure 4B). The relative abundance of symbiotic bacteria involved in carbohydrate metabolism and secondary substance metabolism was significantly higher in the laboratory population than that in the field populations, but the relative abundance of symbiotic bacteria involved in sugar anabolic metabolism, nucleotide metabolism, energy metabolism, and cofactor and vitamin metabolisms was significantly lower in the laboratory population than that in the field populations (P < 0.05). Additionally, the relative abundance of symbiotic bacteria involved in lipid and amino acid metabolism was lower in the laboratory population than in the wild populations. The microbial communities in eggs oviposited by wild cotton bollworm populations were more involved in metabolic pathways. The relative abundance of symbiotic bacteria involved in genetic information processing (including translation, replication, repair, folding, classification, degradation, and transcription functions) was significantly lower in the laboratory population than in the wild population (P < 0.05). However, the relative abundance of symbiotic bacteria involved in membrane transport and signal transduction (belonging to environmental information processing) was higher in the laboratory population than in the wild population (P < 0.05).

 Correlation of environmental factors with the microbial community in <italic>H. armigera</italic>

Correlation analysis results showed that Enterobacter, Morganella, Lactococcus, Asaia, and Apibacter (top 5) were significantly spatially autocorrelated. In addition, the annual mean temperature and precipitation were also significantly spatially autocorrelated (Figure 5). A structural equation model (SEM) explored the relationships between H. armigera symbionts and five putative predictive variables (average temperature, precipitation, mean illumination, latitude, and longitude). The results showed that climate factors (temperature, precipitation, and illumination) significantly negatively correlated with space (latitude and longitude; R = −0.9503). The relative abundance of Apibacter significantly decreased with Asaia (−0.187) but increased with space (0.5546) and climate factors (0.4339), and it is also the case with Asaia. The proportion of Enterobacter significantly decreased with climate factors and space (−0.1722 and −0.9503, respectively), and it was also true with Morganella and Lactococcus. The proportion of Asaia significantly increased with space, climate factors, and other symbionts except Lactococcus. Meanwhile, the proportion of Lactococcus decreased with the other four symbionts, space, and climate factors.

Effects of various factors on the bacterial community compositions of different field populations of Helicoverpa armigera. The path diagram is based on the structural equation model (SEM) for environmental factors and microbial Bray–Curtis dissimilarity. The red solid and blue solid lines represent significant positive and negative paths, respectively. The arrow widths represent the strengths of these relationships. The R2 values of every box indicate the amount of variation in that variable explained by the input arrows. The numbers next to the arrows are unstandardized slopes.

 RT-qPCR analysis of the top 32 bacteria

To evaluate bacterial response to geographical locations, the top 32 bacterial genera (relative abundance) were selected for RT-qPCR. All the bacterial sequences were aligned to the NCBI database, and specific primers of the top 32 bacteria were designed by DNAMAN (Supplementary Table 4). The qPCR results showed that Serratia rarely appeared in samples, and Enterococcus appeared in the egg of H. armigera. Most OTUs were not significantly different between indoor and wild H. armigera larvae populations at the same stage. Copy numbers of OTU1387 (Bartonella), 890 (Citrobacter), 889 (Dysgonomonas), 39 (Acinetobacter), 876 (Serratia), 544 (Rosenbergiella), 497 (Enterobacter), 1,136 (Morganella), 754 (Lactococcus), 723 (Lactococcus), 935 (Bartonella), 559 (Enterococcus), 880 (Pseudomonas), 545 (Rosenbergiella), 961 (Acinetobacter), 561 (Orbus), 362 (Enterococcus), 1,326 (Enterococcus), and 563 (Pseudomonas) in the laboratory population were significantly different from those in the wild population, respectively (Figure 6). Apibacter, Lactococcus, Morganella, Bartonella, Orbus, Dysgonomonas, and Commensalibacter were not detected in the laboratory population, while they were widely present in the wild population. RT-qPCR results indicated significant differences in bacterium copy numbers among different geographical populations.

Thirty-two widespread bacterial species out of total bacterial 16S rDNA gene copies by quantitative PCR. Different lowercase labels above each group and * indicate significant differences (one-way ANOVA, Tukey's HSD test, P < 0.05) of group mean value.

 Discussion

In China, H. armigera is one of the most destructive pests in cotton (Wu and Guo, 2005), and Bt cotton has been an effective measure to control this pest since 1997. The percentage of Bt cotton increased from 50% in 2000 to almost 100% in 2004 in Yellow River, China (Zhang et al., 2011; Lu et al., 2012). The Yellow River region was once the largest cotton-planting area, and Henan, Shandong, and Hebei provinces were the main cotton-growing provinces in this region. In the Yellow River region, wheat is the main host plant of first-generation H. armigera larvae, and subsequent generations mainly feed on cotton, maize, peanuts, and soybeans (Wu and Guo, 2005).

Insects inhabit different ecological niches and face challenges such as adverse food sources, extreme environments, and threats from natural enemies. Interestingly, bacterial symbionts can promote insect fitness, and these beneficial effects of bacterial symbionts on insects' response to environmental stress can be attributed to strong selective overlap between the microbiota and their insect host from an evolutionary perspective (Douglas, 2015). The interaction between plants and phytophagous insects is driven by a hidden interplay between their microbiota (Di Lelio et al., 2023). Helicoverpa armigera is a worldwide migrant pest with highly variable distances traveled under different circumstances (Jones et al., 2019). Hence, the evaluations of the potential effects of bacterial symbionts on its host at different stages from different geographic regions are crucial for this destructive pest.

In this study, four phyla of bacteria (Proteobacteria, Firmicutes, Actinobacteria, and Bacteroidetes) and two genera of bacteria (Enterobacter and Enterococcus) dominated across all life stages of cotton bollworm, which was similar to several other lepidopteran pests (Chen et al., 2018). At the phylum level, there were no significant differences in the microbiota composition of the cotton bollworm across the whole life cycle from different geographical populations, which might be due to the simple intestinal structure of lepidopteran pests (Engel and Moran, 2013). These microbiota have been reported to play a vital role in insect physiology and metabolism (Jones et al., 2013). At the genus level, microbiota composition varied among developmental stages and geographic regions. Environmental habitat, food, and phylogeny of the host affected microbial community structure at different developmental stages (Yun et al., 2014; Staudacher et al., 2016; Zhao et al., 2019). Indeed, food can significantly alter the composition of host microbes (Wu et al., 2011). The insect–microbe interaction may improve the insect's ability to adapt to its environment and explore new niches. Microbial colonization and stable interactions are common in insects and can be mutualistic (Pais et al., 2018). The resistance of cotton bollworm to Bt is related to the contents of Enterobacter and Enterococcus (Broderick et al., 2006, 2009). Enterobacter can reduce the mortality of larvae of Plutella xylostella after exposure to Bt protein (Raymond et al., 2009). Typical intestinal bacteria such as Streptococcus spp. and Staphylococcus spp. have strong influences on the growth of B. thuringiensis in Homona magnanima larvae (Takatsuka and Kunimi, 2000). Hence, the high abundance of Enterococcus and Enterobacter in wild cotton bollworm populations probably has contributed to the increased resistance to Bt in recent 20 years. In addition, the diversity and variability of symbiotic bacteria in cotton bollworm might also explain its strong adaptability, wide distribution, and broad host range.

Enterobacter and Enterococcus are the top symbiotic bacteria in Lepidoptera (Broderick et al., 2004; Teh et al., 2016). Interestingly, these two bacteria are also found in the microbiota of the cotton bollworm eggs, indicating that they might be probiotics and can be vertically transmitted. The high abundance of Enterobacter in egg and pupa might be because neonates or moths often consume eggshells upon hatching or eclosion (Duplouy and Hornett, 2018), and Enterobacter contributes to the increase in the host fitness (Augustinos et al., 2015). In addition, lepidopteran pests such as H. armigera typically consume their chorions after hatching to obtain symbiotic bacteria to form microbiota rapidly to meet their own growth and development needs (Chen et al., 2018).

Insects' symbionts exhibit correlations with environmental changes across diverse geographical regions, responding to shifts in local ecological factors such as temperature, precipitation, altitude, and interactions with natural adversaries. The diversity of bacteria among psyllid species has been associated with their host plant affiliations and geographic dispersion, demonstrating robust preservation of Enterobacteriaceae secondary symbionts across distinct regions (Martoni et al., 2023). Within the realm of stingless bees, as the geographical span between colonies widens, the resemblance between them diminishes for both species, as observed in the study by Liu H. et al. (2023). The application of Mantel tests brought noteworthy associations between geographic proximity and symbiont communities in two specialized aphid species (Xu et al., 2020). An intriguing observation is the reliance of caterpillar microbiomes on soil microbiomes. This enables the transmission of plant-induced soil effects to subsequent insects feeding on different plants, as highlighted by Hannula et al. (2019). Notably, the highest bacterial genetic diversity in topsoil is found in temperate habitats, where the composition of microbial genes exhibits greater variation in response to environmental factors than mere geographical distances, an insight unveiled by Bahram et al. (2018). In our investigation, a positive correlation emerged between temperature and geographical distance. The resemblance in microbial composition within field populations, exemplified by Anyang–Qiuxian and Cangzhou–Binzhou, might find its explanation in climate-related latitude factors, offering a novel perspective on the dynamics of microbial communities.

Insect gut microbial diversity is determined by various factors including its environmental habitat, diet, developmental stage, and host phylogeny (Martinson et al., 2017; Chen et al., 2018). In this study, the highest microbial diversity was observed in the initial and final developmental stages (egg and adult). However, the bacterial diversity was significantly lower in the larva and pupa stages of H. armigera than in its egg and adult stages. This microbial diversity difference associated with development stages is common in lepidopterans since they undergo holometabolous metamorphosis, in which they entirely reshape the body structure and possess an unusually alkaline gut environment with an ~2-h rapid food intake (Paniagua Voirol et al., 2018). The insect population's susceptibility to plant resistance may be attributed to gut microbiota composition (Liu M. et al., 2023). Microorganisms are key players in the modulating of insect physiological pathways required to face the challenge when consuming plants (Di Lelio et al., 2023). The microorganism–host interactions further increase the complexity of host bio-evolution (Rosenberg and Zilber-Rosenberg, 2016). Hence, diet and development stages majorly influence bacterial community variability (Staudacher et al., 2016).

One study of systemic symbiosis in jewel wasps has demonstrated that the composition and function of the microbial communities are closely related to host evolution (Brooks et al., 2016). To complete their life cycle, the insects could utilize the energy supply from their symbiotic bacteria. Bacteria in the lepidopteran gut such as Enterobacter spp. can provide ROS-detoxifying enzymes such as superoxide dismutase or catalase (Xia et al., 2017). However, the roles of highly abundant Enterobacter in cotton bollworm remain to be further explored.

This study indicated that the dominant symbiotic bacteria determined the microbial community composition difference in cotton bollworm eggs of various geographical populations. Gene function prediction analysis showed that the relative abundance of microbiotas involved in translation, replication and repair, folding, classification, degradation, and transcription functions was significantly lower in the laboratory population than in the wild population, indicating that the field population of cotton bollworm could be better adapted to various environment stresses (such as plant defense, pathogenic bacteria, fungi, and virus, and insecticides), which benefited from the symbiotic bacteria. As a return, environmental factors also have a great impact on the microbial community composition. The richness and diversity of the microbial community in the laboratory population were significantly lower than in the field population. Furthermore, the closer the geographical distance was, the more similar the habitat climate and crop composition were, and the greater the similarity in microbial composition among hosts from geographically close areas. Similarly, the microbial communities of hosts from both field population and laboratory population have been investigated in several fruit flies such as Poecilus chalcites, Anopheles stephensi, and Ostrinia nubilalis, and microorganisms in the field population exhibited high diversity but low relative abundance with even distribution (Lehman et al., 2009; Rani et al., 2009; Belda et al., 2011; Chandler et al., 2011). The microbial community composition of Lymantria dispar varies significantly with various foods (Mason and Raffa, 2014). In the larvae of Heliothis virescens, the microbiota in the indoor and wild populations differed greatly and did not share any OTUs (Staudacher et al., 2016). A major reason for the low diversity of microbial communities in the laboratory insect populations might be that artificial foods often contain antibiotics that reduce and/or alter the diversity of microbes in the insects (Tang et al., 2012).

In addition, cotton bollworm is omnivorous and can come into contact with corn, wheat, and other crops in the wild (Wu and Guo, 2005), and the shell of eggs produced by field populations carries more microorganisms from the adult and can transfer to larvae due to the neonates' habit of feeding egg shell after hatching (Koyle et al., 2016; Kietz et al., 2018). Food can significantly influence host microbial diversity (Wu et al., 2011; Yun et al., 2014). Host plant species have been reported to have a greater impact on midgut bacterial diversity than the geographical location of the host plant, and the larval midgut symbionts were similar to those of food leaves (Priya et al., 2012). In this study, the microbial diversity difference in cotton bollworm eggs from different geographic populations may result from the combined effects of habitat environment and the food taken by the mother generation. Commensalibacter, Bartonella, Orbus, Dysgonomonas, Morganella, Pseudomonas, and other bacteria exhibited high abundance in cotton bollworms from different geographic populations. Thus, they may be the dominant bacteria contributing to the field populations' adaptation to the local environment. We speculated that the geographical environment could directly affect the diversity of microorganisms in insects and indirectly affect their progenies.

 Conclusion

We show that variations in food, developmental stages, and geographical locations all affect the composition of bacterial communities of this lepidopteran. The microbiota in H. armigera exhibits dynamic composition. The symbiont communities of the different stages of H. armigera differed greatly from each other. Close geographical proximity leads to more similar microbial composition, likely due to similar climate and crop patterns. Altogether, these discoveries highlight how developmental and environmental cues dramatically reshape symbiotic bacterial communities in H. armigera, a globally destructive invasive species. The findings illuminate the underpinnings of microbial interactions in this pest, laying the groundwork to unravel symbiont functional roles and devise novel microbiome-based management strategies that disrupt associations optimizing its fitness and adaptability. Much remains to be discovered about the tripartite cross-talk between host, microbiota, and habitat sustaining the success of this agricultural nemesis.

 Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://www.ncbi.nlm.nih.gov/, PRJNA591375.

 Author contributions

JC, JL, and JJ carried out lab work. CZ and JC participated in data analysis and drafted the manuscript. LW, KZ, XZ, and DL participated in data analysis. CZ conceived of the study, designed the study, coordinated the study, and helped draft the manuscript. All authors gave final approval for publication.

 Funding

This research was supported by the National Key R&D Program of China (2022YFD1400300) and the Agricultural Science and Technology Innovation Program of the Chinese Academy of Agricultural Sciences.

 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.2023.1251627/full#supplementary-material

Rarefaction curves of bacterial communities across life cycle (A) and from different geographic population (B) of H. armigera.

Summary of high-throughput sequencing read analysis, bacterial community diversity richness (OTUs, 97%), sample coverage (Good's coverage), diversity index (Shannon, ACE, Shannon, Simpson), and estimated OTU richness (Chao1) for community diversity analyses of 27 samples from H. armigera in different development stages.

Summary of high-throughput sequencing read analysis, bacterial community diversity and richness index in different geographic population of H. armigera.

Number of different taxonomic categories of bacterial communities across whole H. armigera life cycle.

Primer sequences used in this experiment.

 References Augustinos A. A. Kyritsis G. A. Papadopoulos N. T. Abd-Alla A. M. Cáceres C. Bourtzis K. . (2015). Exploitation of the medfly gut microbiota for the enhancement of sterile insect technique: use of Enterobacter sp. in larval diet-based probiotic applications. PLoS ONE 10, e0136459. 10.1371/journal.pone.013645926325068 Bahram M. Hildebrand F. Forslund S. K. Anderson J. L. Soudzilovskaia N. A. Bodegom P. M. . (2018). Structure and function of the global topsoil microbiome. Nature 560, 233–237. 10.1038/s41586-018-0386-630069051 Belda E. Pedrola L. Pereto J. Martinez-Blanch J. F. Montagud A. Navarro E. . (2011). Microbial diversity in the midguts of field and lab-reared populations of the European corn borer Ostrinia nubilalis. PLoS ONE 6, e21751. 10.1371/journal.pone.002175121738787 Broderick N. A. Raffa K. F. Goodman R. M. Handelsman J. (2004). Census of the bacterial community of the gypsy moth larval midgut by using culturing and culture-independent methods. Appl. Environ. Microbiol. 70, 293–300. 10.1128/AEM.70.1.293-300.200414711655 Broderick N. A. Raffa K. F. Handelsman J. (2006). Midgut bacteria required for Bacillus thuringiensis insecticidal activity. Proc. Natl. Acad. Sci. U. S. A. 103, 15196–15199. 10.1073/pnas.060486510317005725 Broderick N. A. Robinson C. J. McMahon M. D. Holt J. Handelsman J. Raffa K. F. . (2009). Contributions of gut bacteria to Bacillus thuringiensis-induced mortality vary across a range of Lepidoptera. BMC Biol. 7, 11. 10.1186/1741-7007-7-1119261175 Brooks A. W. Kohl K. D. Brucker R. M. van Opstal E. J. Bordenstein S. R. (2016). Phylosymbiosis: relationships and functional effects of microbial communities across host evolutionary history. PLoS Biol. 14, e2000225. 10.1371/journal.pbio.200022528068336 Caporaso J. G. Kuczynski J. Stombaugh J. Bittinger K. Bushman F. D. Costello E. K. . (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. 10.1038/nmeth.f.30320383131 Chandler J. A. Lang J. M. Bhatnagar S. Eisen J. A. Kopp A. (2011). Bacterial communities of diverse Drosophila species: ecological context of a host-microbe model system. PLoS Genet. 7, e1002272. 10.1371/journal.pgen.100227221966276 Chen B. Du K. Sun C. Vimalanathan A. Liang X. Li Y. . (2018). Gut bacterial and fungal communities of the domesticated silkworm (Bombyx mori) and wild mulberry-feeding relatives. ISME. J. 12, 2252–2262. 10.1038/s41396-018-0174-129895989 Cui J. J. Luo J. Y. Van der Werf W. Ma Y. Xia J. Y. (2011). Effect of pyramiding Bt and CpTI genes on resistance of cotton to Helicoverpa armigera (Lepidoptera: Noctuidae) under laboratory and field conditions. J. Econ. Entomol. 104, 673–684. 10.1603/EC0922821510221 Di Lelio I. Forni G. Magoga G. Brunetti M. Bruno D. Becchimanzi A. . (2023). A soil fungus confers plant resistance against a phytophagous insect by disrupting the symbiotic role of its gut microbiota. Proc. Nat. Acad. Sci. 120, e2216922120. 10.1073/pnas.221692212036848561 Douglas A. E. (2015). Multiorganismal insects: diversity and function of resident microorganisms. Annu. Rev. Entomol. 60, 17–34. 10.1146/annurev-ento-010814-02082225341109 Downes S. Kriticos D. Parry H. Paull C. Schellhorn N. Zalucki M. P. . (2017). A perspective on management of Helicoverpa armigera: transgenic Bt cotton, IPM, and landscapes. Pest Manag. Sci. 73, 485–492. 10.1002/ps.446127753247 Duplais C. Sarou-Kanian V. Massiot D. Hassan A. Perrone B. Estevez Y. . (2021). Gut bacteria are essential for normal cuticle development in herbivorous turtle ants. Nat. Commun. 12, 676. 10.1038/s41467-021-21065-y33514729 Duplouy A. Hornett E. A. (2018). Uncovering the hidden players in Lepidoptera biology: the heritable microbial endosymbionts. PeerJ 6, e4629. 10.7717/peerj.462929761037 Edgar R. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996–998. 10.1038/nmeth.260423955772 Engel P. Moran N. A. (2013). The gut microbiota of insects - diversity in structure and function. FEMS Microbiol. Rev. 37, 699–735. 10.1111/1574-6976.1202523692388 Hamady M. Lozupone C. Knight R. (2010). Fast UniFrac: facilitating high-throughput phylogenetic analyses of microbial communities including analysis of pyrosequencing and PhyloChip data. ISME. J. 4, 17–27. 10.1038/ismej.2009.9719710709 Hammer T. J. Janzen D. H. Hallwachs W. Jaffe S. P. Fierer N. (2017). Caterpillars lack a resident gut microbiome. Proc. Natl. Acad. Sci. USA. 114, 9641–9646. 10.1073/pnas.170718611428830993 Hannula S. E. Zhu F. Heinen R. Bezemer T. M. (2019). Foliar-feeding insects acquire microbiomes from the soil rather than the host plant. Nat. Commun. 10, 1254. 10.1038/s41467-019-09284-w30890706 Jia Y. Jin S. Hu K. Geng L. Han C. Kang R. . (2021). Gut microbiome modulates Drosophila aggression through octopamine signaling. Nat. Commun. 12, 2698. 10.1038/s41467-021-23041-y33976215 Jiang R. X. Shang F. Jiang H. B. Dou W. Cernava T. Wang J. J. . (2023). Environmental factors affect the bacterial community in Diaphorina citri, an important vector of “Candidatus Liberibacter asiaticus”. Microbiol. Spectr. 11, e0529822. 10.1128/spectrum.05298-2236975996 Jones C. M. Parry H. Tay W. T. Reynolds D. R. Chapman J. W. (2019). Movement ecology of pest helicoverpa: implications for ongoing spread. Annu. Rev. Entomol. 64, 277–295. 10.1146/annurev-ento-011118-11195930296859 Jones R. T. Sanchez L. G. Fierer N. (2013). A cross-taxon analysis of insect-associated bacterial diversity. PLoS ONE 8, e61218. 10.1371/journal.pone.006121823613815 Kietz C. Pollari V. Meinander A. (2018). Generating germ-free Drosophila to study gut-microbe interactions: protocol to rear Drosophila under axenic conditions. Curr. Protoc. Toxicol. 77, e52. 10.1002/cptx.5229933523 Koyle M. L. Veloz M. Judd A. M. Wong A. C. Newell P. D. Douglas A. E. . (2016). Rearing the fruit fly Drosophila melanogaster under axenic and gnotobiotic conditions. J. Vis. Exp. 113, e54219. 10.3791/54219-v27500374 Lehman R. M. Lundgren J. G. Petzke L. M. (2009). Bacterial communities associated with the digestive tract of the predatory ground beetle, Poecilus chalcites, and their modification by laboratory rearing and antibiotic treatment. Microb. Ecol. 57, 349–358. 10.1007/s00248-008-9415-618587608 Liu H. Hall M. A. Brettell L. E. Wang J. Halcroft M. Nacko S. . (2023). Microbial diversity in stingless bee gut is linked to host wing size and influenced by the environment. J. Invertebr. Pathol. 198, 107909. 10.1016/j.jip.2023.10790936889457 Liu M. Hong G. Li H. Bing X. Chen Y. Jing X. . (2023). Sakuranetin protects rice from brown planthopper attack by depleting its beneficial endosymbionts. Proc. Natl. Acad. Sci. USA. 120, e2305007120. 10.1073/pnas.230500712037256931 Lu Y. Wu K. Jiang Y. Guo Y. Desneux N. (2012). Widespread adoption of Bt cotton and insecticide decrease promotes biocontrol services. Nature 487, 362–365. 10.1038/nature1115322722864 Lu Y. Wyckhuys K. A. G. Yang L. Liu B. Zeng J. Jiang Y. . (2022). Bt cotton area contraction drives regional pest resurgence, crop loss, and pesticide use. Plant Biotechnol. J. 20, 390–398. 10.1111/pbi.1372134626524 Martinson V. G. Douglas A. E. Jaenike J. (2017). Community structure of the gut microbiota in sympatric species of wild Drosophila. Ecol. Lett. 20, 629–639. 10.1111/ele.1276128371064 Martoni F. Bulman S. R. Piper A. M. Pitman A. Taylor G. S. Armstrong K. F. . (2023). Insect phylogeny structures the bacterial communities in the microbiome of psyllids (Hemiptera: Psylloidea) in Aotearoa New Zealand. PLoS ONE 18, e0285587. 10.1371/journal.pone.028558737186593 Mason C. J. Raffa K. F. (2014). Acquisition and structuring of midgut bacterial communities in gypsy moth (Lepidoptera: Erebidae) larvae. Environ. Entomol. 43, 595–604. 10.1603/EN1403124780292 Mason C. J. Ray S. Shikano I. Peiffer M. Jones A. G. Luthe D. S. . (2019). Plant defenses interact with insect enteric bacteria by initiating a leaky gut syndrome. Proc. Natl. Acad. Sci. USA. 116, 15991–15996. 10.1073/pnas.190874811631332013 Oksanen J. Blanchet F. G. Kindt R. Legendre P. O'Hara R. G. Simpson G. (2009). Vegan: Community Ecology Package. R package version 1. 16–33. Available online at: http://CRAN.R-project.org/package=vegan Pais I. S. Valente R. S. Sporniak M. Teixeira L. (2018). Drosophila melanogaster establishes a species-specific mutualistic interaction with stable gut-colonizing bacteria. PLoS Biol. 16, e2005710. 10.1371/journal.pbio.200571029975680 Paniagua Voirol L. R. Frago E. Kaltenpoth M. Hilker M. Fatouros N. E. (2018). Bacterial symbionts in lepidoptera: their diversity, transmission, and impact on the host. Front. Microbiol. 9, 556. 10.3389/fmicb.2018.0055629636736 Priya N. G. Ojha A. Kajla M. K. Raj A. Rajagopal R. (2012). Host plant induced variation in gut bacteria of Helicoverpa armigera. PLoS ONE 7, e30768. 10.1371/journal.pone.003076822292034 Qin Q. (2015). Rearing the cotton bollworm Helicoverpa armigera in the laboratory. Chin. J. Appl. Entomol. 52, 486–491. 10.7679/j.issn.20951353.2015.055 Rani A. Sharma A. Rajagopal R. Adak T. Bhatnagar R. K. (2009). Bacterial diversity analysis of larvae and adult midgut microflora using culture-dependent and culture-independent methods in lab-reared and field-collected Anopheles stephensi-an Asian malarial vector. BMC Microbiol. 9, 96. 10.1186/1471-2180-9-9619450290 Raymond B. Johnston P. R. Wright D. J. Ellis R. J. Crickmore N. Bonsall M. B. . (2009). A mid-gut microbiota is not required for the pathogenicity of Bacillus thuringiensis to diamondback moth larvae. Environ. Microbiol. 11, 2556–2563. 10.1111/j.1462-2920.2009.01980.x19555371 Rosenberg E. Zilber-Rosenberg I. (2016). Microbes drive evolution of animals and plants: the hologenome concept. MBio 7, e01395. 10.1128/mBio.01395-1527034283 Schloss P. D. Westcott S. L. Ryabin T. Hall J. R. Hartmann M. Hollister E. B. . (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–7541. 10.1128/AEM.01541-0919801464 Shukla S. P. Plata C. Reichelt M. Steiger S. Heckel D. G. Kaltenpoth M. . (2018). Microbiome-assisted carrion preservation aids larval development in a burying beetle. Proc. Natl. Acad. Sci. USA. 115, 11274–11279. 10.1073/pnas.181280811530322931 Staudacher H. Kaltenpoth M. Breeuwer J. A. Menken S. B. Heckel D. G. Groot A. T. . (2016). Variability of bacterial communities in the moth Heliothis virescens indicates transient association with the host. PLoS ONE 11, e0154514. 10.1371/journal.pone.015451427139886 Takatsuka J. Kunimi Y. (2000). Intestinal bacteria affect growth of Bacillus thuringiensis in larvae of the oriental tea tortrix, Homona magnanima diakonoff (Lepidoptera: tortricidae). J. Invertebr. Pathol. 76, 222–226. 10.1006/jipa.2000.497311023751 Tang X. Freitak D. Vogel H. Ping L. Shao Y. Cordero E. A. . (2012). Complexity and variability of gut commensal microbiota in polyphagous lepidopteran larvae. PLoS ONE 7, e36978. 10.1371/journal.pone.003697822815679 Teh B. S. Apel J. Shao Y. Boland W. (2016). Colonization of the intestinal tract of the polyphagous pest Spodoptera littoralis with the GFP-tagged indigenous gut bacterium Enterococcus mundtii. Front. Microbiol. 7, 928. 10.3389/fmicb.2016.0092827379058 Tenenhaus M. Vinzi V. E. Chatelin Y.-M. Lauro C. (2005). PLS path modeling. Comput. Stat. Data Anal. 48, 159–205. 10.1016/j.csda.2004.03.005 Wilson L. J. Whitehouse M. E. A. Herron G. A. (2018). The management of insect pests in australian cotton: an evolving story. Annu. Rev. Entomol. 63, 215–237. 10.1146/annurev-ento-020117-04343229324044 Wu G. D. Chen J. Hoffmann C. Bittinger K. Chen Y. Y. Keilbaugh S. A. . (2011). Linking long-term dietary patterns with gut microbial enterotypes. Science 334, 105–108. 10.1126/science.120834421885731 Wu K. M. Guo Y. Y. (2005). The evolution of cotton pest management practices in China. Annu. Rev. Entomol. 50, 31–52. 10.1146/annurev.ento.50.071803.13034915355239 Xia X. Gurr G. M. Vasseur L. Zheng D. Zhong H. Qin B. . (2017). Metagenomic sequencing of diamondback moth gut microbiome unveils key holobiont adaptations for herbivory. Front. Microbiol. 8, 663. 10.3389/fmicb.2017.0066328491055 Xu N. Tan G. C. Wang H. Y. Gai X. P. (2016). Effect of biochar additions to soil on nitrogen leaching, microbial biomass and bacterial community structure. Eur. J. Soil Biol. 74, 1–8. 10.1016/j.ejsobi.2016.02.00433124402 Xu P. Liu Y. Graham R. I. Wilson K. Wu K. (2014). Densovirus is a mutualistic symbiont of a global crop pest (Helicoverpa armigera) and protects against a baculovirus and Bt biopesticide. PLoS Pathog. 10, e1004490. 10.1371/journal.ppat.100449025357125 Xu S. Jiang L. Qiao G. Chen J. (2020). The bacterial flora associated with the polyphagous aphid Aphis gossypii Glover (Hemiptera: Aphididae) is strongly affected by host plants. Microb. Ecol. 79, 971–984. 10.1007/s00248-019-01435-231802184 Yuan C. Xing L. Wang M. Hu Z. Zou Z. (2021). Microbiota modulates gut immunity and promotes baculovirus infection in Helicoverpa armigera. Insect Sci. 28, 1766–1779. 10.1111/1744-7917.1289433463036 Yun J. H. Roh S. W. Whon T. W. Jung M. J. Kim M. S. Park D. S. . (2014). Insect gut bacterial diversity determined by environmental habitat, diet, developmental stage, and phylogeny of host. Appl. Environ. Microbiol. 80, 5254–5264. 10.1128/AEM.01226-1424928884 Zhang H. N. Yin W. Zhao J. Jin L. Yang Y. H. Wu S. W. . (2011). Early warning of cotton bollworm resistance associated with intensive planting of Bt cotton in China. PLoS ONE 6, e22874. 10.1371/journal.pone.002287421857961 Zhang S. Luo J. Jiang W. Wu L. Zhang L. Ji J. . (2019). Response of the bacterial community of Propylea japonica (Thunberg) to Cry2Ab protein. Environ. Pollut. 254, 113063. 10.1016/j.envpol.2019.11306331454585 Zhao C. Wu L. Luo J. Niu L. Wang C. Zhu X. . (2020). Bt, not a threat to Propylea japonica. Front. Physiol. 11, 758. 10.3389/fphys.2020.0075832903558 Zhao C. Zhao H. Zhang S. Luo J. Zhu X. Wang L. . (2019). The Developmental stage symbionts of the pea aphid-feeding Chrysoperla sinica (Tjeder). Front. Microbiol. 10, 2454. 10.3389/fmicb.2019.0245431736900 Zhu Z. Liu Y. Hu H. Wang G. H. (2022). Nasonia-microbiome associations: a model for evolutionary hologenomics research. Trends Parasitol. 39, 101–112. 10.1016/j.pt.2022.11.00536496327