Drosophila larvae shape topographies and carrying capacities of bacterial community.

(A-D) Representative images of sticky “biofilm-like” formation on the surface of the sugar-corn-yeast medium whereby Drosophila flies with differential robustness were raised. The topographies of surface slick are differentially deconstructed and segmented by flies with different robust flies. (E) Bacterial loads of the diet associated with strong, crowed, weak and infertile flies, respectively. (F) Bacterial loads of the diet associated with male, virgin and aged flies, respectively. (G) Bacterial loads of the diet associated with Drosophila larvae in a dosage-dependent manner. n = 6 for each. Error bars indicate SEM. All variables have different letters, and they are significantly different (p < 0.01). Kruskal–Wallis test followed by Dunn’s multiple comparisons test.

Drosophila larvae outcompete S. marcescens in the diet.

(A) A diagram of a reductionist approach to investigate the role of Drosophila in regulating the physiology and behavior of S. marcescens. Top: Germ-free Drosophila larvae were generated by successive sterilization of fresh eggs with sanitizer Walch, sodium hypochloride (SH), ethanol, and PBS containing 0.01% TritonX-100T (PBST). Bottom: S. marcescens was cultured in a liquid medium and re-inoculated to fly cornmeal food after washing with PBS buffer. In the meantime, GF crawling larvae were transferred to the fly medium in the shared vials with S. marcescens. (B) Representative images of surface slick inoculated with S. marcescens alone and with S. marcescens over time. (C) The prodigiosin production of S. marcescens alone and in coculture at different time points. Prodigiosin production was assessed with the spectrometer at OD534. (D) The bacterial load of S. marcescens alone and in coculture in the time course. (E) The survival rate of adult flies challenged with S. marcescens alone and in coculture. Single and coculturing S. marcescens were obtained after 24-h incubation as described in Figure 1A, and the percentage of living female flies treated with S. marcescens alone and in coculture was calculated to monitor lifespan. n = 180 for each. (F) RT-qPCR analysis of the expression levels of virulence-related genes of with S. marcescens alone and in coculture. n = 3 for each. (G) Transmission electron microscopy of S. marcescens alone and in coculture. Scale bars: 400 nm (left panel) or 200 nm (right panel). (H) RT-qPCR analysis of the expression levels of extracellular polysaccharide production-related genes in the control and larvae groups. n = 3 for each. Error bars indicate SEM. All variables have different letters, and they are significantly different (p < 0.01). Kruskal– Wallis test followed by Dunn’s multiple comparisons test.

Drosophila larvae adjust bacterial global transcriptional adaptation to the host.

(A) Principal coordinate analysis (PCA) of unweighted, jack-knifed UniFrac distances of the transcriptional profile of S. marcescens alone and with larvae. PC1, principal coordinate 1; PC2, principal coordinate 2. Scattered dots in different colors represent samples from different experimental groups. n = 4-5. (B) Volcano plot comparing gene expression profiles of S. marcescens alone and with larvae after 24 h of incubation. X-axis represents the log2-transformed value of gene expression change folds between larvae and control groups. Y-axis represents the logarithmic transformation value of gene expression levels in S. marcescens. Genes belonging to different pathways are represented by different colored shapes as indicated. ▽ depicts genes significantly upregulated in S. marcescens with larvae compared to S. marcescens alone (log2 fold change < 1; adjusted p < 0.01), and △ depicts genes significantly downregulated in S. marcescens with larvae (log2 fold change < 1; adjusted p < 0.01) compared to S. marcescens alone. ○ depicts genes without significant alteration compared to S. marcescens alone. (C, D) Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of the significantly downregulated and upregulated genes in S. marcescens with larvae compared to S. marcescens alone. (E, F) RT-qPCR analysis of the expression levels of downregulated and upregulated genes in the control and larvae groups. n = 3 for each. Error bars indicate SEM.

Drosophila larvae affect the global metabolic profile of S. marcescens.

(A) Principal coordinate analysis (PCA) of unweighted, jack-knifed UniFrac distances of metabolic profile of S. marcescens alone and with larvae. PC1, principal coordinate 1; PC2, principal coordinate 2. Scattered dots in different colors represent samples from different experimental groups. n = 5 for each. (B) Volcano plot comparing metabolic profiles of between control and larvae groups after 24 h of incubation. X-axis represents the log2-transformed value of gene expression change folds between larvae and control groups. Y-axis represents the logarithmic transformation value of gene expression levels in S. marcescens. Red dots depict genes significantly upregulated in S. marcescens with larvae compared to S. marcescens alone (log2 fold change < 1; adjusted p < 0.01), and blue dots depict genes significantly downregulated in S. marcescens with larvae (log2 fold change < 1; adjusted p < 0.01) compared to S. marcescens alone. Grey dots depict genes without significant alteration compared to S. marcescens alone. (C) The distinct clusters of metabolites in S. marcescens alone versus co-culture. (D, E) KEGG pathway analysis of the significantly downregulated and upregulated metabolites in S. marcescens with larvae compared to S. marcescens alone.

Pathogenicity heterogeneity of S. marcescens.

(A) mRNA gene counts per cell for S. marcescens alone, with force and with larvae. Each dot represents a bacterial cell of S. marcescens. (B) Joint UMAP two-dimensional analysis showing that are distinct clusters among S. marcescens alone, with force and with larvae. (C) The cell subpopulation among the control, force and larvae groups. There were three distinct subpopulations in the control and force groups. (D) Mean expression levels of genes involved in ABC transporter, Quorum sensing, Secretion system, Two-component system, LPS and Peptidoglycan biosynthesis and Virulence-related genes in different subclusters. The shape of each dot indicates the proportion of cells in the cluster, while the color indicates the average activity normalized from 0 % to 100 % across all clusters. (E, F) The expression of a representative gene of ABC transporter and Quorum sensing was highlighted on the UMAP. The red color bars represent the normalized expression of a gene across all cells analyzed. (G, H) Violin plots of livI and oppA gene in different subclusters. Each dot represents a single cell and the shapes represent the expression distribution. (I, J) The expression of a representative gene of Secretion system and Virulence-related genes was highlighted on the UMAP. The red color bars represent the normalized expression of a gene across all cells analyzed. (K, L) Violin plots of secY and fp gene in different subclusters. Each dot represents a single cell and the shapes represent the expression distribution.

Growth heterogeneity of S. marcescens.

(A) Mean expression levels of genes involved in Ribsome, DNA replication, Nitrogen metabolism and Carbon metabolism in different subclusters. The shape of each dot indicates the proportion of cells in the cluster, while the color indicates the average activity normalized from 0 % to 100 % across all clusters. (B, C) The expression of a representative gene of Ribsome and DNA replication was highlighted on the UMAP. The red color bars represent the normalized expression of a gene across all cells analyzed. (D, E) Violin plots of rpsL and trpD genes in different subclusters. Each dot represents a single cell and the shapes represent the expression distribution. (F, G) The expression of two representative genes of Nitrogen metabolism and Carbon metabolism was highlighted on the UMAP. The red color bars represent the normalized expression of a gene across all cells analyzed. (H, I) Violin plots of glnK, glnA, pgk and adhE genes in different subclusters. Each dot represents a single cell and the shapes represent the expression distribution. (J) Schematic of the pathogenicity and commensalism regulatory pathway.

Larvae-derived AMPs antagonize S. marcescens.

(A) Representative images of surface slick with S. marcescens alone, with larvae, with secreta and with AMPs. (B) The prodigiosin production of S. marcescens alone, with larvae, with secreta and with AMPs. (C) Bacterial loads of S. marcescens alone, with larvae, with secreta and with AMPs. (D) Representative images of surface slick with S. marcescens alone, with wild-type larvae, with ΔAMP14 larvae, with AMPs and with ΔAMP14 larvae+AMPs. (E) The prodigiosin production of S. marcescens alone, with wild-type larvae, with ΔAMP14 larvae, with AMPs and with ΔAMP14 larvae+AMPs. (F) Bacterial loads of with S. marcescens alone, with wild-type larvae, with ΔAMP14 larvae, with AMPs and with ΔAMP14 larvae+AMPs. (G, H) RT-qPCR analysis of the expression levels of downregulated and upregulated genes in the S. marcescens alone, with wild-type larvae, with ΔAMP14larvae, with AMPs and with ΔAMP14 larvae+AMPs. For B-C and E-F, n = 6 for each. Fo G-H, n = 3 for each. Error bars indicate SEM. All variables have different letters, and they are significantly different (p < 0.01). Kruskal–Wallis test followed by Dunn’s multiple comparisons test.

Drosophila larvae modulate S. marcescens lifestyle switch.

(A) Representative images of surface slick inoculated with S. marcescens in coculture with different numbers of crawling larvae. (B) The prodigiosin production of S. marcescens in coculture with a series of crawling larvae. n = 6 for each. (C) The bacterial load of S. marcescens in coculture with different numbers of crawling larvae. n = 6 for each. (D) S. marcescens promoted the development of GF larvae. Developmental timing of GF, L. plantarum-, and S. marcescens-associated Drosophila was assessed on the poor diet with 0.5% yeast. The cumulative percentage of the pupation emergence is shown over time. n = 60 for each. Error bars indicate SEM. All variables have different letters, and they are significantly different (p < 0.01). Kruskal–Wallis test followed by Dunn’s multiple comparisons test.

Mechanical force didn’t contribute to S. marcescens lifestyle.

(A) The schematic illustration of S. marcescens alone, with force and with larvae. Single and coculturing S. marcescens were generated as described in Figure 1A. In the meantime, S. marcescens with force were agitated using sterile glass sticks at the 2-h interval. (B) Surface slicks associated with S. marcescens alone, with force and with larvae. (C) The prodigiosin production of S. marcescens alone, with force and with larvae. n = 6 for each. (D) Bacterial loads of S. marcescens alone, with force and with larvae. n = 6 for each. (E) RT-qPCR analysis of the expression levels of the pigA, pigC, pigM and pigI genes of S. marcescens alone, with force and with larvae. n = 3 for each. Error bars indicate SEM. All variables have different letters, they are significantly different (p < 0.01). Kruskal–Wallis test followed by Dunn’s multiple comparisons test.

Interaction network analysis of transcriptome and metabolome.

(A) Interaction network analysis of transcriptome and metabolome of S. marcescens alone and with larvae. 42 differentially expressed metabolites (red circles) are related to differentially expressedgenes involved in the ribosome (light brown), transcription (dark brown), DNA replication (light purple), energy metabolism (dark green), ABC transporters (dark blue), phosphotransferase system (light blue), quorum sensing (dark purple) and exosome (light blue). (B) The prodigiosin production of S. marcescens alone and with predicted metabolites. (C) Bacterial loads of S. marcescens alone and with predicted metabolites. n = 6 for each. (D) The survival rate of adult flies challenged with S. marcescens alone and with metabolites. (E) RT-qPCR analysis of the expression levels of virulence-related genes of with S. marcescens alone and with metabolites. n = 3 for each. Error bars indicate SEM. All variables have different letters, and they are significantly different (p < 0.01). Kruskal–Wallis test followed by Dunn’s multiple comparisons test.

Heat-induced phenotypic heterogeneity of S. marcescens.

(A) mRNA gene counts per cell for the control and heat-shock groups. Each dot represents a bacterial cell of S. marcescens. (B) Joint UMAP two-dimensional analysis of S. marcescens showing that are distinct clusters between the control and heat-shock groups. (C) The cell subpopulation typing of the control and heat-shock groups. (E, F) The expression of two representative genes of heat shock and stress response was highlighted on the UMAP. The red color bars represent the normalized expression of a gene across all cells analyzed.

(A) The survival rate of adult flies challenged with S. marcescens alone, with agitation and in coculture. S. marcescens with alone, with agitation and in coculture were obtained after 24-h incubation as described in Figure 1A, and the percentage of living female flies was calculated to monitor lifespan. n = 180 for each. (B, C) The expression of two representative genes of Nitrogen metabolism and Carbon metabolism was highlighted on the UMAP. The red color bars represent the normalized expression of a gene across all cells analyzed. (D, E) Violin plots of glnK, glnA, pgk and adhE genes in different subclusters. Each dot represents a single cell and the shapes represent the expression distribution.