Quorum Sensing Signals Alter in vitro Soil Virus Abundance and Bacterial Community Composition

Cell-density dependent quorum sensing (QS) is fundamental for many coordinated behaviors among bacteria. Most recently several studies have revealed a role for bacterial QS communication in bacteriophage (phage) reproductive decisions. However, QS based phage-host interactions remain largely unknown, with the mechanistic details revealed for only a few phage-host pairs and a dearth of information available at the microbial community level. Here we report on the specific action of eight different individual QS signals (acyl-homoserine lactones; AHLs varying in acyl-chain length from four to 14 carbon atoms) on prophage induction in soil microbial communities. We show QS autoinducers, triggered prophage induction in soil bacteria and the response was significant enough to alter bacterial community composition in vitro. AHL treatment significantly decreased the bacterial diversity (Shannon Index) but did not significantly impact species richness. Exposure to short chain-length AHLs resulted in a decrease in the abundance of different taxa than exposure to higher molecular weight AHLs. Each AHL targeted a different subset of bacterial taxa. Our observations indicate that individual AHLs may trigger prophage induction in different bacterial taxa leading to changes in microbial community structure. The findings also have implications for the role of phage-host interactions in ecologically significant processes such as biogeochemical cycles, and phage mediated transfer of host genes, e.g., photosynthesis and heavy metal/antibiotic resistance.


INTRODUCTION
Bacteriophages may infect bacterial host cells via lytic and lysogenic cycles, both of which have shown ecological significance. For instance, lytic cycles of reproduction can impact population and community dynamics through lysis of host cells effectively re-routing dissolved organic carbon and other nutrients back to the dissolved pool, a process referred to as the "viral shunt" (Danovaro et al., 2008;Sullivan et al., 2017). Lysogenic cycles in which the phage genome is inserted into the host cell genome without killing the host, may promote host fitness and regulate metabolic functions through selective expression of certain phage encoded genes and transcriptional regulators without production of progeny phage particles (Hurwitz and U'Ren, 2016;Howard-Varona et al., 2017; the QS-mediated prophage-induction mechanisms are likely operative in natural systems. However, the potential significance of this phenomena at the microbial community level has only been demonstrated by Ghosh et al. (2009).
Communication among phages through phage-encoded arbitrium peptides was first described by Erez et al. (2017), and the following studies (Dou et al., 2018;Wang et al., 2018;Gallego del Sol et al., 2019;Stokar-Avihail et al., 2019) revealed the molecular basis for the production, detection, and consequences of the short signaling peptides on phages lysogenylysis decisions. Notably, these reports also showed that phages communicate only with their close relatives using a very specific arbitrium peptide, which suggests that phage communication peptides act in a taxon-specific manner just as bacterial QS signals. Inspired by the above studies, especially the phagebacterium QS connections, we hypothesized that any single QS signal should only induce prophages within a small subset of closely related host bacteria. Toward that end, we tried to determine the impacts of the addition of individual AHL signaling molecules on prophage-induction and further assess the resulting impact of phage-mediated host cell lysis on bacterial community composition in vitro using microbial communities extracted directly from soil.

Sample Collection and Bacterial Extraction
Soil samples were collected from an agricultural field at the East Tennessee Agricultural Research and Education Center (Latitude = 35.899166; Longitude = -83.961120) on March 5, 2019. Soil samples were collected from one location and thoroughly mixed onsite. The extraction of bacterial cells from soil was performed as described elsewhere (Williamson et al., 2005). For sufficient microbial biomass to complete all the induction assays described below, a larger quantity of soil (300 g) was extracted using cold (stored at 4 • C) extraction buffer consisting of 10 g/L potassium citrate, 1.44 g/L Na 2 HPO 4 · 7H 2 O, and 0.24 g/L KH 2 PO 4 . The soil to extraction buffer ratio was maintained as described in Williamson et al. (2005), and the mixture was blended in a sterilized blender vessel at the maximum speed for 3 min. The extracted bacteria were concentrated by centrifuging the slurries on a cushion of 60% (w/v) Nycodenz solution (Accurate Chemical & Scientific, Westbury, NY, United States) at 4,000 g for 20 min at 4 • C. Supernatant including the Nycodenz phase that contains the bacterial cells was transferred to 50 ml centrifuge tubes and was then centrifuged at 5,000 × g for 20 min. The bacterial pellets were washed twice using sterile extraction buffer and resuspended in autoclaved 0.2 µmfiltered soil extract (Ghosh et al., 2009) after decanting the supernatant to discard extracellular viruses. The resuspended bacterial extracts were pooled and homogenized by gently shaking. The bacteria and viruses were enumerated after incubation via epifluorescence microscopy direct counting as described below.

Epifluorescence Microscopy Counting
Epifluorescence microscopy was used for enumeration of bacteria and viruses in the solution (Williamson et al., 2005;Liang et al., 2019b). For removal of extracellular DNA, Deoxyribonuclease I (DNase I, 2.5 units/µl, Thermo Scientific) was used to treat all suspensions. After DNase I treatment, each bacterial solution was filtered through a set of filters, a 0.2-µm-pore-size Isopore membrane filter (Merck Millipore Ltd., Cork, Ireland) on top of a glass microfiber filter (Whatman International Ltd., Maidstone, United Kingdom). For viral enumeration, each DNase-treated solution was filtered through a 0.22-µm Millex syringe filter (Merck Millipore Ltd., Tullagreen, Co., Cork, Ireland) to remove bacteria. Viruses were captured by filtering each viral filtrate through a 0.02-µm-pore-size Whatman Anodisc filter (Whatman International Ltd., Maidstone, United Kingdom). The filtration of bacteria and viruses was performed in a Millipore vacuum manifold (Sigma-Aldrich Corporation, St. Louis, MO, United States) under a pressure of less than 62 kPa. The filters containing bacteria and viruses were stained using SYBR gold in a final 2× concentrate (5,000-fold dilution of the original stock solution). The filters were immediately analyzed using a Nikon Eclipse E600 epifluorescence microscopy equipped with FITC filter set and Retiga EX-i CCD camera (Qimaging, Surrey, BC, Canada). Enumeration of bacteria and viruses were performed using IPlab software (BD Biosciences, Franklin Lakes, NJ, United States). For counting bacteria or viruses in each sample, at least 10 fields were digitally photographed at a magnification of 1,000X and an average of bacterial/viral counts was used for further calculation of their density in the sample. The bacterial and viral abundances in each treatment were derived from triplicate samples.

In vitro Prophage Induction With AHLs as Inducing Agents
The pooled and homogenized bacterial extracts was distributed into 30 aliquots in tubes each containing the respective inducing agents for induction assay or controls (Supplementary Figure S1). Sterile soil extract was used as the medium for induction assays to provide native growth substrates and trace elements to support bacterial metabolism and viral reproduction in induced lysogenic cells.
Eight different AHLs, (N-butyryl-, hexanoyl-, βketocaproyl-, heptanoyl-, octanoyl-, 3-oxododecanoyl-, tetradecanoyl-, and 3-oxotetradecanoyl-homoserine lactones, Wellington and Greenberg, 2019) varying in molecular weight and structure ( Table 1 and Supplementary Figure S2) and mitomycin C were selected for induction assays. The presence and production of these AHL molecules in soil has been demonstrated by in situ characterization of AHL production from both soil bacterial isolations and indigenous soil communities (Burmølle et al., 2005;Huang et al., 2013;Erdönmez et al., 2017). The selected AHLs were designated AHL1 to 8 based on their molecular weight from lowest to highest (acyl C chain length refer to Table 1). Each AHL was dissolved in ethyl acetate acidified with acetic acid (0.1%, vol/vol) as stock solution and applied at a final concentration of 1 µM in the bacterial suspension. The required amount of each AHL compound in stock solution was added into glass tubes, and the tubes were gently shaken for evaporation of the solvents so that the AHL compounds bonded to the tube bottom as a film (Ghosh et al., 2009). The pooled bacterial extract was distributed into AHLs-and mitomycin C-coated glass tubes, each containing a 10 ml aliquot of bacterial suspension. Ten ml of the exact same bacterial suspension was also distributed to clean glass tubes lacking any inducing agent to serve as control. Each treatment and control were prepared in triplicate. All suspensions were incubated in the dark for 18 h at room temperature (Ghosh et al., 2009). The viruses and bacteria in the suspensions were enumerated using epifluorescence microscopy to determine the induction response due to each inducing agent. Viral and bacterial abundance in each sample was estimated by epifluorescence microscopy enumeration as previously described (Williamson et al., 2007;Liang et al., 2019b).

Bacterial 16S rRNA Genes Sequencing and Statistical Analysis
After 18 h incubation, 1 ml of each cell suspension from the induction assays was transferred to a new sterile centrifuge tube and treated with DNase I for 20 min. The reaction was terminated by addition of EDTA prior to centrifugation at 5,000 × g for 20 min at 4 • C. The bacterial pellets were washed twice with sterilized extraction buffer to remove all lysed bacterial cells and any residual undigested free DNA. The genomic DNA of un-lysed bacterial cells that survived prophage induction from each sample was extracted using PowerLyser PowerSoil DNA isolation kit (Qiagen, Hilden, Germany) and quantified with using a Nanodrop one spectrophotometer (Thermo Scientific). The DNA samples were sent to the Genomics Core Laboratory at University of Tennessee (Knoxville, TN, United States) for sequencing. The V3-V4 region of 16S rRNA genes were amplified using PCR primer set (341F_CCTACGGGNGGCWGCAG, and 785R_GACTACHVGGGTATCTAATCC) for library construction, and finally sequenced via 300PE (paired-end) on the Illumina MiSeq platform (Illumina, United States) by using the manufacturers' protocol.
The obtained sequence data was demultiplexed depending on the barcodes. The raw 16S rRNA gene sequence data with all sequence reads was processed using the MOTHUR v.1.40.0 pipeline according to the MiSeq SOP (Kozich et al., 2013). In general, two sets of sequence reads (forward and reverse) for each sample were merged into contigs in which process the pairs of sequences were aligned with quality score of each base calculated. The total sequence reads from each treatment ranged from 57,608 to 213,551. To reduce sequencing and PCR errors, the resulting 3,821,672 sequences were trimmed of sequences containing ambiguous bases and long polymers. As many of the quality-filtered sequences were identical to each other, the unique sequences were identified and sorted with abundance. Sequences were classified by alignment to the customized SILVA reference (Quast et al., 2013), which pre-clustered sequences by 99% similarity. Following reference alignment, the removal of chimeric sequences was performed by VSEARCH algorithm. The divided sequences classified into OTUs (operational taxonomic unit) at 97% nucleotide identity by using the Bayesian classifier, which resulted in a total of 56,996 OTUs. Statistical analyses of the processed sequence data were performed using software R version 3.6.1 packages phyloseq (McMurdie and Holmes, 2014), vegan (version 2.5-2, Oksanen et al., 2016), ggplot2 (Wickham, 2016), and DESeq2 (Love et al., 2014). Bacterial taxonomic composition and alpha-diversity were calculated, and beta-diversity was also assessed based on Bray-Curtis dissimilarity matrix. The quantification and statistical inference of systematic differences of the bacterial taxonomic composition between each induction assay and the control samples were performed using the package DESeq2 (Love et al., 2014). Variance mean was estimated over all treatments (each was based on three replicate samples and treated independently and identically). The files of raw sequences were archived at the National Center for Biotechnology Information Databases (Sequence Read Archive) and can be obtained under accession number SRR10238150.

Prophage Induction
Bacteria were extracted, purified and concentrated from the agricultural soil in a way that eliminated most extracellular viruses as described above for bacterial extraction. In this way, the viral background counts were reduced from 6.3 × 10 8 to 5.3 × 10 5 ml −1 in the pooled bacterial suspensions. The resulting bacterial suspensions were pooled to create one homogeneous mixture of native soil bacteria for use in prophage induction assays. The pooled cell concentrates (before the 18 h incubation assay period) had a bacterial abundance of 1.28 × 10 8 cells ml −1 , and the density of viruses was 5.33 × 10 5 particles ml −1 . After the 18 h induction assays, the induced cell suspensions were compared directly to the uninduced control suspensions for variations in viral and cell abundance and bacterial community composition. Viral and bacterial abundance was quantified to assess the prophage induction response of host cells. The mitomycin C-treated suspensions had significantly lower bacterial abundance and notably higher viral abundance than the control suspensions (P < 0.01, t-test, Figure 1). No significant decrease in bacterial abundance was observed in the eight AHL-treated suspensions compared to the control samples. The viral abundance in AHL1-, AHL2-, AHL3-, AHL4-, AHL5-, or AHL7-treated suspensions was significantly higher compared with that in the control suspensions (P < 0.05, t-test, Figure 1B).

Bacterial Community Diversity and Composition
To evaluate the impacts of prophage induction at the community level, the bacterial species richness and alpha diversity were estimated based on species number and Shannon index, respectively as determined by analysis of 16S rRNA gene sequencing of un-lysed cells remaining following the induction assays. The results of 16S rRNA gene sequencing of each treatment were based on three replicate assays, which were treated independently and identically. The bacterial community in mitomycin C-treated suspensions had significantly higher species richness and diversity than that in the control suspensions (P < 0.01, t-test, Figure 2). In contrast, the bacterial species diversity in AHLs-treated suspensions was lower than that in the control suspensions (P < 0.05, t-test, Figure 2A). Similar results were obtained for inverse Simpson index (data not shown). There was no statistical difference in bacterial species richness between each AHL treated cell suspension and the control suspensions ( Figure 2B).
Principal component analysis (PCA) was used to visualize the dissimilarities of bacterial community structure between suspensions after prophage induction and showed that the communities were separated according to the specific AHL signal compound added to the cell suspensions ( Figure 2C). The bacterial community structure in the induction assays using AHL1-5 were notably different from those in the uninduced controls. Some dissimilarity was also observed in bacterial community structure among samples treated with different AHLs. For example, communities resulting from treatment with AHL1-5 were distinctly different from communities treated with AHL6-8 ( Figure 2C). mitomycin C-treated suspensions were clearly separated from controls and AHL-treated suspensions, and the first two principal coordinates explained 93.3% of the variation of bacterial community structure ( Figure 2C).

Bacterial Taxonomic Profiles Following Prophage Induction
To further examine the differences in bacterial community structure resulting from AHL-mediated prophage induction, we compared the abundance of each bacterial taxonomic group in every induction treatment to that in the control samples at both class and genus level. A stacked bar chart was constructed to show the community composition of the dominant bacterial classes across all groups, and notable differences were revealed between each treatment and the control Frontiers in Microbiology | www.frontiersin.org Control indicates unamended control cell suspensions without AHLs or Mitomycin C treatment. Each bar represents the mean of triplicate assays (n = 3), and the error bars show one standard deviation. Statistical significance between the treatment and the control group is indicated by * p < 0.05, * * p < 0.01, and * * * p < 0.001. group (Figure 3). The AHL1 treatment led to a decreased relative abundance of the greatest number of taxa including: Actinobacteria, Alpha-Proteobacteria and Acidobacteria Gp6 at the class level (P < 0.01; Figure 3 and Supplementary Figure S3) and Microbacterium, Actinophytocola, Mamoricola, Nocardioides, Novosphingobium, Lysobacter, Pandoraea, and Sphingomonas at the genus level (P < 0.01; Figure 4). AHL2 treatment led to a decreased abundance of two bacterial classes, Actinobacteria and Alpha-Proteobacteria, and genera, Brevundimonas, Mesorhizobium, Lysobacter, and Sphingomonas (P < 0.01; Figure 4). AHL5 treatment caused a decreased abundance of the second greatest number of bacterial taxa including Acidobacteria Gp6, Gp17 and Actinobacteria at the class level (P < 0.01; Figure 3 and Supplementary Figure S3) and Mamoricola, Gp17, Bosea, Nocardioides, Bradyrhizobium, Novosphingobium, and Acidobacteria Gp6 at the genus level (P < 0.01; Figure 4). Treatment with AHL3 and 4 resulted in a decreased abundance of only one class, Chlamydiae. However, AHL3 and 4 had distinct induction effects at the genus level where the abundance of Streptomyces, and unclassified Anaerolinaceae and Parachlamydiaceae declined with AHL3 treatment and Bosea, Verrucomicrobium, and Pedobacter in AHL4 exposure (P < 0.01). The AHLs discussed above also resulted in an increased abundance of some bacterial taxa, e.g., classes of Gamma-Proteobacteria, Bacilli and Flavobacteria (P < 0.01). It is interesting to note that an increased abundance of Pseudomonas was observed in the treatment with AHL1-5 that had significant shifts in the density of specific bacterial groups. The bacterial genera having increased relative abundance after treatment with some AHLs also included Filimonas, Cellvibrio, Duganella, Pelomonas, and Flavobacterium and certain other unclassified genera (P < 0.01). The bacterial taxonomic profiles in the treatment of AHL6, 7, and 8 had no statistical difference compared with the control cell suspensions. Mitomycin C treatment stimulated much broader impacts on the taxonomic composition of the bacterial community than the AHLs suggesting that the prophage induction response brought about by exposure to mitomycin C was less specific than any of the AHLs used as inducing agents. This is consistent with the more general prophage induction response likely brought about by activation of DNA-repair systems associated with mitomycin C-mediated prophage induction. As mitomycin C is an antibiotic, it likely also directly inhibits bacterial growth and causes more cell death leading to the more extensive apparent impacts on the bacterial community composition compared with AHLs treatment. The abundance of six bacterial classes, including Flavobacteria, Bacilli, Gamma-Proteobacteria, Acidobacteria Gp7, Opitutae and Verrucomicrobiae, declined in mitomycin C treated compared to the control suspensions (P < 0.01; Figure 5). Up to 53 genus-level bacterial groups, such as Aeromonas, Flavobacterium, Albidiferax, Chitinimonas, Arthrobacter, Chryseobacterium and Paenibacillus, had decreased relative abundance after mitomycin C treatment compared with the control suspensions (P < 0.01; Supplementary Figure S4). Cellvibrio, Filimonas, Pelomonas, Flavobacterium, and Pseudomonas, the bacterial genera that had increased relative FIGURE 4 | Genus-level differences of the bacterial community composition between each induction assay of N-(Butyryl, Hexanoyl, β-Ketocaproyl, Heptanoyl, Octanoyl, 3-Oxododecanoyl, Tetradecanoyl, and 3-Oxotetradecanoyl) homoserine lactones (indicated by AHL1-8, respectively) and the control suspensions. Only statistically significant differences (P < 0.01) were shown. The direction of bars represents decreases (left) or increases (right) in relative abundance of the specific bacterial taxonomic group after the induction assays relative to the unamended controls. Differential analysis of bacterial taxonomic abundances between the treatment and control group was performed using the DESeq2 package for R (Love et al., 2014). Variance mean was estimated over all treatments (each was based on three replicate assays and treated independently and identically). The color of the bar represents taxonomic group shared by two or more different treatments, and the bars without fill represent the taxonomic groups with significant differences compared with the control only detected in the current treatment.
abundance after the treatment of some AHLs, had decreased relative abundance in the mitomycin C treatment (P < 0.01). The common bacterial genera having decreased relative abundance (comparison between the treatment and control group) and shared by mitomycin C and some AHL treatments included Novosphingobium and Pedobacter. Mitomycin C exposure also resulted in a slightly increased proportion of 44 bacterial genera (P < 0.01; Supplementary Figure S4).

DISCUSSION
This study revealed the lysogeny-lysis switch of some temperate phages is responsive to QS autoinducers and the phenomenon may be more widespread than the handful of well-characterized phage-host systems reported. The molecular basis of a host QS autoinducer controlling a phage lysogeny-lysis decision has been recently characterized (Silpe and Bassler, 2019a).
FIGURE 5 | Systematic differences of the bacterial taxonomic composition between induction assay of mitomycin C and the control samples at Genus levels. Only statistically significant differences (P < 0.01) are shown. The direction of bars represents decreases (left) or increases (right) in relative abundance of the specific bacterial taxonomic groups after the induction assays.
Frontiers in Microbiology | www.frontiersin.org Subsequent studies reported phage responses to other types of host autoinducers (Laganenka et al., 2019;Silpe and Bassler, 2019b). Bacteria communicate with their close relatives using specific QS signals, thus phage-bacterium QS connections may also be molecular structure dependent (Liang and Radosevich, 2019). Building on this idea, we hypothesized that any single QS signal should only induce prophages within a small subset of closely related host bacteria. In this study, eight AHLs of varying molecular weight and structure were selected for evaluation of phage responses of lysogeny-lysis switching and its significance in structuring the exposed well mixed bacterial communities extracted and purified from a single soil source in vitro.
Significant differences in viral abundance were observed in the treated microbial cell suspensions exposed to AHL1, 2, 3, 4, 5, and 7 compared with the control suspensions. The increase of viral abundance in AHL-treated suspensions was consistent with a burst of viral production from prophage induction, especially in a relatively few host taxa, and was in agreement with our hypothesis. Other recent studies also reported prophage induction by QS molecules in bacterial hosts, such as Enterococcus faecalis (Rossmann et al., 2015) and E. coli (Laganenka et al., 2019), leading to viral production and bacterial lysis. Interestingly, Oh et al. (2019) discovered that though different prophages were induced by environmental cues they achieved differential induction responses indicating a unique governing system for each prophage. A significant decrease in total cell abundance was typically not observed upon exposure to AHLs. This too is consistent with the hypothesis if each AHL triggered prophage induction in a relatively narrow range of perhaps less numerically abundant taxa. In order to test this hypothesis in future studies, investigating bacterial community responses under addition of cocktails of multiple AHLs to the same bacterial suspensions should be considered. Additionally, meta-genomic and meta-transcriptomic analyses during the induction assay could reveal specific prophagehost responses to AHL-exposure. Two other possibilities may also be consistent with the observed results of the treatmentdriven changes in viral and bacterial abundances. One possibility is that some of the increased abundance of fluorescence particles from AHL-treated cell suspensions may have been DNA-carrying membrane vesicles secreted/released by AHLtreatment (Toyofuku et al., 2017). Previous reports have indicated that membrane vesicles play important roles in QS and can increase epifluorescence-based viral counts (Li et al., 2016;Biller et al., 2017). The other possibility may be that the produced phages during treatment were released without host lysis (Rakonjac et al., 1999;Xue et al., 2012;Loh et al., 2019). Therefore, future work is needed to elucidate the effects of QS on viral and bacterial populations and advance our understanding of this important phenomenon which may have important ecological consequences in nature. Other experimental approaches that will be useful in confirming the findings reported here might include qPCR quantification of bacterial 16S rRNA gene abundance and phage marker genes, transmission electron microscopy enumeration and morphological description of viral particles.
Even a small collection of lysogenic bacterial taxa triggered to enter the lytic cycle could result in changes in taxonomic composition. We inspected the bacterial community diversity based on species richness and Shannon index. Significantly lower Shannon indices but no significant changes in species richness were observed after AHL treatment. These results suggest that AHL treatment decreased species evenness potentially by suppressing a subset of bacterial species and resulting in increased relative abundance of some other species. Further analysis of Pielou's evenness across the samples showed that the species evenness was significantly lowered in AHL1 (P < 0.05), AHL3 (P < 0.05), AHL5 (P < 0.01), and AHL7 (P < 0.01) compared to the control suspensions which further supported the proposed mechanism of AHL treatment influencing the bacterial community structure. In contrast, mitomycin C treatment, as a broad-spectrum inducing agent, also possessing broad toxicity, commonly used for prophage induction (Williamson et al., 2007;Knowles et al., 2017;Liang et al., 2020), prompted an increase in both species richness and diversity relative to untreated controls. Mitomycin C treatment brought about a significant increase in viral abundance (P < 0.001) and notable decrease in bacterial abundance (P < 0.01) through a wide range of viral lysis, thus contributing to higher richness and evenness in bacterial taxonomic profiles. The measurement of Pielou's evenness in these suspensions demonstrated that mitomycin C treatment significantly increased the species evenness compared to the control group (P < 0.01).
Since viral production has been shown to be correlated to nutrient cycling and bacterial metabolism, the prophage induction can also contribute to resource redistribution and bacterial community composition (Liang et al., 2019c,b;Oh et al., 2019). To further examine the putative effects of AHL-dependent prophage induction on bacterial community structure, we determined changes in relative abundance (expressed as log2fold changes relative to control cell suspensions) of the affected taxa. Up to 10 bacterial genera decreased in relative abundance after any single AHL treatment, and a total of 23 different genera for all AHL treatments combined. Decreased relative abundance of Lysobacter, Novosphingobium, Sphingomonas, Bosea, and Nocardioides was observed in at least two AHL treatments. The significant decrease in density of the AHL-targeted bacterial genera, e.g., Lysobacter, Novosphingobium, and Sphingomonas, suggests AHL-directed transition of prophages from lysogeny to lysis in these bacteria. While AHLs caused reduced relative abundances of some bacterial taxonomic groups, AHL treatment also resulted in an increased relative abundance of some bacterial taxa, e.g., classes of Gamma-Proteobacteria, Bacilli and Flavobacteria although the positive affect on some taxa may be indirectly attributed to prophage induction and lysis of other susceptible taxa and/or some growth of the non-susceptible groups. However, if the observed increases in relative abundance were strictly due to growth of certain taxa on the extractable soil dissolved organic carbon (DOC) that was in the induction assays then the patterns would have been similar in all the treatments instead of the differing patterns that were observed. Interestingly, the relative abundance of Nocardioides and Bacilli, gram-positive genera, varied significantly in the presence of AHL, a QS autoinducer in gram-negative bacteria which may be hard to explain. However, previous studies have shown that many of Nocardioides and Bacilli members can degrade AHLs leading to quorum quenching (Dong et al., 2002;Yoon et al., 2006;Hong et al., 2012;Vinoj et al., 2014;Kusada et al., 2019) thus likely having undetermined important associations with AHLproducing bacterial groups. Microbial interactions within the community are also an important factor to be considered for variations in community composition and diversity (Tian et al., 2018). However, the explanations of abundance variations of Nocardioides and Bacilli cannot be demonstrated by the current results in this research.
Quorum sensing has been demonstrated to be widespread among bacteria (Polkade et al., 2016;Liao et al., 2018;Ling et al., 2019), however, AHLs have not been shown as QS signals for many of the AHL-impacted bacterial groups in the present study, such as Nocardioides, Streptomyces, Pedobacter, and Verrucomicrobium. While we consider the AHL-mediated prophage induction as the main driving force of differences in the bacterial community composition we observed, other factors should also be considered. For example, QS autoinducers can regulate bacterial collective behaviors such as virulence and biofilm formation and influence inter-and intra-population interactions within the community which thus may influence bacterial community structure (Kim et al., 2016). There are also reports showing that AHLs can be utilized by some bacteria for growth (Huang et al., 2003). So, the selective impact of AHLs themselves needs to be further considered.
Mitomycin C treatment resulted in decreased abundance of 53 bacterial genera. Though the observed decreases may have resulted from virus-mediated host cell lysis upon prophage induction or simply from direct toxicity of the mitomycin C, these broad-spectrum changes clearly contributed to the observed increases in evenness of bacterial species profiles and thus increases in the community diversity. Increased proportion of specific bacterial groups observed in AHL or mitomycin C treatment might be derived from competitive release (Loudon et al., 2014). Growth of some rare bacterial species might also be promoted which contributed to the increased species richness in mitomycin C treated suspensions. Susceptible bacterial species were lysed by chemical induction of prophages allowing the remaining competitors to utilize the resources more fully, and the remaining members of the bacterial community may also have access to the nutrients released through the viral shunt (Danovaro et al., 2008;Kuzyakov and Mason-Jones, 2018).
Lysogeny has been demonstrated to be widespread and a common viral life strategy in nature and shown to have links with the dynamics of the nutrient regime and host density (Howard-Varona et al., 2017). Chemical induction assays, like mitomycin C, have been adopted to assess lysogeny among viral and bacterial communities (Ghosh et al., 2009;Knowles et al., 2017;Kronheim et al., 2018). Though the phage-bacterium connections through QS was recently discovered, the QS based prophage-induction are largely unknown except for a few phage-host pairs none of which were derived from soil, with little known about the influence of prophage induction on microbial community dynamics. In this study, we focused on the inducing effects of eight AHLs among different well-known bacterial autoinducers on microbial community structure. Our findings revealed that a broad range of bacterial taxonomic groups were putatively susceptible to prophage induction by these host autoinducers (i.e., AHLs), and we also demonstrated how transitions from lysogeny to lysis of temperate phages responsive to different host autoinducers can have pivotal roles in influencing bacterial community structure. This research provides theoretical and methodological foundation for future study of phage-bacterium communication and the lysogeny-lysis switch of soil viral communities. For future study, the results reported here should be further investigated by including analysis of lysed bacteria as template for 16S rRNA sequence analyses. In addition, metagenomic and metatranscriptomic analyses before during and after induction assays may reveal more specific and direct evidence supporting QS-controlled lysogeny-lysis switching and the hypothesis addressed in this study may resolve more statistically robust relationships and provide unique highresolution views of virosphere responses to host autoinducers.

DATA AVAILABILITY STATEMENT
The datasets generated for this study can be found in the National Center for Biotechnology Information Databases and can be obtained under accession number SRP224599.