- Research
- Open access
- Published:
Gut dysbacteriosis induces expression differences in the adult head transcriptome of Spodoptera frugiperda in a sex-specific manner
BMC Microbiology volume 23, Article number: 388 (2023)
Abstract
Mounting evidence indicates that the gut microbiota influences the neurodevelopment and behavior of insects through the gut-brain axis. However, it is currently unclear whether the gut microbiota affect the head profiles and immune pathway in pests. Here, we find that gut bacteria is essential for the immune and neural development of adult Spodoptera frugiperda, which is an extremely destructive agricultural pest worldwide. 16 S rRNA sequencing analysis showed that antibiotics exposure significantly disturbed the composition and diversity of gut bacteria. Further transcriptomic analysis revealed that the adult head transcripts were greatly affected by gut dysbacteriosis, and differently expression genes critical for brain and neural development including A4galt, Tret1, nsun4, Galt, Mitofilin, SLC2A3, snk, GABRB3, Oamb and SLC6A1 were substantially repressed. Interestingly, the dysbacteriosis caused sex-specific differences in immune response. The mRNA levels of pll (serine/threonine protein kinase Pelle), PGRP (peptidoglycan-sensing receptor), CECA (cecropin A) and CECB (cecropin B) involved in Toll and Imd signaling pathway were drastically decreased in treated male adults’ heads but not in female adults; however, genes of HIVEP2, ZNF131, inducible zinc finger protein 1-like and zinc finger protein 99-like encoding zinc-finger antiviral protein (ZAP) involved in the interferon (IFNα/β) pathway were significantly inhibited in treated female adults’ heads. Collectively, these results demonstrate that gut microbiota may regulate head transcription and impact the S. frugiperda adults’ heads through the immune pathway in a sex-specific manner. Our finding highlights the relationship between the gut microbiota and head immune systems of S. frugiperda adults, which is an astonishing similarity with the discoveries of other animals. Therefore, this is the basis for further research to understand the interactions between hosts and microorganisms via the gut-brain axis in S. frugiperda and other insects.
Background
The bidirectional interactions of the nervous and immune systems have been widely demonstrated in vertebrates [1, 2]. Some compounds released by immune cells cause changes in neuronal activities [3]. In some cases, it can lead to changes in behavioral coordination. On the other hand, the nervous system changes immune function through the endocrine system [4]. Therefore, the immune and nervous systems are involved in the crosstalk related to function and homeostasis. It is the main function of the immune system, serving as a sensory organ and sending signals to the central nervous system, monitoring immune status and microbial challenges [5]. More studies have suggested that interactions between the nervous and immune systems are no longer unique to vertebrates. The physiological behavioral phenomena in insects are similar to those in vertebrates [6]. The stimulation of the immune response can cause insect anorexia, changes in behavior and reproduction, and a decrease in learning abilities [7, 8].
Recently, many studies have demonstrated that gut microbes have underestimated potential to contribute to brain development and function. The metabolites short-chain fatty acids (SCFAs) and others derived from the gut microbiome regulate the behavior of immune cells throughout the body, including those in the brain [9]. Neuroactive compounds produced by gut bacteria influence brain function and behavior [10], leading to nervous disorders [11, 12]. Two-way regulation of the gut microbiota-brain axis in mammalian models has revealed the contribution of microbiota to the occurrence of neurodegenerative diseases [13, 14] and emotional regulation, including anxiety and depression [15]. Recent studies have also shown a link between gut dysbiosis and social barriers, such as autism spectrum disorder (ASD) [16] and schizophrenia [17]. However, studies on the connections between the gut microbiota and brain have mainly focused on vertebrate model organisms, such as mice and rats [18]. As one of the most diverse species in the world, insects possess more diverse gut microorganisms affecting development, life span, reproduction and immune responses [19]. The insect brain controls the chemical sense, physiological state and behavioral ability. The gut microbiota has shown profound effects on odorant profiles, olfactory behavior [20], and neurophysiological development [21] in insects. For example, the lower termite Reticulitermes speratus can easily recognize and attack invaders that are colonized by foreign gut bacteria, promoting unfamiliar scents [22]. Gut dysbacteriosis interferes with the hive behaviors of nurse bees, affecting neurotransmitter concentrations, protein levels in the brain and circulating metabolic profiles [23]. Compared with fruit flies fed traditional gut microbiota, axenic Drosophila have worse learning and memory in aversive phototactic assays [24]. To date, the comprehensive understanding of these interactions is still limited due to only a few model organisms being used, resulting in the inability to truly understand the proximate mechanisms involved in this phenomenon.
In Lepidoptera, Spodoptera frugiperda (Lepidoptera: Noctuidae) is a typical agricultural pest that feeds on more than 300 host plant species and has strong adaptability, flight capabilities and reproductive capability [25, 26]. Due to the influence of the environment, development stage and host plants, S. frugiperda has diverse gut microbial species. The phyla Firmicutes and Bacteroidetes dominate across S. frugiperda life stages, and the genera Enterococcus, Enterobacter, Acinetobacter, Bacillus and Klebsiella are significantly enriched and more frequent in different samples, especially Enterococcus, which is a core species [27,28,29,30,31]. There are only a few reports on the specific regulation and molecular mechanism of the gut microbiota on S. frugiperda physiological functions. Antibiotic exposure could affect energy and metabolic homeostasis in S. frugiperda [32], and the larval food intake and body weight of S. frugiperda decreased, and the larval stage was prolonged due to changes in the gut microbiota [33]. The field-collected larvae showed greater potential for insecticides than those from laboratory-selected resistant strains, which is an important reason for S. frugiperda resistance to various pesticides in the field [34]. Therefore, further studies are necessary to explore the functional roles of the gut microbiota in S. frugiperda.
To fully appreciate the role of bacterial symbionts in the evolution and brain behavior of S. frugiperda, the potential impact of gut microbiota on the head transcriptome was investigated in this study. First, S. frugiperda with gut dysbacteriosis was established, and RNA-seq results proved that gut bacteria could impact the head transcriptional program, especially immune pathways, and the immune response to gut microbiome changes occurred in a sex-specific manner. Finally, the expression levels of genes related to the immune response in the brain were determined by qRT‒PCR. This information will help to further investigate the roles of gut microbes in the behavior of S. frugiperda and may be a feasible strategy for future insect biological management by adjusting insect behavioral traits.
Materials and methods
Insect rearing and antibiotic treatment
The population of S. frugiperda originally collected from corn fields in Dali city (Yunnan Province, China) was established and maintained in the laboratory. After hatching from eggs, S. frugiperda larvae were fed an artificial diet as previously described [35]. All larvae and adults were reared under the conditions of temperature 26 ± 1 °C, relative humidity 65%±5%, and a photoperiod of 14:10 h (light: dark). The emerged adults were supplied with 10% honey solution. For the same batch of egg masses, parts of the egg masses were raised routinely, while others were collected 48 h after laying and dechorionated for 4 min in 4% formaldehyde solution, immersed in 8% sodium hypochlorite solution (containing 4% sodium hydroxide) for 3 min, and rinsed twice with sterile water [36]. The eggs were then transferred to beef-protein medium and allowed to develop. After confirmation of free bacterial contamination, the newly hatched larvae were transferred to a sterile artificial diet that was mixed with a combination of antibiotics (gentamycin, penicillin, streptomycin, ciprofloxacin hydrochloride, rifampicin and vancomycin, 100 µg/mL each). To maintain sterile conditions, all culture tubes and tools underwent autoclave sterilization, and all manipulations were performed in a biosafety cabinet. The pupae that emerged from the treated larvae were transferred to sterile tissue culture bottles. Adults were fed a 10% sterile honey solution mixed with a combination of antibiotics. Eggs were adhered to sterile spawning papers in bottles. Larvae and adults of the next generation continued to be reared through the same process. Then, the second generation of treated adults and untreated adults were collected for subsequent experiments.
Extraction of gut samples and 16 S rRNA sequencing
Conventional female adult (C-female) and male adult (C-male) gut samples, and treated female adult (T-female) and male adult (T-male) gut samples (3 days old) were collected for microbiome analysis. On a super clean bench, all procedures were performed with sterile scalpel tweezers on ice. Gut samples were rinsed in 75% ethanol for 30 s and rinsed again in sterile 0.1mol PBS (pH 7.2). A total of 15 moths from each group were used for three biological replicates. Using HiPure Soil DNA Kits (Magen, Guangzhou, China), bacterial DNA was extracted according to the manufacturer’s protocols. The V3-V4 region of the 16S rRNA gene was amplified using the primers 341F (5’-CCTACGGGNGGCWGCAG-3’) and 806R (5’-GGACTACHVGGGTATCTAAT-3’) [37]. Amplicons were extracted from 2% agarose gels and purified using an AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, U.S.) according to the manufacturer’s instructions and quantified using an ABI StepOnePlus Real-Time PCR System (Life Technologies, Foster City, USA).
The purified amplicons were pooled in equimolar amounts and paired-end sequenced (PE250) on an Illumina NovaSeq 6000 platform. DNA sequences for this experiment are available at the NCBI Sequence Read Archive (SRA) database. Raw 16 S rRNA fastq files were quality filtered according to the standard protocols described previously [38]. OTUs were clustered using a cutoff of 97% similarity with the UPARSE (version 9.2.64) pipeline [39] and clustered into organisms through the naïve Bayesian model using the RDP classifier (version 2.2) [40] based on the SILVA database (version 132) [41] and UNITE database (version 8.0) [42] with a confidence threshold value of 0.8. The stacked bar plot of the community composition at the phylum or genus level was visualized in R project ggplot2 package (version 2.2.1) [43]. Chao1 and Shannon index were all measured based on OTUs at the species level with QIIME (version 1.9.1) [44]. Principal coordinate analysis (PCoA) of weighted UniFrac distances was used to compare bacterial diversity between groups with the Vegan package of R program (version 2.5.3) [45]. Linear discriminant analysis coupled with effect size (LEfSe) analysis (LDA score > 2 and P < 0.05) was used to identify distinguishing genus of each group [46].
Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis was performed using Tax4Fun (version 1.0) [47]. Functional differences between groups were calculated by Welch’s t test and Kruskal‒Wallis H test or Tukey’s HSD test in the Vegan package of R program (version 2.5.3) [48].
RNA extraction and transcriptome sequencing
Virgin adult head samples were collected at 72 h after emergence for transcriptome analysis, and all procedures were performed with sterile tweezers and petri dishes on ice. Moreover, corresponding brain and antenna samples were also extracted for subsequent qRT-PCR analysis. Samples were rinsed in 75% ethanol for 30 s and rinsed again in sterile 0.1 mol PBS (pH 7.2), followed by tissue dissection. Three biological replicates (total 15 moths per group) of frozen samples were ground in liquid nitrogen. Using the TRIzol reagent kit (Invitrogen, Carlsbad, CA, USA), total RNA was isolated according to the manufacturer’s protocol. RNA quality was assessed by RNase-free agarose gel electrophoresis in an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA).
RNA sequencing analysis
After total RNA was extracted, eukaryotic mRNA was enriched by Oligo (dT) beads. Then the enriched mRNA was fragmented into short fragments using fragmentation buffer and reversely transcribed into cDNA by using NEBNext Ultra RNA Library Prep Kit for Illumina (NEB #7530, New England Biolabs, Ipswich, MA, USA).The purified double-stranded cDNA fragments were end repaired, a base added, and ligated to Illumina sequencing adapters. The ligation reaction was purified with the AMPure XP Beads(1.0X).And polymerase chain reaction (PCR) amplified.The resulting cDNA library was sequenced using Illumina Novaseq6000 by Gene Denovo Biotechnology Co. (Guangzhou, China). Reads were further filtered by fastp [49] (version 0.18.0) to get high quality clean reads and short reads alignment tool Bowtie2 [50] (version 2.2.8) was used for mapping reads to ribosome. Paired-end clean reads were further used in assembly and mapped to the reference assembly genome ZJU_Sfru_1.0 (Assembly accession: GCF_011064685.1) using HISAT2. 2.4 [51]. The mapped reads of each sample were assembled using StringTie v1.3.1 [52] in a reference-based approach. For each transcription region, a FPKM (fragment per kilobase of ranscript per million mapped reads) value was calculated to quantify its expression abundance and variations using RSEM [53] software. RNAs differential expression analysis was performed by DESeq2 [54] software between two different groups. Analysis of differential expression between two different groups was performed using DESeq2 software [55]. Transcripts with a false discovery rate (FDR) below 0.05 and absolute fold change ≥ 2 were considered differentially expressed.
Gene Ontology (GO) and KEGG pathway enrichment analyses were conducted to analyze the biological function. All DEGs were mapped to GO terms in the Gene Ontology database (http://www.geneontology.org/), gene numbers were calculated for every term, significantly enriched GO terms in DEGs comparing to the genome background were defined by hypergeometric test. Significantly enriched GO terms and pathways of the genes were defined by a hypergeometric test and a threshold of FDR less than 0.05 [56]. In addition, gene set enrichment analysis (GSEA) and MSigDB software programs [57] were used to identify specific GO terms/pathways and show significant differences between the two groups. Enrichment scores and p values were calculated with the default parameters.
Quantitative real-time PCR
Quantitative real-time PCR (qRT-PCR) was performed using SYBR Green I Master Mix in a standard 96-well block to measure gene expression levels. The relative expression level of genes was determined by the 2−ΔΔCt method with the gene GAPDH (glyceraldehyde-3-phosphate dehydrogenase, GenBank ID: 118,271,716) as an internal reference. Three biological replicates were performed in the same sample. The primers used in this study are listed in Table S1.
Statistical analyses
Statistical analysis between the control and treated groups was conducted using Student’s t test for unpaired comparisons. When p < 0.05, the difference was considered statistically significant. All statistics were completed using GraphPad 5.0. All data are displayed as the mean ± standard error (SEM). Welch’s t test in R Project was used to analyze the sequencing data of gut microbiota.
Results
Effect of antibiotic treatment on gut microbiota composition and function
To test whether the gut microbiota impacts the head transcriptome of S. frugiperda adults, the amplicons of 16 S rRNA genes isolated from the gut microbiota of 12 adult samples were sequenced first. A total of 1,578,994 raw reads with an average length of 472.58 bp were obtained. The reads with sequence similarity greater than 97% at the species level were clustered into 2656 OTUs. Among them, OTUs were distributed in four groups, as follows: C-female (308 OTUs), C-male (532 OTUs), T-female (903 OTUs), and T-female (913 OTUs) (Table S2).
Composition analysis of the adult gut microbiome demonstrated that the microbial species in the gut were similar between virgin female and male adults (Fig. 1A-B). Proteobacteria (female, 50.52%; male, 45.07%) and Firmicutes (female, 43.62%; male, 51.00%) (average value across all of the samples) were the dominant phyla in the control adult guts, followed by Bacteroidetes and Actinobacteria (Fig. 1A). In antibiotic-treated adults, a decrease in Firmicutes (female, 5.73%; male, 8.90%) and an increase in Bacteroidetes (female, 22.59%; male, 25.46%) were observed (Fig. 1A). Similarly, the relative abundances of Enterococcus (female, 42.48%; male, 52.53%) and Enterobacter (female, 31.48%; male, 28.97%) in the control groups were dramatically greater than those in the treated groups (female, 3.62%; male, 2.13%; female, 8.74%; male, 3.67%) (Fig. 1B). In addition, Ralstonia (female, 29.82%; male, 37.80%) and Sediminibacterium (female, 19.90%; male, 25.30%) were shown to be the most obviously increased genera (Fig. 1B). These results suggest that antibiotic treatment results in dysbacteriosis of the adult gut. Furthermore, the antibiotic-treated adults showed a significant increase in bacterial α-diversity, as evidenced by the Chao1 and Shannon indices (Fig. 1C). The weighted UniFrac PCoA showed two independent clusters with the bacterial composition between conventional and antibiotic-treated groups, profound along the PCO1 axis (reaching 71.17% of overall variation), while no significant difference was found between the female and male adults (Fig. 1D). To identify the dominant bacterial species in these two groups, linear discriminant effect size (LEfSe) analysis was performed, which showed 21 differential bacterial taxa in females and 31 differential bacterial taxa in males (Fig. 1E-F). Firmicutes and Enterococcus had higher relative abundances in both the untreated females and males. However, the most differentially abundant bacterial taxa in antibiotic-treated female adults were Proteobacteria and Rubrobacter, while the key bacterial taxa Bacteroidetes and Ralstonia were more abundant in the antibiotic-treated male groups (Fig. 1E-F). In brief, these results suggest that antibiotic exposure contributes to gut microbiota dysbiosis.
Furthermore, the Tax4Fun program was used to predict gut microbiota functions. As shown in Fig. 2A, the top 20 pathways at the third level of KEGG, such as ABC transporters, purine metabolism, pyrimidine metabolism, amino sugar and nucleotide sugar metabolism, peptidoglycan biosynthesis, fructose and mannose metabolism, and phosphotransferase system (PTS), were enriched differently between the treated female and the control female adults (Fig. 2A). The top 20 pathways such as amino sugar and nucleotide sugar metabolism, arginine and proline metabolism, porphyrin and chlorophyll metabolism, oxidative phosphorylation, fructose and mannose metabolism, phosphotransferase system (PTS), glycine, serine and threonine metabolism were significantly affected in male adults after antibiotic exposure (Fig. 2B). In particular, the most related KEGG pathways were metabolism which were less enriched in all treated groups (Fig. 2A-B). These data reveals that gut microbiota affects host through their biological functions regulated by the genes of gut microbiota, causing female and male adults to have different metabolic responses.
Sex differences in the head transcriptome between female and male adults
To further investigate whether there was sexual dimorphism in the female and male heads of S. frugiperda, a comparative transcriptomic analysis of 3-day-old female and male adult heads was performed. There were 34 differentially expressed genes (DEGs) (|log2-fold change|>2, false discovery rate [FDR] < 0.05, Table S3) identified, of which 12 DEGs were upregulated and 22 DEGs were downregulated in the comparison group of conventional male heads and female heads (Fig. 3A). Furthermore, GO enrichment analysis indicated that the DEGs were mainly enriched in odorant binding, extracellular region, and mRNA cap binding complex (Fig. 3B). Among them, four genes encoding odorant-binding proteins (OBPs) (ncbi_118265739, ncbi_118267992, ncbi_118275388, ncbi_118279193) exhibited significantly reduced expression in the male head (Fig. 3B). Genes related to the RNA cap binding complex were also decreased in male heads (Fig. 3B). This implies that male adults may be more responsive to the environment than female adults.
To further explore the potential biological processes associated with immune response activity, we used gene set enrichment analysis (GSEA). The results of GSEA-KEGG showed that gene sets of neurodegenerative diseases and the Toll and Imd signaling pathways were more significantly upregulated in male adults than in female adults (Fig. 3C), whereas Toll-like receptor (TLR) signaling gene sets were enriched in female adults through GSEA-GO analysis (Fig. 3D). These data indicate that male and female adults perform immune activity in sexually distinct patterns.
Physiological response of the adult head to gut dysbacteriosis
To explore whether the gut microbiota influences gene expression in the S. frugiperda adult head, transcriptomic analysis of female and male adult heads was performed. Gut dysbacteriosis resulted in 20 upregulated and 27 downregulated genes in female adult heads (Table S4), and 40 upregulated and 44 downregulated genes in male adult heads (Fig. S1; Table S5). Additionally, the data revealed that 27 DEGs in total were common in female and male adult heads, while 44 DEGs were unique in differentially expressed gene sets of female adult heads and 150 DEGs were unique in male groups (Fig. 4A). Interestingly, DEGs critical for brain and neural development, including the lactosylceramide 4-alpha-galactosyltransferase gene (A4galt), facilitated trehalose transporter Tret1 (Tret1), 5-methylcytosine rRNA methyltransferase nsun (nsun4), probable galactose-1-phosphate uridylyltransferase (Galt), MICOS complex subunit Mic60 (Mitofilin), solute carrier family 2-facilitated glucose transporter member 1-like (SLC2A3), granzyme-like protein 1 (snk) and gamma-aminobutyric acid receptor subunit delta (GABRB3), were significantly downregulated in treated female adult heads. However, only two genes encoding the major heat shock 70 kDa protein Ba (Hsp68) and sialin (SLC17A5) essential for brain development and immunity, were significantly upregulated in treated female heads. Moreover, the mRNA levels of genes including Tret1, Mitofilin, SLC2A3 and GABRB3 were also inhibited in treated male heads, and there were two unique downregulated genes in male adult heads: octopamine receptor Oamb (Oamb) and sodium- and chloride-dependent GABA transporter 1 (SLC6A1) (Table S4-S5). These data showed that gut dysbacteriosis might also affect the head nervous system of S. frugiperda.
For female heads, KEGG enrichment demonstrated that the most affected pathway was linked to the interferon (IFNα/β) response within the herpes simplex infection pathway (p < 0.05), whereas treated male adult heads showed a closer association with the Toll and Imd signaling pathways (p < 0.05) (Fig. 4B). In addition, GO enrichment analysis indicated that a few DEGs from male heads were involved in biological processes of the humoral immune system, but the DEGs in female samples were concentrated in the inflammatory response of biological processes (p < 0.05) (Fig. 4C). Furthermore, two pathways in the immune system were enriched in control female adult heads according to GSEA-KEGG analysis: complement and coagulation cascades and hematopoietic cell lineage (Fig. 4D-E). Altogether, these results suggest that gut dysbacteriosis has significant and sex-specific impacts on the head transcriptome of S. frugiperda male and female adults, especially in the immune system.
Gut dysbacteriosis diminishes the head immune response
To further identify downregulated immune effectors after gut dysbacteriosis, the expression patterns of DEGs related to immunity in the adult heads, brain and antennae were analyzed. The results showed that the mRNA levels of coding zinc finger antiviral protein (ZAP) genes, including HIVEP2 (ncbi_118280991), GDNF-inducible zinc finger protein 1-like (ncbi_118281138), zinc finger protein 99-like (ncbi_118281047) and ZNF131 (ncbi_118281200) in the interferon pathway, was downregulated in treated female heads (Fig. 5A), while mRNA levels of HIVEP2 and GDNF-inducible zinc finger protein 1-like also showed significant decrease upon dysbacteriosis in male adult heads but zinc finger protein 99-like and ZNF131 had no differences (Table S6; Fig. S2). The mRNA levels of serine/threonine protein kinase Pelle (ncbi_118270757), peptidoglycan-sensing receptor (PGRP, ncbi_118268611), cecropin A (CecA, MSTRG.11940; MSTRG.12130) and cecropin B (cecropin B, MSTRG.11938; MSTRG.12127) in the Toll and Imd signaling pathways were drastically reduced in treated male adult heads (Fig. 5A), whereas all of these genes had no significant differential expression in female adult with dysbacteriosis (Table S7; Fig. S3). The expression of these immune response-related DEGs was confirmed by qRT-PCR analysis (Fig. 5B-B”). In brief, gut dysbacteriosis results in remarkable expression differences in immune genes in the S. frugiperda head in a sex-manner.
Discussion
Antibiotic-induced gut microbiota perturbation
The α-diversity and LEfSe analyses revealed dramatic alterations in the gut microbiota composition and diversity under conditions of antibiotic exposure (Fig. 1C-D). The abundances of the main species within the phylum Firmicutes and the genera Enterococcus and Enterobacter were significantly decreased in treated adult guts (Fig. 1A-B). Firmicutes have been reported to play key roles in nutritional supplementation, energy absorption and host immunity [58]. Moreover, Enterococcus and Enterobacter contribute to the synthesis of vitamins and pheromones, the degradation of plant compounds and nitrogen fixation [59, 60]. Additionally, the PTS system is used for bacteria to absorb carbohydrates and their derivatives through the phosphorylation cascade into the cell [61], and the bacterial PTS was remarkably repressed (Fig. 1E-F). The existence of gut-dominant microorganisms is one of the most important reasons that S. frugiperda are able to feed on various plants and have strong adaptability and resistance. LEfSe analysis revealed that the genera Enterococcus and the family Lactobacillales have lower abundance in the antibiotic-treated groups, which are associated with serotonin, dopamine, acetylcholine, histamine and SCFAs [62]. Therefore, we speculated that the differences in metabolites caused by changes in microbial composition might affect immune cells by accessing the brain or around the brain, and Enterococcus might be a core bacterium that links the gut-brain relationship in S. frugiperda.
Expression difference between untreated female heads and male heads
A comparison between male and female adult heads revealed only 50 differentially expressed genes (Fig. 3A). In contrast, we found that the expression of four key coding genes of OBPs was greater in female adult heads (Fig. 3A-B). Insect chemosensory systems play important roles in reproductive success and survival, while OBPs are essential for detecting and distinguishing specific odors by transporting odorant molecules to olfactory receptors [63]. In the natural environment, the olfaction, gustation, vision and tactile sensing systems of female adult insects are simultaneously used to recognize suitable egg-laying locations and reduce competitive pressure [64]. Furthermore, DEGs of the RNA cap binding complex and the Toll-like receptor signaling pathways associated with inflammatory factors all showed higher expression in female heads (Fig. 3B and D). TLRs play significant roles in the host innate immune system, which can recognize various classes of pathogens, including bacteria, fungi, viruses and protozoa [65], and TLRs are necessary for immune cell survival, regulation and proliferation [66]. These data suggest that female S. frugiperda adults have a stable immune system, which helps female adults to better resist external environmental invasion for reproduction and egg laying. However, male adults had higher expression within the Toll and Imd signaling pathways and a different humoral immune response to combat microbial infection (Fig. 3C). Male and female innate immune responses are different. For example, the neuronal diseases Parkinson’s disease [67] and Alzheimer’s disease [68] show sex differences. Additionally, gender differences in brain immune function have also been reported across diverse species, with females generally having more dynamic immune function than males [69]. The Tolllike receptor 7 (TLR7) gene is located on the X chromosome and can escape X inactivation. Therefore, the expression level of TLR7 in females is higher than that in males [70]. Interferonα production in cells from women is higher than that in cells from men when peripheral blood mononuclear cells (PBMCs) are exposed to TLR7 ligands in vitro [71]. Our results also indicate that the Tolllike receptor signaling pathway shows higher expression in the female adult heads (Fig. 3D), which is consistent with reports in mammals. The reasons for this biological phenomenon might be genetic mediators, age and reproductive status, hormonal mediators and environmental mediators [72].
Female and male heads respond differently to gut dysbacteriosis
RNA-seq analysis of adult heads showed that gut dysbacteriosis altered the expression of more transcripts in males than that in females (Fig. 5A). Gut dysbacteriosis induced the transcriptional downregulation of antimicrobial peptides (AMPs) in the treated male S. frugiperda adult heads, and the repression of ZAP involved in the interferon IFNα/β pathway in the treated female S. frugiperda adult heads (Fig. 5). ZAP, as a kind of restriction factor, mediates the degradation of viral mRNA, inhibits viral translation and ultimately suppresses viral replication [73]. Thus, the gut microbiome likely contributes to the immune system in the head from S. frugiperda adults in a sex-specific manner. The gut microbiota affects the proportions, migration and functions of various immune cells [74]. Numerous examples have illustrated that gut bacteria can control both innate and adaptive immune responses on the mucosal surface during infection and inflammation [75]. However, how gut microbes control microglial maturation and innate immune reactions in the insect head remains incompletely understood. Importantly, both adult male and female heads of gut dysbacteriosis showed marked decreases in the expression of genes linked to brain and neural development compared with the levels in conventionally raised S. frugiperda (Table S4-S5). Neurotransmitters that transmit information between neurons are essential for brain functions and behaviors. Recent studies have indicated that microorganisms regulate the stimulation of glial cells by producing microbial metabolites. For example, SCFAs, neurotransmitters and gut hormones [76, 77] enter the circulatory system and signal to the host brain in this way [78]. Some studies have also found a cross-regulatory mechanism between the gut nervous system and the immune system. This once again confirms that gut microbiota can affect the brain through the immune system [79]. For example, compared to conventional flies, Drosophila are hyperactive under axenic conditions [80]. Additionally, the results of evaluating gene expression differences through qRT‒PCR indicated that DEGs related to the immune pathway changed more significantly in the brain than in other tissues (Fig. 5B-B”). According to reports, the toll channel in the central nervous system promotes immune responses by providing key factors, but excessive immune responses can cause damage to neurons [81]. Our study showed that adult brains of S. frugiperda are particularly affected when complex microorganisms are lacking. Sex-specific differences reveal that male adults are more sensitive to changes in microbial communities than females. We speculate that the immune signaling proteins of female adults are of the stable and constitutive type, while the immune proteins of male adults are of the susceptible and inducible type. Furthermore, the differences in microbial composition between female adults and male adults resulted in different regulatory functions of metabolic metabolites. More research is needed to confirm these findings and better understand sex, microbiota manipulation, microbiota-linked metabolites, and immune responses. Therefore, we will further explain the relationship between gut bacteria and the brain by supplementing the microbial community and conducting a comprehensive analysis. We will also detect behavioral changes in the treated S. frugiperda.
Conclusion
In summary, we demonstrate the head transcriptomic differences between female and male adults of S. frugiperda. In particular, we found the potential impact of gut bacteria on the head immune pathway. Gut microbiota dysbiosis may drive inhibition of the IFN pathway in female adult heads, but Toll and Imd signaling pathway in male adult heads, thereby decreasing the expression level of genes encoding antibacterial protein. Gut microbiota dysbiosis also resulted in the gene expression changes involved in the neurodevelopment. Furthermore, S. frugiperda adults showed sex-specific response differences to gut dysbiosis. Further research involved in the underlying mechanism by how gut bacteria affect head immune pathways and brain neurodevelopment need to be substantiated. In brief, our study suggests a probable linkage between gut microbiota dysbiosis and the brain immunity in S. frugiperda. These findings benefit our understanding of the impact of gut symbionts on the ecology and evolution of insects and provide some model references for further research on psychiatric disorders.
Data Availability
All data needed to evaluate the conclusions in this paper are presented in the paper and the Supplementary Materials. All data of the raw reads have been deposited into the NCBI Sequence Read Archive (SRA) database (accession numbers PRJNA986514, PRJNA986662 and PRJNA983288). The datasets generated and/or analyzed during the current study are available from the corresponding author on reasonable request.
References
Castellani G, Croese T, Peralta Ramos JM, Schwartz M. Transforming the understanding of brain immunity. Science. 2023;380:eabo7649.
Wrona D. Neural-immune interactions: an integrative view of the bidirectional relationship between the brain and immune systems. J Neuroimmunol. 2006;172(1–2):38–58.
Dantzer R. Cytokine-induced sickness behavior: a neuroimmune response to activation of innate immunity. Eur J Pharmacol. 2004;500(1–3):399–411.
Udit S, Blake K, Chiu IM. Somatosensory and autonomic neuronal regulation of the immune response. Nat Rev Neurosci. 2022;23(3):157–71.
Ziemssen T, Kern S. Psychoneuroimmunology-cross-talk between the immune and nervous systems. J Neurol. 2007;254:8–11.
Liberti J, Engel P. The gut microbiota-brain axis of insects. Insect Sci. 2020;39:6–13.
Adamo SA. Comparative psychoneuroimmunology: evidence from the insects. Behav Cogn Neurosci Rev. 2006;5:128–40.
Zhang Z, Mu X, Cao Q, Shi Y, Hu X, Zheng H. Honeybee gut Lactobacillus modulates host learning and memory behaviors by regulating tryptophan metabolism. Nat Commun. 2022;13(1):2037–49.
Seo DO, O’Donnell D, Jain N, Ulrich JD, Herz J, Li Y, Holtzman DM. ApoE isoform-and microbiota-dependent progression of neurodegeneration in a mouse model of tauopathy. Science. 2023;379:eadd1236.
Johnson KVA, Foster KR. Why does the microbiome affect behavior? Nat Rev Microbiol. 2018;16(10):647–55.
Quigley EM. Microbiota-brain-gut axis and neurodegenerative Diseases. Curr Neurol Neurosci Rep. 2017;17:1–9.
Davidson GL, Cooke AC, Johnson CN, Quinn JL. The gut microbiome as a driver of individual variation in cognition and functional behavior. Philosophical Trans Royal Soc B. 2018;373:1–12.
Shen L, Ji H-F. Associations between gut microbiota and Alzheimer’s Disease: current evidence and future therapeutic and diagnostic perspectives. J Alzheimers Dis. 2019;68(1):25–31.
Sampson TR, Debelius JW, Thron T, Janssen S, Shastri GG, Ilhan ZE, et al. Gut microbiota regulate motor deficits and neuroinflammation in a model of Parkinson’s Disease. Cell. 2016;167(6):1469–80.
Vuong HE, Yano JM, Fung TC, Hsiao EY. The microbiome and host behavior. Annu Rev Neurosci. 2017;40:21–49.
Sgritta M, Dooling SW, Buffington SA, Momin EN, Francis MB, Britton RA, Costa-Mattioli M. Mechanisms underlying microbial-mediated changes in social behavior in mouse models of autism spectrum disorder. Neuron. 2019;101:246–59.
Zhu F, Guo R, Wang W, Ju Y, Wang Q, Ma Q, et al. Transplantation of microbiota from drug-free patients with schizophrenia causes schizophrenia-like abnormal behaviors and dysregulated kynurenine metabolism in mice. Mol Psychiatry. 2020;25(11):2905–18.
Sherwin E, Bordenstein SR, Quinn JL, Dinan TG, Cryan JF. Microbiota and the social brain. Science. 2019;366:1–17.
Philipp E, Nancy AM. The gut microbiota of insects-diversity in structure and function. FEMS Microbiol Rev. 2013;37(5):699–735.
Carthey AJR, Gillings MR, Blumstein DT. The extended genotype: microbially mediated olfactory communication. Trends Ecol Evol. 2018;33(11):885–94.
Heijtz RD, Wang S, Anuar F, Qian Y, Björkholm B, Samuelsson A, et al. Normal gut microbiota modulates brain development and behavior. Proc Natl Acad Sci. 2011;108(7):3047–52.
Matsuura K. Nestmate recognition mediated by intestinal bacteria in a termite, Reticulitermes Speratus. Oikos. 2011;92(1):20–6.
Zhang Z, Mu X, Cao Q, Shi Y, Hu X, Zheng H. Honeybee gut microbiota modulates host behaviors and neurological processes. Cold Spring Harbor Laboratory; 2020.
Denieu M, Mounts K, Bhojwani D, Manier MK. Two gut microbes are necessary and sufficient for normal cognition in drosophila melanogaster. BioRxiv. 2019;593723:1–18.
Johnson SJ. Migration and the life history strategy of the fall armyworm, Spodoptera frugiperda in the Western Hemisphere. Int J Trop Insect Sci. 1987;8(4–5–6):543–9.
Early R, González-Moreno P, Murphy ST, Day R. Forecasting the global extent of invasion of the cereal pest Spodoptera frugiperda, the fall armyworm. BioRxiv, 2018; 391847.
Lv D, Liu X, Dong Y, Yan Z, Zhang X, Wang P, et al. Comparison of gut bacterial communities of fall armyworm (Spodoptera frugiperda) reared on different host plants. Int J Mol Sci. 2021;22(20):1–16.
Chen Y-P, Li Y-H, Sun Z-X, Du E-W, Lu Z-H, Li H, Gui FR. Effects of host plants on bacterial community structure in larvae midgut of Spodoptera frugiperda. Insects. 2022;13(4):1–13.
Higuita Palacio MF, Montoya OI, Saldamando CI, García-Bonilla E, Junca H, Cadavid-Restrepo GE, Moreno-Herrera CX. Dry and rainy seasons significantly alter the gut microbiome composition and reveal a key Enterococcus sp. (Lactobacillales: Enterococcaceae) core component in Spodoptera frugiperda (Lepidoptera: Noctuidae) corn strain from Northwestern Colombia. J Insect Sci. 2021;21(6):1–11.
Li D-D, Li JY, Hu Z-Q, Liu T-X, Zhang S. Fall armyworm gut bacterial diversity associated with different developmental stages, environmental habitats, and diets. Insects. 2022;13(9):762–77.
Chen J, Ma Y, Huang S, Li J, Zhang Y, Wang H, Pang R. The dynamics of the microbial community in fall armyworm Spodoptera frugiperda during a life cycle. Entomol Exp Appl. 2023;0:1–12.
Chen Y, Zhou H, Lai Y, Chen Q, Yu X-Q, Wang X. Gut microbiota dysbiosis influences metabolic homeostasis in Spodoptera frugiperda. Front Microbiol. 2021;12:1–13.
Lü D, Dong Y, Yan Z, Liu X, Zhang Y, Yang D, et al. Dynamics of gut microflora across the life cycle of Spodoptera frugiperda and its effects on the feeding and growth of larvae. Pest Manag Sci. 2023;79(1):173–82.
Gomes AFF, Omoto C, Cônsoli FL. Gut bacteria of field-collected larvae of Spodoptera frugiperda undergo selection and are more diverse and active in metabolizing multiple insecticides than laboratory-selected resistant strains. J Pest Sci. 2020;93(2):833–51.
Jia B, Liu Y, Zhu Y-C, Liu X, Gao C, Shen J. Inheritance, fitness cost and mechanism of resistance to tebufenozide in Spodoptera exigua (Hübner) (Lepidoptera: Noctuidae). Pest Manage Science: Former Pesticide Sci. 2009;65(9):996–1002.
Elgart M, Stern S, Salton O, Gnainsky Y, Heifetz Y, Soen Y. Impact of gut microbiota on the fly’s germ line. Nat Commun. 2016;7(1):11280–91.
GuoM,WuF,HaoG,QiQ,LiR,LiN,etal.Bacillus subtilis improves immunity and disease resistance in rabbits.Frontiers in immunology,2017,8:354.
Chen S, Zhou Y, Chen Y, Gu J. Fastp: an ultrafast all-in-one FASTQ preprocessor. Bioinformatics. 2018;34:884–90.
Edgar RC. Uparse: highly accurate otu sequences from microbial amplicon reads. Nat Methods. 2013;10(10):996–8.
Wang Q, Garrity GM, Tiedje JM, Cole JR. Naive bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007;73:5261–7.
Pruesse E, Quast C, Knittel K, Fuchs BM, Ludwig W, Peplies J, Glöckner. F O. SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucleic Acids Res. 2007;35:7188–96.
Nilsson RH, Larsson KH, Taylor AFS, Bengtsson-Palme J, Jeppesen TS, Schigel D, et al. The UNITE database for molecular identification of fungi: handling dark taxa and parallel taxonomic classifications. Nucleic Acids Res. 2019;47(1):259–64.
Wickham H. ggplot2. Wiley Interdisciplinary Reviews. Comput Stat. 2011;3(2):180–5.
Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Knight R. Caparose Jgkj, stombaugh j, bittinger k, bushman fd. Qiime allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.
Oksanen J, Kindt R, Legendre P, O’Hara B, Oksanen M. VEGAN: community ecology package.2010.
Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, Huttenhower C. Metagenomic biomarker discovery and explanation. Genome Biol. 2011;12(6):R60.
Aßhauer KP, Wemheuer B, Daniel R, Meinicke P. Tax4Fun: predicting functional profiles from metagenomic 16S rRNA data. Bioinformatics. 2015;31:2882–4.
Oksanen J. Vegan: community ecology package. R package version 1.17-4. http://CRAN.R-project.org/package=vegan, 2011.
Chen S, Zhou Y, Chen Y, Gu J. Fastp: an ultra-fast all-in-one FASTQ preprocessor. BioRxiv. 2018; 274100.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9(4):357–9.
Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12(4):357.
Pertea M, Pertea GM, Antonescu CM, Chang TC, Mendell JT, Salzberg SL. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat Biotechnol. 2015;33(3):290.
Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12:323.
Love MI, Huber W, Anders S. Moderated estimation of Fold change and dispersion for RNAseq data with DESeq2. Genome Biol. 2014;15(12):550.
Love MI, Huber W, Anders S. Moderated estimation of Fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:1–21.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci. 2005;102(43):15545–50.
Colston TJ, Jackson CR. Microbiome evolution along divergent branches of the vertebrate tree of life: what is known and unknown. Mol Ecol. 2016;25(16):3776–800.
Lilburn TG, Kim KS, Ostrom NE, Byzek KR, Leadbetter JR, Breznak JA. Nitrogen fixation by symbiotic and free-living spirochetes. Science. 2001;292(5526):2495–8.
Morales-Jiménez J, Zúñiga G, Ramírez-Saad HC, Hernández-Rodríguez C. Gut-associated bacteria throughout the life cycle of the bark beetle Dendroctonus rhizophagus Thomas and Bright (Curculionidae: Scolytinae) and their cellulolytic activities. Microb Ecol. 2012;64:268–78.
Tanaka Y, Teramoto H, Inui M, Yukawa H. Regulation of expression of general components of the phosphoenolpyruvate: carbohydrate phosphotransferase system (PTS) by the global regulator SugR in Corynebacterium glutamicum. Appl Microbiol Biotechnol. 2008;78:309–18.
Benakis C, Martin-Gallausiaux C, Trezzi JP, Melton P, Liesz A, Wilmes P. The microbiome-gut-brain axis in acute and chronic brain Diseases. Curr Opin Neurobiol. 2020;61:1–9.
Gong D-P, Zhang H-J, Zhao P, Xia Q-Y, Xiang Z-H. The odorant binding protein gene family from the genome of silkworm, Bombyx mori. BMC Genomics. 2009;10(1):1–14.
Zhang L, Yu J, Guo X, Wei J, Liu T, Zhang W. Parallel mechanosensory pathways direct oviposition decision-making in Drosophila. Curr Biol. 2020;30(16):3075–88.
Kawai T, Akira S. The role of pattern-recognition receptors in innate immunity: update on toll-like receptors. Nat Immunol. 2010;11(5):373–84.
Li X, Jiang S, Tapping RI. Toll-like receptor signaling in cell proliferation and survival. Cytokine. 2010;49(1):1–9.
Bakeberg MC, Gorecki AM, Kenna JE, Jefferson A, Byrnes M, Ghosh S, et al. Differential effects of sex on longitudinal patterns of cognitive decline in Parkinson’s Disease. J Neurol. 2021;268:1903–12.
Ferretti MT, Iulita MF, Cavedo E, Chiesa PA, Schumacher Dimech A, Santuccione Chadha A, et al. Sex differences in Alzheimer disease-the gateway to precision medicine. Nat Reviews Neurol. 2018;14(8):457–69.
Nunn CL, Lindenfors P, Pursall ER, Rolff J. On sexual dimorphism in immune function. Philosophical Trans Royal Soc B: Biol Sci. 2009;364(1513):61–9.
Pisitkun P, Deane JA, Difilippantonio MJ, Tarasenko T, Satterthwaite AB, Bolland S. Autoreactive B-cell responses to RNA-related antigens due to TLR7 gene duplication. Science. 2006;312(5780):1669–72.
Berghöfer B, Frommer T, Haley G, Fink L, Bein G, Hackstein H. Tlr7 ligands induce higher ifn-alpha production in females. J Immunol. 2006;177(4):2088–96.
Klein SL, Flanagan KL. Sex differences in immune responses. Nat Rev Immunol. 2016;16(10):626–38.
Gao G, Guo X, Goff SP. Inhibition of retroviral RNA production by ZAP, a CCCH-type zinc finger protein. Science. 2002;297(5587):1703–6.
Rooks MG, Garrett WS. Gut microbiota, metabolites and host immunity. Nat Rev Immunol. 2016;16(6):341–52.
Kamada N, Seo SU, Chen GY, Núñez G. Role of the gut microbiota in immunity and inflammatory Disease. Nat Rev Immunol. 2013;13(5):321–35.
Erny D, Hrabě de Angelis AL, Jaitin D, Wieghofer P, Staszewski O, David E, et al. Host microbiota constantly control maturation and function of microglia in the CNS. Nat Neurosci. 2015;18(7):965–77.
Spichak S, Donoso F, Moloney GM, Gunnigle E, Brown JM, Codagnone M et al. Microbially derived short-chain fatty acids impact astrocyte gene expression in a sex-specific manner. Brain, behavior, & immunity-health. 2021; 16:100318.
Sarkar A, Lehto SM, Harty S, Dinan TG, Cryan JF, Burnet PW. Psychobiotics and the manipulation of bacteria-gut-brain signals. Trends Neurosci. 2016;39(11):763–81.
Muller PA, Koscsó B, Rajani GM, Stevanovic K, Berres ML, Hashimoto D, et al. Crosstalk between muscularis macrophages and enteric neurons regulates gastrointestinal motility. Cell. 2014;158(2):300–13.
Schretter CE, Vielmetter J, Bartos I, Marka Z, Marka S, Argade S, Mazmanian SK. A gut microbial factor modulates locomotor behavior in Drosophila. Nature. 2018;563(7731):402–6.
Mao D, Lai T, Chen Y, Song B, Song Y. Roles of toll signaling pathway in response to injuries of the Central Nervous System. Chin J Biochem Mol Biology. 2010;26(2):121–7.
Acknowledgements
Not applicable.
Funding
This study was supported by a grant from the Laboratory of Lingnan Modern Agriculture Project (Grant number NZ2021020), Natural Science Foundation of Guangdong Province-Surface project (Grant number 2023A1515012458), Forestry Science and Technology Innovation Project of Shaanxi Province (SXLK2021-0101) and Young Teachers Research and Cultivation Foundation of South China Normal University (Grant number22KJ28).
Author information
Authors and Affiliations
Contributions
Author Contributions: H.D. and Q.F. designed and monitored the project. J. F., J.W., and X. H. performed the experiments. J.F. interpreted the data. X.H., R. Z. and X.L. prepared some of the experimental materials. J.F. wrote the manuscript. H.D. and Q.F. finalized the manuscript.All authors reviewed the manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Junrui-Fu, Rong, Z., Huang, X. et al. Gut dysbacteriosis induces expression differences in the adult head transcriptome of Spodoptera frugiperda in a sex-specific manner. BMC Microbiol 23, 388 (2023). https://doi.org/10.1186/s12866-023-03089-0
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12866-023-03089-0