The global anaerobic metabolism regulator fnr is necessary for the degradation of food dyes and drugs by Escherichia coli

ABSTRACT The microbiome is an underappreciated contributor to intestinal drug metabolism with broad implications for drug efficacy and toxicity. While considerable progress has been made toward identifying the gut bacterial genes and enzymes involved, the role of environmental factors in shaping their activity remains poorly understood. Here, we focus on the gut bacterial reduction of azo bonds (R-N = N-R’), found in diverse chemicals in both food and drugs. Surprisingly, the canonical azoR gene in Escherichia coli was dispensable for azo bond reduction. Instead, azoreductase activity was controlled by the fumarate and nitrate reduction (fnr) regulator, consistent with a requirement for the anoxic conditions found within the gastrointestinal tract. Paired transcriptomic and proteomic analysis of the fnr regulon revealed that in addition to altering the expression of multiple reductases, FNR is necessary for the metabolism of L-Cysteine to hydrogen sulfide, enabling the degradation of azo bonds. Furthermore, we found that FNR indirectly regulates this process through the small noncoding regulatory RNA fnrS. Taken together, these results show how gut bacteria sense and respond to their intestinal environment to enable the metabolism of chemical groups found in both dietary and pharmaceutical compounds. IMPORTANCE This work has broad relevance due to the ubiquity of dyes containing azo bonds in food and drugs. We report that azo dyes can be degraded by human gut bacteria through both enzymatic and nonenzymatic mechanisms, even from a single gut bacterial species. Furthermore, we revealed that environmental factors, oxygen, and L-Cysteine control the ability of E. coli to degrade azo dyes due to their impacts on bacterial transcription and metabolism. These results open up new opportunities to manipulate the azoreductase activity of the gut microbiome through the manipulation of host diet, suggest that azoreductase potential may be altered in patients suffering from gastrointestinal disease, and highlight the importance of studying bacterial enzymes for drug metabolism in their natural cellular and ecological context.

W hile it has long been appreciated that the gut microbiota, the trillions of microor ganisms found in the gastrointestinal (GI) tract, and its aggregate genomes (the gut microbiome) contribute to the digestion and metabolism of dietary macronutrients, the broader role of the microbiome in the metabolism of xenobiotics (diet-derived and pharmaceutical small molecules) is less well understood.Recently, work in cell culture, mice, and humans has emphasized that both excipients (food and drug additives) and pharmaceuticals are extensively metabolized by the gut microbiota (1)(2)(3)(4)(5)(6).The presence of an azo bond (R-N = N-R') is shared between both pharmaceuticals and excipients (7,8), which is notable given the ability of diverse human gut bacteria to reduce azo bonds (9).While azoreductase activity may be a core metabolic function of all human gut microbiotas (10)(11)(12)(13)(14), the enzymatic and nonenzymatic mechanisms responsible for this activity and their sensitivity to environmental factors remain poorly understood.
Mechanistic insights into this process are essential given the broad impact of azo reduction for antibiotics (15) and antiinflammatory drugs (7,16).Furthermore, the consumption of food, drug, and cosmetic (FD&C) dyes is increasing (17), providing additional substrates for gut bacterial metabolism.We recently discovered that FD&C dyes are potent inhibitors of the mammalian influx transporter OATP2B1, interfering with drug absorption in mice (18).This effect was rescued by human gut bacterial metabolism due to an inability of the downstream microbial metabolites to inhibit OATP2B1 (18).Additional work in rodent models has implicated FD&C dyes in carcinogenesis (19)(20)(21) and inflammatory bowel disease (22).Thus, the ability to predict or control azoreductase potential within the human gut microbiome could have broad implications for host health and disease.
The canonical enzyme implicated in this metabolic activity is azoreductase (AzoR) (9,11,14,23).Work on the model human gut bacterium Escherichia coli has demon strated that the purified AzoR protein is sufficient to reduce azo bonds (24).Other bacteria, for example, Pseudomonas aeruginosa (25)(26)(27)(28) and Enterococcus faecalis (29)(30)(31) encode multiple azoreductases.Alternative mechanisms have been proposed for azoreductase activity, including the electron transport chain of Shewanella oneidensis, nicotinamide adenine dinucleotide (NADH), and hydrogen sulfide (H 2 S) (32)(33)(34).While these studies provide valuable mechanistic insight, a major limitation is their focus on in vitro biochemistry, neglecting to address the cellular and genetic mechanisms that impact azoreductase activity or the potential confounding effects of the complex physiological and microbiological interactions within the GI tract.
We sought to address this knowledge gap through the mechanistic dissection of a representative member of the human gut microbiota.Surprisingly, we found that the azoR gene is dispensable for the azoreductase activity of E. coli.In our search for alternative azoreductase enzymes, we discovered that E. coli azo reduction is regulated by the oxygen-sensing dual-transcriptional regulator, FNR.In turn, we show that FNR, via small regulatory RNA fnrS, regulates L-Cysteine metabolism, which produces H 2 S. Surprisingly, E. coli-produced H 2 S is sufficient to reduce azo dyes and drugs.Taken together, we demonstrate that environmental oxygen regulates the metabolism of L-Cysteine to H 2 S, which degrades commonly consumed azo-bonded compounds.Our results highlight the importance of studying and developing an understanding of the environmental context and regulation of key microbial metabolisms beyond annotated enzyme function.

E. coli AzoR is dispensable for azo dye depletion under anaerobic conditions
Prior work on azoreduction in the context of the human gut and industrial settings has focused on anaerobic conditions (11,23).We compared the depletion of a representative azo bond-containing food coloring (FD&C Red No. 40) by E. coli K-12 BW25113 under aerobic and anaerobic conditions, revealing that oxygen interferes with this activity (Fig. 1A).L-Cysteine, a commonly used supplement in anaerobic microbiology, was also required for FD&C Red No. 40 depletion by E. coli (Fig. 1B).Unless otherwise noted, these optimized culture conditions [Luria Broth (LB) medium plus 0.05% L-Cysteine] were used for our whole-cell azo dye depletion experiments.
Previously, we (18,35) and others (36) had identified genes homologous to E. coli azoreductase (AzoR) in the genomes of diverse human gut bacteria capable of depleting drugs and other azo bond-containing compounds.We had originally sought to establish a tractable system to study these heterologously expressed genes by first abolishing this activity in E. coli.Much to our surprise, the E. coli Keio collection azoreductase knockout strain (ΔazoR::kan R ) had activity indistinguishable from wild-type (wt) during growth in both solid (Fig. 1C) and liquid (Fig. 1D) media containing FD&C Red No. 40.We validated this finding by constructing and testing a clean deletion of azoR (Fig. 1E), which was confirmed by gel electrophoresis (Fig. S1A) and Sanger sequencing (Fig. S1B).The ΔazoR strain did not differ from wt in whole-cell dye depletion kinetics (Fig. 1F).
Given these unexpected results, we wondered if AzoR might be inactive in the tested E. coli strains, potentially due to mutations that could have occurred during the creation of our lab's culture collection (37).Sanger sequencing of the azoR gene verified that the coding sequence was 100% identical to the deposited E. coli BW25113 genome.To further validate the activity of this enzyme, we purified the Histidine-tag ged AzoR protein from E. coli AG1 with the pCA24N plasmid for overexpression (ASKA collection) (38).Cells were grown to mid-exponential growth phase, azoR expression was induced with isopropyl β-D-1-thiogalactopyranoside (IPTG), then cells were harvested and protein was purified using a nickel-charged affinity column.The fraction containing the most protein was determined using an SDS-PAGE gel and the size of the protein was consistent with AzoR (Fig. S2A).As expected, the purified protein was sufficient to clear a panel of azo bond-containing dyes and drugs (Fig. S2B).Additionally, we were able to create a Michaelis-Menten curve describing the enzymatic activity against FD&C Red No. 40, where we calculated a K M of 67.7 µM (Fig. 1G).Finally, we tested our whole-cell assay using a panel of 10 food additives and pharmaceutical compounds containing azo bonds (Fig. 1H), demonstrating that azoR is unnecessary for the depletion of multiple dyes using the growth conditions established in Fig. 1A and B. These results confirm that the canonical AzoR enzyme in E. coli (24) is sufficient for azoreductase activity; however, azoR is not required for azo bond reduction by E. coli cells.These results emphasize that biochemical activity of purified proteins is not necessarily predictive of metabolism in whole cells and motivate a renewed search for the mechanism(s) responsible.

FNR is necessary for the depletion of azo dyes by E. coli
To identify genes necessary for azo dye depletion, we screened 113 candidate gene deletions from the Keio collection using a high-throughput liquid media assay for the depletion of the widely used food coloring FD&C Red No. 40.We focused on anaerobic reductases based on the requirement for anaerobic growth (Fig. 1A).We also included multiple genes based on prior studies linking them to azo dyes: genes upregulated by azo dyes in the presence of a quinone (39) or S. oneidensis MR-1 electron transport chain that is involved in azo dye depletion (40).Gene annotations were cross referenced on EcoCyc (41) (Table S1).Consistent with our previous data (Fig. 1), ΔazoR was comparable to wt (Fig. 2A; Table S1).A stringent cutoff of 95% confidence interval (CI) of all genes tested revealed 22 gene deletions that had significantly impaired activity (Fig. 2A).The most dramatic loss-of-function was found for Δfnr::Kan R (fumarate and nitrate reduction regulator) resulting in a 2.3-fold decrease in azo dye depletion.We validated the role of fnr in FD&C Red No. 40 depletion by constructing a clean deletion (Fig. 2B) and confirming the deletion by gel electrophoresis (Fig. S1C) and Sanger sequencing (Fig. S1D).Activity was restored through complementation of fnr on the pCA24N plasmid from the ASKA collection (38) (Fig. 2C).
FNR plays a key role in oxygen sensing and the transition to anaerobic metabolism (42,43), explaining the previously observed requirement for anaerobic conditions (Fig. 1A).The requirement for fnr was not unique to FD&C Red No. 40; significantly decreased dye depletion was observed across a panel of seven food dyes and three azo bondcontaining drugs (Fig. S3).Overexpression of AzoR was insufficient to rescue Δfnr (Fig. 2C), suggesting that other aspects of the vast fnr regulon (44)(45)(46)(47) are responsible for the observed phenotype.Permeability issues likely result in inadequate dye reaching AzoR in the cytosol (48).

Paired transcriptomics and proteomics reveals additional genes necessary for azo dye depletion
FNR has a large regulon that has been extensively characterized by transcriptomics (RNAseq) and chromatin immunoprecipitation sequencing (ChIP-Seq) (43,44,46,49).However, only 14 of the 55 genes found within the FNR regulon included in our genetic screen led to a significant loss-of-function at our 95% CI (Fig. 2A; Table S1), indicating that only a subset of the FNR regulon is necessary for azo dye depletion.Prior to testing additional candidate genes, we sought to confirm and extend prior work on the FNR regulon.Our culture conditions differ from past experiments in three ways that could alter gene expression and/or metabolic activity: (i) the presence of L-Cysteine; (ii) the inclusion of the azo dye FD&C Red No. 40; and (iii) the use of the vehicle DMSO, which is metabolized by E. coli under anaerobic conditions (50,51).We also wanted to assess gene expression during stationary phase, given that FD&C Red No. 40 levels do not decrease until the exponential phase of growth (Fig. S4).Finally, we sought to build upon prior work on the FNR regulon by also assessing protein levels using proteomics.
We grew wt and Δfnr E. coli in LB media with 0.05% L-Cysteine and inoculated the cultures with 250 µM FD&C Red No. 40 or DMSO vehicle at either mid-exponential or stationary growth phases (Fig. 3A).After a 40 min incubation period, samples were collected and split for paired transcriptomics (RNA-seq) and proteomics analysis.In total, we generated 12.4 ± 6.2 million high-quality reads/sample (RNA-seq; Table S2) and 881 ± 23 unique peptides after matching spectra through the MS-GF+ database (52) and adjusting for false discovery rate (FDR) of 0.01 against the MS-GF+ generated decoy database (proteomics; Table S3).
Principal coordinates analysis revealed a clear separation of profiles between strains at the RNA (Fig. 3B) and protein (Fig. 3C) levels.As expected (53)(54)(55), we also saw a clear difference in expression profile between growth phases; however, the separation between strains was maintained at both timepoints (Fig. 3B and C).There was a more minimal impact of FD&C Red No. 40 (Fig. 3B and C).
The number of differentially expressed genes was consistent with the overall trends in these datasets, revealing marked strain-and growth-phase differences in transcript and protein levels, with a minimal impact of FD&C Red No. 40 (Fig. 3D and E; Table S4).Only a single gene was significantly upregulated by FD&C Red No. 40 (log 2 fold change >1 and p adj <0.1): ydeN, a putative sulfatase (56,57).The ΔydeN::KanR strain had comparable azoreductase activity to wt (Fig. S5).In the proteomics data set three features were enriched in the vehicle group (Fig. 3E).
Given the marked difference in azo dye depletion between the wt and Δfnr strains, we next focused on the 588 transcripts and 216 proteins significantly different between strains with a log 2 fold change>|1| (p adj <0.1; Fig. 3F through H).Eighty-three features were both significant and had a fold change in the same direction in both datasets (Fig. 3H).Notably, these consistent features included 62 enzymes and 18 [oxido]reductases (Table S5).azoR was slightly significantly upregulated in the wt strain (Fig. S6A), with a trend toward increased peptide intensity (Fig. S6B).As expected, fnr transcripts were undetectable in the Δfnr strain (Fig. S6C).Furthermore, 96/113 genes used in our azo dye depletion screen (Fig. 2A) had detectable expression (Table S1).
Given the size of the FNR regulon in both prior studies (43,44,46,49) and our datasets, including numerous enzymes that could feasibly impact azo dyes, we sought to further narrow down our candidate gene list by focusing on FNR-dependent regulators.Features that were significantly increased in both RNA-seq and proteomics in wt compared to Δfnr (log 2 fold change>|1|, p adj <0.1) were analyzed for direct regulators using EcoCyc (see Materials and Methods).This analysis revealed 11 regulators whose regulons were significantly enriched in our gene set (Fig. 3I).As expected, fnr was the most significantly enriched regulator; however, the proportion of the FNR regulon represented was lower than some of the other regulators due to the much larger size of the FNR regulon and our stringent criteria for identifying FNR-dependent genes.

L-Cysteine metabolism is required for depletion of azo dyes
L-Cysteine metabolism by E. coli generates hydrogen sulfide (H 2 S) as a dead-end metabolite (Fig. 4A) (63,66), which has been previously implicated in azo dye degrada tion (33,34,65).Thus, we hypothesized that fnr controls H 2 S production and thus azo bond degradation by increasing the expression of the cyuR-regulated cyuA and cyuP.First, we tested the requirements for L-Cysteine and H 2 S in azoreduction.To facilitate rapidly establishing anaerobic conditions, we routinely supplement our media with L-Cysteine (67)(68)(69).Removal of L-Cysteine impaired the ability of E. coli to clear FD&C Red No. 40 (Fig. 1B).L-Cysteine promoted FD&C Red No. 40 depletion in a dose-dependent manner (Fig. 4B), suggesting that L-Cysteine is a key substrate for azoreductase activity during in vitro growth under anaerobic conditions.Anaerobic conditions in LB media lacking L-Cysteine, equilibrated in the anaerobic chamber for 1 wk were confirmed by eye using the indicator resazurin.Notably, L-Cysteine alone was unable to deplete FD&C Red No. 40 (Fig. 4C), while addition of sodium sulfide demonstrated a dose-dependent reduction in sterile media (Fig. 4D).To test the relationship between L-Cysteine and H 2 S, we quantified H 2 S in whole-cell culture using the methylene blue assay (see Materials and Methods).H 2 S levels were below our limit of detection (1.3 µM) in the Δfnr strain, significantly less than the triple-digit µM levels produced by wt (Fig. 4E).We also quantified the amount of FD&C Red No. 40 depleted in these samples and saw the inverse trend (Fig. 4F).Consistent with the importance of L-Cysteine metabolism for dye depletion, peak H 2 S concentration corresponded to the decrease in FD&C Red No. 40 for wt with Δfnr stable over time (Fig. 4G).
Next, we investigated if H 2 S alone was the fnr-driven azoreductase mechanism for FD&C Red No. 40 depletion.We utilized single and multi-gene knockouts for the cyuR regulon.As previously mentioned, CyuR regulates both the L-Cysteine transporter (cyuP) and L-Cysteine desulfidase (cyuA) (63) (Fig. 5A).Consistent with our hypothesis, cyuA, cyuP, and cyuR were significantly decreased at the transcript level in Δfnr E. coli relative to wt controls (Fig. 5B; Fig. S6D).We validated these findings using RT-qPCR (Fig. S6E through G).The cyuP and cyuR knockout strains produced significantly less sulfide (Fig. 5C) and depleted significantly less FD&C Red No. 40 (Fig. 5D) in L-Cysteine-containing media.The combined deletion of cyuR and azoR did not significantly impact sulfide levels (Fig. 5E); however, there was a slight but statistically significant impairment in azo dye depletion (Fig. 5F).The addition of 0.05% L-Cysteine at mid-exponential phase increased the expression of cyuR, cyuA, and cyuP; cyuR and cyuP were induced signifi cantly higher in wt relative to Δfnr (Fig. S7A through D).

FNR indirectly regulates the cyu genes
Finally, we sought to better understand the mechanism through which FNR impacts the expression of cyuR and its target genes cyuAP.FNR was not detected at the cyuR promoter in chromatin immunoprecipitation experiments (44), consistent with the lack of the FNR binding motif TTGATnnnnATCAA (44).These results led us to hypothesize that FNR indirectly regulates the cyu genes.A review of the prior literature and the EcoCyc database identified multiple potential intermediates that are both under the control of FNR and affect cyu: the MarA and SoxS transcription factors (70) and the regulatory small noncoding RNA fnrS (71, 72) (Fig. 6A).MarA and SoxS both regulate cyuR (70), while fnrS is a proposed positive regulator of cyuP (72).Neither marA nor soxS complemented the Δfnr phenotype (Fig. S7E and F).Furthermore, the ΔmarA::Kan R and ΔsoxS::Kan R strains were not significantly different from wt (Fig. S7G).Together, these experiments suggest that marA and soxS are not involved in the fnr-dependent FD&C Red No. 40 depletion or H 2 S production.In contrast, we found that the small noncoding RNA fnrS is critical for these pheno types.We were able to complement the Δfnr dye and sulfide phenotypes with fnrS (Fig. 6B and C).While complementation was sufficient to rescue the Δfnr phenotype, the ΔfnrS::Kan R carrying the empty plasmid did not show a loss-of-function phenotype at our usual sampling time of 24 h (Fig. 6B and C).We hypothesized that this is likely due to alternative regulators activating cyuAP and other L-Cysteine metabolism genes.Thus, we measured FD&C Red No. 40 degradation and sulfide production in the knockout strain over time.fnrS::Kan R demonstrated a delayed onset of dye degradation and lower sulfide production compared to wt (Fig. 6D and E).Furthermore, fnrS expression was induced significantly more in response to L-Cysteine in wt relative to Δfnr E. coli (Fig. 6F).

DISCUSSION
Our results provide mechanistic insight into the degradation of azo dyes by human gut bacteria.Surprisingly, the canonical azoreductase enzyme of E. coli is dispensable for whole-cell azoreductase activity.Instead, the shift to anaerobic growth is critical.We dissected a pathway through which oxygen is sensed by FNR, leading to upregulation of the cyu operon, increased uptake of L-Cysteine, enhanced production of H 2 S, and nonenzymatic azo bond reduction (Fig. 7).These findings provide a cautionary tale for the ability to extrapolate results from biochemical studies of purified proteins to infer the metabolic activity of cells and complex microbial communities.Despite a rich literature on E. coli AzoR (24,48,73), we report that the azoR gene is dispensable for azo reduction by E. coli during anaerobic growth.These findings emphasize the importance of pairing biochemical and heterolo gous expression studies of enzyme function with methods for gene deletion in the original strain.Future work should extend this approach to the other bacterial species with biochemically characterized azoreductase enzymes (25,29,74,75).For species where genetic tools do not exist (the vast majority of human gut bacteria), one could leverage natural strain-level variation to test the association between gene presence and sequence variation with metabolic activity (37).
The global regulator FNR is essential for E. coli azoreduction, highlighting the importance of environmental oxygen for the gut bacterial metabolism of drugs and excipients.Confirming the prior work on FNR in different growth conditions (42,61), we confirmed that the FNR regulon is vast, including hundreds of genes involved in different aspects of anaerobic metabolism.Interestingly, efforts to utilize bacterial azoreduction to degrade dye runoff in water systems also require anaerobic conditions, suggesting that our findings may be generalizable to microbial communities outside of the human gut (23,76,77).Importantly, E. coli exhibits features of both aerobic and anaerobic growth within the GI tract (61,78), depending upon its physical niche and host pathophysiology (79,80).Our results suggest that inter-individual variations in the environment in which E. coli grows could have downstream consequences for its ability to act upon azo bondcontaining compounds.More work is necessary to assess how the redox potential of the gut lumen and other environmental factors shape gut bacterial azo reduction in the context of host health and disease.Surprisingly, FNR controls azoreductase potential through the reduction of azo bonds by H 2 S. Given that L-Cysteine is the major source of H 2 S for E. coli, these results provide an additional layer of environmental control for azoreduction.Differences in dietary protein, host proteoglycans, and/or L-Cysteine biosynthesis by other members of the gut microbiota could all potentially influence the activity of E. coli and its ability to act on azo bonds.Future studies leveraging gnotobiotic mice, defined dietary interventions, and isotope tracing could be informative in this regard.
Our results indicate that the L-Cysteine utilization regulator cyuR, among others, is under indirect FNR control, emphasizing the complex multi-level regulatory network contributing to bacterial azoreduction.The cyu operon is the main source for L-Cysteine import and metabolism by E. coli (63,81).We also found that H 2 S production requires cyuR.Interestingly, the L-Cysteine uptake transporter, CyuP, was more important than the enzyme CyuA, suggesting that the import of L-Cysteine could provide substrate for additional pathways for H 2 S production.CyuA only functions anaerobically (63), consistent with our observations that azoreduction only occurs anaerobically and is regulated by FNR and CyuR.Moreover, the key role of L-Cysteine import suggests that amino acid availability from diet and/or host tissues could be critical for gut bacterial azoreduction.Future studies with mouse models and various dietary amino acid interventions could provide an opportunity to modulate microbiome azoreductase activity.
Deletion of both cyuR and azoR did not fully abolish E. coli azoreduction.Thus, additional pathways for L-Cysteine metabolism and/or as-of-yet undiscovered enzyme(s) may contribute to the metabolism of azo dyes by E. coli.In our initial screen, we found 21 deletion strains with a partial loss-of-function phenotype, all of which encode enzymes theoretically capable of azo reduction.Of these, four were also identified in our paired proteomic and transcriptomic data set (Table S6), which could be viable candidates for further biochemical characterization and/or complementation screens in the Δfnr background.
H 2 S production is widespread throughout gut bacterial species (82), suggesting that far more bacterial taxa are capable of azo reduction than previously appreciated.L-Cysteine can be made from a variety of different sources, thus methionine and other sulfur starting points could impact the gut microbiota's H 2 S production capacity.Multiple Desulfovibrio species found in the GI tract can reduce sulfate to H 2 S (83-86).Thus, GI levels of H 2 S represent the net effect of multiple distinct pathways, which will be important to consider when translating these results into humans.As a first step, studies in gnotobiotic mice with defined microbial communities representing distinct assemblages of H 2 S producers would help to better understand the relative importance of each species for azo dye degradation under different dietary and host selection pressures.
There are multiple limitations of this study.We focused entirely on the model gut proteobacterium E. coli due to its genetic tractability and long history of molecular research; however, the generalizability of the pathway we described to other gut bacteria remains unclear.The development of genetic tools that permit similar mechanistic dissections into these more exotic gut bacteria will be essential for studying the conservation of enzymatic and nonenzymatic azo bond degradation in diverse gut bacterial species.The physiological relevance of our findings remains unclear, requiring more work in mice or other preclinical models.It will also be important to test the clinical relevance of these findings in the context of diseases intentionally treated with azo bond-containing drugs (e.g., inflammatory bowel disease) or diseases such as cancer that could potentially be exacerbated by dietary exposures to azo dyes.Our results empha size that in addition to metagenomics, it will be critical to measure key environmental factors in human clinical cohorts such as oxygen and amino acid levels.
Additionally, our results favor the hypothesis that FNR indirectly regulates cyuR and in turn cyuAP.While our data indicate that fnrS is responsible for the majority of this indirect regulation, these results do not address the combinatorial impact of all intermediate regulators.Ideally, a quadruple knockout of marA, soxS, fnrS, and cyuR would be necessary to confirm these are the only pathways leading toward cyuAP expression and thus H 2 S production.These results again highlight the complex regulation that underlies metabolic pathways and the importance of taking a nuanced approach to understanding not only gene product function but regulatory factors.
Despite these limitations, our results clearly emphasize the importance of environ mental factors in controlling a clinically relevant metabolic activity of the human gut microbiome.These experiments are particularly timely given the rapid increase in food dye intake due to their incorporation in processed foods (87) and emerging evidence that they have unintended consequences for gut epithelial (18) and immune (22) cells.Continued progress toward dissecting the mechanisms responsible for this metabolic activity in E. coli and other experimentally tractable gut bacterial species is a key step toward predicting and controlling the location and extent of azo dye metabolism.

Bacterial strains, media, and chemicals
Bacterial strains used in this study are listed in Table S7.E. coli was grown in lysogeny broth (LB) supplemented with 0.05% L-Cysteine under anaerobic conditions at 37°C unless otherwise noted.L-Cysteine was purchased from Sigma-Aldrich (St. Louis, MO).

Assaying whole-cell E. coli depletion of azo compounds
Cultures were grown in LB medium supplemented with 250 µM of azo dye and 0.05% L-Cysteine.Azo dye stocks were prepared in DMSO at a concentration of 25 mM and used at a final concentration of 1% (v/v).Cultures were incubated at 37°C.All incubations were performed in a COY Laboratory Products Inc anaerobic chamber (Grass Lake, MI) with the following atmosphere: 5% H 2 , 20% CO 2 , and 75% N 2 .All reagents were equilibrated in the anaerobic chamber for at least 24 h before use.After 24 h of growth, samples were removed, spun down to remove cells, and analyzed for residual dye concentration by measuring absorbance at appropriate wavelengths.Wavelengths used: 450 nm (FD&C Red No. 40, FD&C Yellow No. 6, D&C Orange No. 4), 360 nm (balsalazide, sulfasalazine, olsalazine, phenazopyridine), 420 nm (D&C Brown No. 1), 500 nm (D&C Red No. 6, FD&C Red No. 4), 535 nm (D&C Red No. 33).Residual dye concentration was determined by creating a standard curve of dye concentrations in LB and measuring absorbance at appropriate wavelengths.Concentration was calculated using the equation of the fitted standard curve line.

Clean deletion creation from Keio collection strains
Competent cells were created from cultures of the Keio collection knockouts for azoR::Kan R and fnr::Kan R knockouts.Cells were electroporated with vector pSIJ8 (88) and transformants were selected for on LB supplemented with 25 µg/mL ampicillin at 30°C.Removal of the kanamycin cassette at FRT sites was done on liquid culture of transform ant using 50 mM L-rhamnose for 4 h, then patch streaked to identify loss of kanamycin resistance.Transformants were streakpurified two times, and then PCR was performed for confirmation of kanamycin cassette removal.Strains were cured of plasmid at 37°C and then patch-streaked to confirm loss of ampicillin resistance.

Purified enzyme preparation
The ASKA collection (38) strain expressing azoR was streaked from glycerol stock on LB with 30 µg/mL chloramphenicol.A single colony was used to inoculate an overnight culture of 25 mL of LB with 30 µg/mL chloramphenicol aerobically at 37°C.The overnight was used to inoculate 1 L of LB with 30 µg/mL chloramphenicol at a 1:10 dilution factor.Sample was grown with shaking aerobically to an optical density of 0.5 at 600 nm.Expression of azoR was induced with isopropyl β-D-1-thiogalactopyranoside (IPTG) at 1 mM concentration for 3 h.Cells were then cooled on ice water for 10 min and centrifuged at 3,500 × g for 15 min, washed with ice cold PBS, spun again, then cell pellets were frozen overnight at −20°C.Pellets were thawed on ice, then resuspended in buffer A (50 mM HEPES, 300 mM NaCl, 10 mM imidazole, pH 7.5) and sonicated with the following program: 5 min, 4 s on, 4 s off, at 30% amplitude.Protease tablets (Sigma Aldrich) were crushed in buffer A and added to the sonicated lysate.Lysate was centrifuged for 20 min at 16,000 rpm.Supernatant was applied to a 5 mL nickel-nitrilotri acetic acid column (Qiagen, Valencia, CA) that had been equilibrated with buffer A at a rate of 1 mL/min.The column was washed with five column volumes of buffer A, then five column volumes of buffer B (50 mM HEPES, 300 mM NaCl, 20 mM imidazole, pH 7.5).Protein was eluted with buffer C (50 mM Hepes, 300 mM NaCl, 300 mM imidazole) and collected in 5 mL portions.Fractions were analyzed for protein content using a sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) gel.Fractions containing proteins corresponding to AzoR molecular weight were pooled and dialyzed with 3 L of buffer D (50 mM Hepes, 300 mM NaCl, 15% glycerol by weight, in cold water, pH 7.5) overnight at 4°C.Protein concentration was determined using Nanodrop (Thermo Scientific).

Purified enzyme experiments
AzoR reactions with azo dyes were carried out in 20 mM sodium phosphate buffer, 100 µM azo dye (unless otherwise noted), 2,000 µM NADH, and 20 µM FMN with varying concentrations of enzyme as appropriate.Reactions were activated by the addition of NADH.Dye concentration was monitored at absorbance 450 nm every 30 s after a linear shake of 15 s.All reactions were performed aerobically at 37°C.

Keio collection screening for alternative azo reducers
Homologs to genes implicated in azoreduction by S. oneidensis MR-1 (40) were identified by first gene name and second gene function (tested or predicted) in E. coli.Other strains of interest were identified based on gene name and function.All strains were cross-referenced on EcoCyc (41).Strains were revived from glycerol stocks on LB agar with 25 µg/mL kanamycin.Single isolates were picked from agar plates and grown overnight in 1 mL of LB supplemented with 0.05% L-Cysteine in a 96-deep well plate under anaerobic conditions.For azoreduction assays, previously anaerobically equilibra ted liquid media was supplemented with 250 µM FD&C Red No. 40 and 0.05% L-Cysteine. 1 mL of media was aliquoted to 2 mL deep 96-well plates and strains were inoculated at 1:100 from overnight cultures in triplicate.Plates were sealed with TempPlate Sealing Foil (USA Scientific, Cat #2923-0110) and allowed to grow for 24 h.After 24 h, cells were spun down and 100 µL of supernatant was used to measure dye absorbance at 450 nm.Proportion of dye remaining was calculated by dividing the absorbance of each well by the average absorbance of uninoculated control wells.Where standard curves were used, a 1:2 dilution series of FD&C Red No. 40 was made in LB with 0.05% L-Cysteine in triplicate and dye absorbance measured at 450 nm.Concentrations were calculated using the equation of the standard curve line in Excel.

Validation of Keio collection knockouts
To ensure the knockouts used from the Keio collection were indeed the annotated gene, primers were designed up-and down-stream of the gene of interest.PCR and Sanger sequencing of the kanamycin insert and surrounding regions were performed to ensure the expected band size and gene alignment.Primers are listed in Table S8.

Complement strain construction and assay
Complement strains were made by preparing competent cells for the background strain of interest.Plasmids were prepared by growing up overnight cultures from the ASKA strain of interest, then using Qiagen Plasmid Mini Kit (Qiagen; Cat #12125) to extract the pCA24N plasmid with the gene of interest.The plasmid preparation was desalted using a membrane filter on water, then electroporated to the background strain of interest.Correct insertion was ensured by growth of the transformant on chloramphenicol.A mini prep was also made from the transformant, and PCR and Sanger sequencing done to ensure the expected gene insertion.The same procedure was used to introduce the pBD-pLac and pBD-pLac-fnrS into strains of interest.

RNA sequencing and proteomics sample growth and exposure
100 mL of overnight cultures of both wt and Δfnr E. coli were grown in LB supplemented with 0.05% L-Cysteine under anaerobic conditions.Overnight culture OD 600 nm was measured for both cultures and normalized to a starting OD 600 nm of 0.05 in 100 mL of LB with 0.05% L-Cysteine in 12 replicates for each strain.OD 600 nm was monitored until mid-exponential phase for each strain (previously determined).At mid-exponential three replicates for each strain were treated with a final concentration of 250 µM FD&C Red No. 40 in DMSO, or 1% DMSO.After 45 min of exposure, 15 mL of the sample was removed for RNA extraction, while the remaining volume was used for proteomics.Both volumes of samples were spun down to pellet cells, supernatant removed, and cell pellets were flash frozen in liquid nitrogen then stored at −80°C.The remaining samples were allowed to reach stationary phase, then they were treated with a final concentration of 250 µM FD&C Red No. 40 in DMSO, or 1% DMSO.Samples were processed the same way as for mid-exponential.

RNA sequencing sample preparation and analysis
1 µL of TRI Reagent (Sigma Aldrich catalog number T9424) was added to bacterial pellets and incubated for 10 min, samples were then transferred to 2 mL Lysing Matrix E tubes (MP Biomecials, catalog number 116914050).Cells were lysed for 5 min in the bead beater at room temperature.200 µL of chloroform was added.Samples were vortexed for 15 s and incubated at room temperature for 10 min.Samples were then centrifuged at 16,000 × g for 15 min at 4°C. 500 µL of the upper aqueous phase was transferred to a new tube, 500 µL of 100% ethanol was added, and vortexed to mix.Mixture was transferred to a spin column (PureLink RNA Mini Kit; Life Technologies; catalog number 12183025) and centrifuged at ≥12,000 × g for 30 s, discarding flowthrough, until all the material had been added to the column.To the spin column, 350 µL Wash Buffer I (PureLink RNA Mini Kit; Life Technologies; catalog number 12183025) was added, then the column was centrifuged at ≥12,000 × g for 30 s, discarding flowthrough.80 µL of PureLink DNase mix was added to the column and incubated at room temperature for 15 min.350 µL of Wash Buffer I (Purelink RNA mini kit) was added, and the column spun at >12,000 × g for 30 s.The column was transferred to a new collection tube, and 500 µL Wash Buffer II was added, followed by centrifugation at ≥12,000 × g for 30 s, discarding flowthrough.The column was centrifuged at ≥12,000 × g for 60 s and dried, then moved to a collection tube.50 µL RNase-free water was added, and the column was incubated at room temperature for 1 min.Finally, the column was centrifuged for 1 min at ≥12,000 × g, retaining the flowthrough, which contained total RNA.
Samples were DNase treated again using TURBO-DNase (Ambion; ThermoFisher catalog number AM2238) and incubated at 37°C for 30 min.RNA Ampure XP Beads were used to clean up the reactions.1.8 volumes of the RNA Ampure XP beads were added to 1 vol of RNA sample and allowed to sit for 5 min at room temperature.Tubes were placed on a magnetic stand until the liquid cleared.Then liquid was removed, and beads washed with 200 µL of 100% ethanol.Samples were incubated for 30 s, then removed.This was repeated twice.Samples were allowed to dry for 5 min, removed from the magnetic stand, and 30 µL of RNase-free water was added.This was incubated for 2 min at room temperature then placed back on the magnetic rack where liquid was collected after turning clear.
rRNA was removed from total RNA using Ribominus Transcriptome Isolation Kit for Bacteria and Yeast (Invitrogen; catalog number K155004, LOT: 2116711), following manufacturer's protocol.RNA fragmentation, cDNA synthesis, and library preparation were performed using the NEBNext Ultra RNA library Prep Kit for Illumina and NEBNext Multiplex Oligos for Illumina (Dual Index Primers) (Ipswich, MA).Samples were dual-end sequenced (2 × 75 bp) using the NextSeq Mid Output platform (Table S2).Reads were mapped to the E. coli K12 BW25113 genome sequence (NCBI Reference Sequence: GCA_000750555.1) using Bowtie2 (89), and HTSeq (90) was used to count the number of reads to E. coli genes.Differential gene expression was analyzed using limma (91).Differentially expressed genes were defined as transcripts exhibiting an absolute log 2 fold change ≥1 and p adj <0.1.Detection of expression was based on limma's default filtering parameters in the "filter genes by expression level" function.

Proteomics sample preparation
Cell pellets were washed three times with 5 mL of cold PBS and resuspended in 1 mL, transferred to 2 mL Lysing Matrix E tubes (MP Biomecials, catalog number 116914050) and lysed via bead beating for 1 min, then resting tubes on ice for 2 min.This was repeated twice.Samples were centrifuged for 15 min at 4°C at 16,000 × g.A BCA assay was performed to quantify protein concentration in the lysate supernatant.Urea was added to the samples to a target concentration of 8 M followed by an appropriate volume of dithiothreitol to obtain a 5 mM concentration.Samples were then incubated at 60°C for 30 min.While samples were incubating, the trypsin was preactivated for 10 min at 37°C.Samples were diluted 10-fold with 100 mM NH 4 HCO 3 .1 M CaCl 2 was added at an appropriate volume to create a final sample concentration of 1 mM CaCl 2 .Samples were then digested for 3 h with trypsin at 37°C at a concentration of 1 µg trypsin/50 µg protein.Samples were snap frozen.
Samples were cleaned with a C 18 column on a vacuum manifold.The column was conditioned with 3 mL of methanol, then rinsed with 2 mL of 0.1% TFA acidified water.Sample was run though column, then column was rinsed with 4 mL of 95:5 H 2 O:acetoni trile (ACN), 0.1% trifluoroacetic acid (TFA).The column was allowed to go to dryness, then eluted slowly to dryness with 1 mL of 80:20 ACN:H 2 O, 0.1% TFA into a collection tube.Samples were concentrated in the speed-vac to a volume of approximately 50-100 μL.Protein quantification was done using a BCA test, then samples were stored at −80°C until analysis.

Proteomics
A Waters nano-Acquity M-Class dual pumping UPLC system (Milford, MA) was configured for online trapping of a 5 µL injection at 5 µL/min with reverseflow elution into the analytical column at 300 nL/min.The trapping column was packed in-house (PNNL) using 360 µm o.d.fused silica (Polymicro Technologies Inc., Phoenix, AZ) with 5 mm Kasil frits for media retention and contained Jupiter C 18 media (Phenomenex, Torrence, CA) in 5 µm particle size for the trapping column (150 µm i.d.× 4 cm long) and 3 µm particle size for the analytical column (75 µm i.d.× 70 cm long).Mobile phases consisted of (A) water with 0.1% formic acid and (B) acetonitrile with 0.1% formic acid.The following gradient profile was performed (min, %B): 0, 1; 2, 8; 25, 12; 85, 35; 105, 55; 110, 95; 115, 95; 117, 50; 119, 95; 121, 95; 123, 1. MS analysis was performed using a Velos Orbitrap Elite mass spectrometer (Thermo Scientific, San Jose, CA) outfitted with an in-house made nano-electrospray ionization interface.Electrospray emitters were prepared using 150 µm o.d.× 20 µm i.d.chemically etched fused silica (92).The ion transfer tube temperature and spray voltage were 325°C and 2.2 kV, respectively.Data were collected for 100 min following a 20 min delay from sample injection.FT-MS spectra were acquired from 400 to 2000 m/z at a resolution of 35 k (AGC target 3e6) and while the top 12 FT-HCD-MS/MS spectra were acquired in data-dependent mode with an isolation window of 2.0 m/z and at a resolution of 17.5 k (AGC target 1e5) using a normalized collision energy of 30 s and a 30 s exclusion time.Extracted MS data were run through the MS-GF+ database and filtered by a FDR of 0.01 using the MS-GF+ generated decoy database.Differential peptide abundance was analyzed using DEP (93) and limma (91) using a cutoff of an absolute log 2 fold change ≥1 and p adj < 0.1.

Regulator enrichment analysis
Features found significantly differentially expressed (log 2 fold change > |1|, p adj < 0.1) in both the proteomics and transcriptomics data sets were imported to EcoCyc (41) for enrichement analysis.The gene list was analyzed using the "Genes enriched for transcriptional/translational regulators (direct only)".The following settings were used: p adj < 0.1, Fisher Exact Statistic, Benjamini-Hochberg Correction).Proportion of regulon represented was determined using the list of genes "directly regulated by gene" feature of EcoCyc for total regulon size.Then, the number of genes regulated by each regulator was divided by number of total genes regulated by that particular regulator.

Hydrogen sulfide quantification
H 2 S was measured using the Cline reaction (94) with modifications.Samples were put into a zinc acetate buffer (16.7 mM) on ice at a ratio of 1:3 sample:buffer.Then 180 µL of zinc acetate solution was transferred and mixed with 20 µL of line reagent (2 g of N,N-dimethyl-p-phenylenediamine sulfate with 3 g of FeCl 3 in 500 mL of cold 50% v/v HCl).The reaction was allowed to proceed for 20-30 min at room temperature in the dark.Then absorbances were measured at 670 nm.A standard curve of sodium sulfide was created in zinc acetate and reacted with Cline reagents for sulfide concentration calculations.

RT-qPCR assays
Samples were grown under anaerobic conditions with 250 µM FD&C Red No. 40.At mid-exponential phase, 4.2 mM L-Cysteine was added to cultures for 15 min.Then samples were removed from anaerobic conditions, centrifuged for 5 min at ≥12,000 × g, supernatant removed, and pellets were frozen at −80°C.Samples were thawed and RNA extracted as previously described for RNA sequencing.RT-qPCR was done using RioRad's 1-step Quantitative Reverse Transcription PCR kit (Bio-Rad, Hercules, CA, catalog number 1725150) according to kit protocol.
To analyze the differences between wt and Δfnr samples, fold change of sample ΔC t minus the maximum ΔC t for each gene across all samples in wt or Δfnr E. coli was calculated.Within each replicate the Ct of the housekeeping gene was subtracted from the C t of the gene of interest (ΔC t ).The maximum ΔC t across samples from both strains was subtracted from all samples to give ΔΔC t .Fold change was calculated as 2^(-ΔΔC t ).To analyze the impact of L-Cysteine on gene expression, fold change was calculated using expression of the sample with L-Cysteine compared to the vehicle.Within each replicate the C t of the housekeeping gene was subtracted from the C t of the gene of interest (ΔC t ).The ΔC t from the paired replicate control (vehicle exposure) samples was subtracted from the exposed sample ΔC t to give ΔΔC t .Fold change was calculated as 2^(-ΔΔC t ).

Double gene knockout creation
P1 lysates were generated of each strain of interest carrying the kanamycin resistance cassette (Keio collection) adapting methods from previously described techniques (95).Briefly, 150 µL of overnight culture in LB supplemented with 12.5 µg/mL kanamycin was mixed with 1 to 25 µL of P1 phage (previously propagated from ATCC on MG1655).This mixture was incubated for 10 min at 37°C to aid absorption, added to 3 mL of 0.7% agar, and overlaid on prewarmed LB agar supplemented with 25 µg/mL kanamycin and 10 mM MgSO 4 .Plates were incubated overnight at 37°C, and phage was harvested by adding 5 mL of SM buffer, incubating at room temperature for 10 min, and breaking the top agar for phage harvest.The mixture was briefly centrifuged to pellet agar, then supernatant was passed through a 100 µm cell straining, then 0.45 µm syringe filter.Lysates were stored at 4°C.
A clean deletion of the recipient strain of interest was created using the earlier described method with pSIJ8 (88).To transduce the clean recipient strain, 1 mL of an overnight culture of recipient strain was pelleted and resuspended in ⅓ volume of LB 10 mM MgSO 4 , 5 mM CaCl 2 .100 µL of cells were mixed with 1 μL to 10 µL of P1 lysate and incubated for 60 min at 37°C.200 µL of 1M sodium citrate was added with 1 mL of LB to minimize secondary infection.This mixture was incubated at 37°C for 2 h, then plated on LB supplemented with 10 mM sodium citrate and 25 µg/mL kanamycin to select for transductants.Transduction was confirmed using PCR and Sanger sequencing.Azoreduction and sulfide activity were tested as described earlier.

Statistical analysis
For multiple comparisons of data, an ordinary one-way ANOVA test with Dunnet's correction was chosen, unless otherwise specified.For the targeted knockout screen, a 95% CI of averages from two experiments was calculated using GraphPad Prism v9.Michaelis-Menten curves were created in GraphPad Prism v9.For proteomics and transcriptomics analysis, R package limma was used to determine the fold change and adjusted P values.Feature count differences and Mann-Whitney significance testing were done in R. For the regulator enrichment analysis, EcoCyc's enrichment analysis for "Genes enriched for transcriptional/translational regulators (direct only)" was used with a Fisher's exact test with Benjamini-Hochberg correction.

FFIG 1
FIG 1 E. coli AzoR is sufficient but dispensable for azoreductase activity.(A) Azoreduction activity is significantly greater under anaerobic conditions.(B) Significant differences in FD&C Red No. 40 dye depletion is observed when E. coli is grown with or without L-Cysteine in the media.E. coli azoR is dispensable for azoreductase activity during anaerobic growth on (C) solid and (D) liquid media.(E) The knockout phenotype holds with a clean azoR deletion.(A through F and H) 250 µM FD&C Red No. 40.(C) 72 h growth on LB agar with 4.13 mM L-Cysteine.(D and E) 24 h growth in LB with 4.13 mM L-Cysteine.(A and B, D and E) Supernatant was removed from samples and analyzed for residual dye concentration spectrophotometrically.Concentrations calculated based on a standard curve, limit of detection = 1.2 µM.(F) Azoreduction activity over a 24 h time course.Media with 4.13 mM L-Cysteine and 250 µM FD&C Red No. 40 were inoculated with wt or ΔazoR and dye depletion was monitored over time.(G) Michaelis-Menten curve for purified AzoR enzyme with FD&C Red No. 40.Michaelis-Menten curve fit using PRISM.(H) Dye depletion in liquid media by wt or ΔazoR across a panel of azo drugs and dyes.Grown with 250 µM azo dye/drug, for 24 h in LB with 4.13 mM L-Cysteine.Fold changes are relative to average blank concentration.(A and B, D and E) Bars are mean.(F) Lines are mean ± stdev.(A and B) *P < 0.05, ****P < 0.0001, (A and F) two-way ANOVA, (B) Mann-Whitney test, (D and E) one-way ANOVA.(A and B, D through H) n = 3 biological replicates/strain.

FIG 2 FFIG 3
FIG 2The fnr regulator is necessary for azoreductase activity in E. coli.(A) A targeted screen of 113 deletions of genes that encode enzymes and regulators revealed 22 loss-of-function strains (red bars, more dye remaining than the 95% confidence interval of the mean; n = 6 biological replicates/strain across two independent experiments).The Δfnr (fumarate and nitrate reduction regulator) strain had the most extreme phenotype (complete loss-of-function).No significant gain-of-function phenotypes were observed.(B) Removal of the kanamycin cassette from Δfnr still leads to a complete loss-of-function (n = 4 biological replicates/strain).(C) Complementation of fnr but not azoR rescues the azoreduction phenotype.Samples all induced with 0.1 mM IPTG at mid-exponential growth phase (n = 3 biological replicates/strain).All experiments: 24 h of growth in LB media, with 4.13 mM L-Cysteine and 250 µM FD&C Red No. 40, under anaerobic conditions.(B and C) Values are mean ± stdev.****P < 0.0001, one-way ANOVA.Supernatant was removed from samples and analyzed for residual dye concentration spectrophotometrically.Concentrations calculated based on a standard curve.Limit of detection = 1.2 µM.

FIG 3 (
FIG 3 (Continued) FD&C Red No. 40 or vehicle control.After 45 min, cultures were spun down and fractions taken for proteomics or transcrip tomics processing.Created with BioRender.com.(B and C) Principal components analysis (PCA) of Euclidean distances for transcriptomics (B) and proteomics (C) datasets.Each point represents one sample from the experiment.(D and E) Bar plots representing the number of features enriched or depleted (p adj <0.1, log 2 fold change>|1|) for each comparison for transcriptomics (D) and proteomics (E) datasets.(F and G) Volcano plots for transcriptomics (F) and proteomics (G) datasets.Each point represents the average expression of one transcript or protein, respectively.Differential expression is based on the strain-level comparison.Significantly differentially expressed features (p adj <0.1, log 2 fold change≥|1|) are indicated by dashed lines, comparing Δfnr to wt. (H) Venn diagram of the number of differentially expressed features (p adj <0.1, log 2 fold change >|1|) in the two datasets.(I) Regulons enriched in the overlapping features of transcriptomics and proteomics datasets (p adj <0.1, with Benjamini-Hochberg correction) of differentially expressed genes (p adj <0.1, log 2 fold change>|1|).(J) Azo depletion activity by enriched regulator knockout strains from the Keio Collection.Values are mean ± stdev.

FIG 4
FIG 4 The fnr regulator is necessary for the metabolism of L-Cysteine to hydrogen sulfide, enabling the reduction of azo bonds.(A) Schematic depicting L-Cysteine metabolism pathway in E. coli, which metabolizes L-Cysteine to pyruvate, releasing hydrogen sulfide.Created with BioRender.com.(B) wt and Δfnr E. coli were incubated in various concentrations of L-Cysteine.Statistics relative to blanks.(C) Increasing L-Cysteine concentrations in LB with 250 µM FD&C Red No. 40.(D) Increasing sodium sulfide concentrations in LB with 250 µM FD&C Red No. 40 with 0.05% L-Cysteine.(E) Hydrogen sulfide production by wt and Δfnr E. coli.(F) FD&C Red No. 40 depletion by wt and Δfnr E. coli.(G) FD&C Red No. 40 depletion and sulfide measurements over time in wt (blue/green) and Δfnr (red/orange) E. coli.(B, E and F) *P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001, one-way ANOVA.All experiments: values are mean ± stdev.24 h of growth in LB media, with 4.13 mM L-Cysteine and 250 µM FD&C Red No. 40 unless otherwise noted.Supernatant was removed from samples and analyzed for residual dye concentration spectrophotometrically.Concentrations calculated based on a standard curve.Limit of detection = 1.2 µM (n = 3 biological replicates).

FIG 5 L
FIG 5 L-Cysteine import and metabolism genes impact bacterially produced sulfide and FD&C Red No. 40 depletion.(A) Model of the cyuR regulon's role in L-Cysteine metabolism.CyuP acts as a L-Cysteine transporter into the cell, while CyuA is a L-Cysteine desulfidase.Created with BioRender.com.(B) Volcano plot for transcriptomics data: cyuA and cyuP are highlighted in pink.Each point represents the average expression of one transcript.Differential expression is based on the strain-level comparison.(C and D) Paired sulfide (C) and FD&C Red No. 40 (D) were measured from cyuR regulon knockouts and cyuR (Keio collection).(E and F) Paired sulfide (E) and FD&C Red No. 40 (F) were measured from ΔazoR, ΔcyuR::Kan R , and ΔazoRΔcyuR::Kan R strains.(C through F) Values are mean ± stdev.*P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001, one-way ANOVA (n = 3 replicates/strain).

FIG 6
FIG 6 FNR indirectly regulates sulfide production and FD&C Red No. 40 depletion predominantly via small RNA fnrS.(A) Model of the potential regulation cascade from FNR to the expression of cyuAP.Created with BioRender.com.(B and C) Paired FD&C Red No. 40 (B) and sulfide (C) were measured from wt, Δfnr::Kan R complemented with pBD-lac or pBD-lac-fnrS, and ΔfnrS::Kan R complemented with pBD-lac or pBD-lac-fnrS.(D and E) Paired FD&C Red No. 40 (D) and sulfide (E) were measured from wt, Δfnr::Kan R , and ΔfnrS::Kan R over time.Samples inoculated at time point zero.(E) Plotted on a log 10 scale y-axis.(F) Fold change of fnrS in the presence of L-Cysteine compared to vehicle at mid-exponential growth.All values are mean ± stdev.*P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001.(B and C) One-way ANOVA (n = 4 replicates/strain).(D and E) Two-way ANOVA for overall comparisons and two-way ANOVA comparing wt to each group at the final time point (n = 3 replicates/strain).(F) Mann-Whitney test (n = 4 biological replicates averaged from three technical replicates/strain).

FIG 7
FIG 7 Working model describing how the host environment shapes gut bacterial azo reduction.Under anaerobic conditions, FNR activates the transcription of fnrS, a small noncoding regulatory RNA.fnrS in turn activates expression of the cyuAP operon.Alternatively, fnrS may upregulate the expression of marA which can drive expression of both cyuR and soxS.SoxS can also induce cyuR expression.CyuP transports L-Cysteine into the cell, where it is metabolized to H 2 S through the CyuA desulfidase.Released H 2 S can reduce the azo bond of azo compounds.Created with BioRender.com.